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This  is  a  triple-exposure  photograph  showing  a  10X  microscope 
objective  (to  the  left),  a  gallium  arsenide  rectangular  substrate  with  a 
semi-circular  film  of  zinc  sulfide,  and  a  penny  for  size  comparison. 
Also  shown  is  a  red  (A  =  6328  A)  laser  beam  in  two  positions.   The  upper 
beam  (in  air)  is  visible  due  to  scatter  from  surface  dust  particles  on 
the  ZnS-GaAs  sample  as  well  as  from  smoke  used  to  help  indicate  the 
path  of  the  beam.   This  beam  (slightly  off-axis  with  respect  to  the  10X 
objective  axis)  represents  light  rays  undeviated  by  the  film.   The  lower 
beam  represents  light  refracted  away  from  the  undeviated  beam  according 
to  Snell's  Law  into  the  zinc  sulfide  when  the  film  is  placed  at  the  focal 
point  of  the  objective.   The  beam  is  guided  by  the  film  and  propagates 
until  it  reaches  the  other  end  of  the  semi-circle,  at  which  point  the 
light  radiates  brightly.   That  this  lower  beam  is  visible  may  be 
attributed  to  surface  irregularities  as  well  as  to  internal  scattering 
from  ZnS  crystal  imperfections. 


To  my  wife,  Ingrid 
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PREFACE 


On  June  26,  1972  Dr.  Stanley  S.  Ballard,  Professor  of  Physics, 
University  of  Florida,  visited  the  Optical  Systems  Division  of  Itek 
Corporation  and  met  with  this  author  and  with  Dr.  Julius  Feinleib, 
Manager  of  the  Physics  Scientific  Staff  (part  of  Itek's  Central 
Research  Laboratories).   Dr.  Feinleib  described  Itek's  capability  to 
fabricate  cubic-structure  zinc  sulfide  films  on  gallium  arsenide 
substrates.   As  a  result  of  that  meeting,  it  was  agreed  that  ZnS-on- 
GaAs  samples  would  be  made  available  for  electro-optical  waveguide 
experiments . 

Subsequently,  preliminary  experiments  and  analyses  were  performed, 
and  a  proposed  line  of  research  was  submitted  to  the  author's  Ph.D. 
supervisory  committee.   That  proposal  suggested  both  a  theoretical 
and  an  experimental  research  effort.   The  theoretical  part  was  to 
concentrate  on  the  propagation  of  light  in  inhomogeneous  waveguides. 
The  experimental  part  was  to  examine  the  electro-optic  effect  of  zinc 
sulfide  waveguides.   This  dissertation  presents  the  results  of  research 
encompassing  both  aspects  of  the  proposed  research. 
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The  basic  physics  of  homogeneous  asymmetric  waveguides  is  reviewed, 
allowing  ready  comparison  of  similar  concepts  between  homogeneous  and 
inhomogeneous  optical  waveguides.   It  is  shown  that  a  quasi-geometrical 
ray-optics  approach  leads  to  the  same  results  obtained  from  more  rigorous 
mathematics,  permitting  a  practical  understanding  of  the  phenomenon 
of  total  internal  refraction  and/or  reflection  in  general-index-gradient 
planar  optical  waveguides.   The  nhase  change  of  tt/2  (not  0  or  tt)  occurring 
upon  total  internal  refraction  is  explained  as  being  due  to  a  cylindrical 
focussing  of  the  light  within  the  guide.   Vaveguiding  and  electro-ODtic 
modulation  of  light  in  epitaxial  zinc  sulfide  films  on  gallium  arsenide 
substrates  are  also  discussed.   Two  mechanisms  for  light  guiding  at  the 
ZnS/GaAs  boundary  are  described:   total  internal  refraction  or  partial 
internal  reflection.  From  the  discovery  of  a  very  thin,  transparent, 
low-index  inhomogeneous  layer  of  zinc  sulfide  near  the  boundary,  it  is 
concluded  that  1)  90-95%  of  the  zinc  sulfide  thin  film  is  homogeneous, 
and  2)  due  to  this  thin  inhomogenous  laver,  the  zinc  sulfide  guide  is 
isolated  from  the  gallium  arsenide  substrate;  therefore  the  guide  is 
not  a  leaky  one  but  instead  is  a  true  dielectric  waveguide.   Finally, 
electro-optic  measurements  are  presented  which  agree  with  those  pre- 
dicted from  theory. 


CHAPTER  I 
INTRODUCTION 

This  dissertation  begins  with  a  discussion  of  homogeneous  waveguides 
(even  though  they  have  been  treated  by  other  authors)  and  develops  two 
systematic  means  for  understanding  the  phenomenon  of  total  internal 
refraction  and/or  reflection  in  inhomogeneous  waveguides;  a  quasi- 
geometrical  ray-optics  approach,  and  a  more  rigorous  approach  using 
electric  and  magnetic  vector  potentials  to  solve  Maxwell's  equations. 
Having  done  this,  the  theoretical  and  physical-insight  discussions  of  the 
inhomogeneous  case  should  present  no  problems,  for  although  the  mathe- 
matics are  more  complicated  and  more  tedious,  the  approach  is  basically 
the  same  as  in  the  homogeneous  case. 

Because  of  the  complexity  of  the  mathematics  in  the  theoretical 
part,  the  long  derivations  have  been  relegated  to  several  appendices. 
Then,  in  the  main  body,  summaries  of  the  important  results  are 
presented  with  physically  based  explanations,  permitting  a  practical 
understanding  of  the  propagation  of  light  in  general-index-gradient 
planar  optical  waveguides. 

After  homogeneous  and  inhomogeneous  waveguides  are  discussed,  the 
results  of  waveguiding  and  electro-optic  light  modulation  experiments 
with  epitaxial  zinc  sulfide  films  are  presented. 


CHAPTER  II 
BASIC  PHYSICS  OF  HOMOGENEOUS  ASYMMETRIC  PLANAR  OPTICAL  WAVEGUIDES 

Introduction 
A  planar  or  slab  (as  opposed  to  a  circular  fiber)  optical  wave- 
guide is  one  that  confines  light  in  one  dimension  only,  as  shown  in 
Figure  II-l.   Light  propagating  in  the  +  z  direction  is  guided  or 
confined  by  the  two  surfaces  at  x  =  0,  -W  by  total  internal  reflec- 
tion.  There  is  no  confinement  in  the  y  direction.   A  homogeneous 
optical  waveguide  is  one  that  has  an  index  of  refraction  n  which  is 
constant  throughout  the  guide.   An  asymmetric  guide  has  a  substrate 

refractive  index  n  which  is  not  equal  to  the  superstrate  index  n  , 
s  a 

i.e.,  n  ^  n  .   With  no  loss  of  generality,  n   is  assumed  to  be 
s    a  s 

greater  than  n  ,  i.e.,  n   >  n  .   A  large  number  of  optical  waveguides 

a         s    a  ° 

can  be  classified  as  homogeneous  asymmetric  planar  waveguides. 
Vector  Potentials 

In  order  to  describe  the  propagation  of  light  in  optical 

1  2 
waveguides,  the  method  of  dual  vector  potentials    is  reviewed  in 

Appendix  A.   Beginning  with  Maxwell's  equations,  the  wave  equations 

for  E  and  F  are  obtained  and  the  form  of  solution  is  shown.   Then, 

for  the  case  where  no  charge  density  or  currents  are  present,  it  is 

shown  that  both  E  and  H  can  be  expressed  in  terms  of  a  magnetic  vector 

potential  A  and  an  electric  vector  potential  F.   The  wave  equations 

for  A  and  F  for  both  the  homogeneous  and  inhomogeneous  cases  are 


x=0 


x=-W 


*s^ 


6— ** 


Figure  II-l.   Cross-sectional  view  of  light  propagating  in  a 
homogeneous  asymmetric  planar  optical  waveguide, 


v-F  +  A,"p-o  u_(2) 

2 
Of  course,  for  the  homogeneous  case,  Vn  =  0  in  Eq.  (1).   Here 

k     =   2it/A      =  co/c,    A      is   the   free-space  wavelength,    co   is    the 
ooo  &      > 

angular  frequency,  and  c  is  the  speed  of  light  in  free  space.   The 
following  expressions  result  for  E  and  H: 


=   VxA    ~  (C/*.)    V(7>P)  -  <>o0aA 


ti-(3) 


u-(4) 


For  the  choice  of  potential 


II-(5) 


II- (6) 
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(where  x,  v,  and  z  are  Cartesian  unit  normals),  it  is  shown  that 

x  H-(7) 

a  **                                                   n-(8) 
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H  =  ~  ^/A^iJfi.  ii- (12) 


For  the  choice 


A-0 

II-(13) 


II-(14) 


then  E  =  0,  F  =0,  and  H  =  0,  in  which  case,  because  there  is  no 

x       z  y 

z-component  of  E,  a  pure  transverse  electric  (TE)  mode  of  propagation 
exists.   For  the  choice 


-J      A   A 

II-C15) 


^=0  II-(16) 


then  E  =  0,  H  =0,H  =0,  and  a  pure  transverse  magnetic  (TM) 
y       x       z 

mode  of  propagation  exists.   It  mav  be  seen  that  a  field  of  any 
polarization  in  an  asymmetric  planar  optical  waveguide  can  be  expressed 
as  the  sum  (i.e.,  the  superposition)  of  a  TE  field  and  a  TM  field. 

TE  Modes 
In  Appendix  B,  homogeneous  asymmetric  waveguides  are  discussed  in 
in  detail  for  both  TE  and  TM  modes.   After  an  analysis  involving  boun- 
dary conditions,  simultaneous  equations,  and  complex  algebra,  the 


following  solutions  for  E  and  H  are  obtained  for  TE  modes: 

A  E. 


V    a  fc3 


it_=   x  Vzl, 


j    '  J  '      °y7  n-(i9) 


Cu>t -f2 z) 


p-*.»v 


II-C17) 


Te  II- (18) 


II-(20) 


Kix,Ztt)=    Hfl200e  !  ii- (21) 


Here  (3  is  the  wavevector  component  along  the  z-axis,  i.e., 


II-(22) 


and  has  a  constant  value  in  all  three  of  the  n  ,  n  ,  and  n  regions. 

a   g       s 

Eq .  (22)  implies  that  the  wave  travels  at  the  same  velocity  v  in  all 

three  regions  since  v  =  eg.   Finally, 
z 

-*,(P«H         .   rJTT  N  -OJC 

01  J  II- (23) 


1     l  °  II- (24) 

r  /     n-(25) 


and,  from  Eqs.  (8),  (10),  and  (12),  H  =  (-  i/ko)3E  /3z  and 

H  =  (i/k  3E  /3x).  The  exact  form  of  H   and  H   is  given  in 
z       o  y  ox      oz 

Appendix  B,  Eqs.  B-(83)  through  B-(86) .   These  values  for  E   ,  H   , 
■  '  >  o  oy   ox 

and  H   show  that  the  wave  is  confined  to  the  n    region  (the  wave- 
oz  g 

guide)  and  that  the  fields  go  to  zero  as  x  ->  ±  °°.   The  coefficient  F„ 
is  arbitrary,  being  determined,  for  instance,  by  how  much  light  is 
coupled  into  the  guide.   The  remaining  parameters  as  determined  in 
Appendix  B  are  given  below: 


ay/a 


«--  M0  (n3as^q9—n^) 
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II- (30) 


In  the  course  of  obtaining  the  solutions  for  E   and  H   ,  two 
useful  TE  equations  were  obtained: 

=     oLmTY 

II-C31) 


W, 


rnt/i  /1*F 


S  K^JJi 


WIT  "rtaw 


'3 


II-(32) 


where  (W  .  )   is  the  minimum  waveguide  thickness  required  in  order  to 

mm  TE 

propagate  a  mode  of  order  m.   Eq .  (31)  is  called  the  characteristic 

waveguide  equation. 

Eqs.  (31)  and  (32)  show  that,  depending  on  A  ,  n  ,  n  ,  and  n  , 

o   a   g       s 

only  certain  modes  can  propagate  in  a  film  of  thickness  W.   Further, 

for  n  /  n  ,  there  is  a  cutoff  thickness,  W  .  ,  below  which  no  guided 
s    a  mm 

wave  can  be  propagated,  even  for  m  =  0.   It  may  be  noted  that  Eq .  (31) 
is  a  transcendental  equation  in  0. 

Fq.  (32)  is  derived  from  Eq.  (31)  using  the  following  argument. 
In  order  for  the  waves  to  leak  out  or  escape  (i.e.,  for  propagation 
to  cease) ,  the  waves  must  escape  either  into  the  n  -  region  or  into 

the  n  -  region.   For  such  a  situation,  at  least  one  of  the  functions 

s 

e     [for  x  =  0]  or  e    [for  x  =  -  W]  in  Eqs.  (23)  and  (25)  must  no 

longer  be  decaying  exponentials   (which  are  real  functions)   but  must 

instead  be  purely  imaginary.   From  Eqs.  (26)  and  (28),  it  may  be  seen 

that  when  n  sinQ  <  n  or  when  n  sinG  <  n  ,  then  e    or  e 
g        a  g        s 

respectively  is  imaginary.   Thus  when  either 

Sua  9   ~  n.  lr\a 

'     3  II- (33) 


^UA  9-=  f).  IV\ 


Sin^  II-(34) 


then  waves  will  begin  to  leak  out.   When  sinQ  =  n  /n  ,  then  Q  is  at 
B  s   g 

the  critical  angle  for  total  internal  reflection  at  the  n  -  n 
&  g    s 

boundary.   Since  n   >  n  ,  total  reflection  still  occurs  at  the  n   -  n 
s    a'  a    g 

boundary.   When  sinQ  <  n  /n  ,  plane  waves  escape  into  the  n  medium, 

s   g'  P  s 


and  propagation  ceases,  i.e.,  the  mode  is  cut  off.   Using  Eq .  (34) 
in  Eq.  (31)  results  in  Eq .  (32). 

Referring  to  Eq .  (31),  it  may  be  seen  that  0  can  never  reach  the 
value  of  tt/2  (propagation  straight  down  the  guide),  since  in  the 
limit  as  0  -*■   tt/2,  W  -*■  °°,  which  is  physically  unrealistic.   Put 
another  way,  no  longer  is  there  really  any  waveguiding  as  0  -»■  tt/2. 
Combining  this  fact  with  Eq .  (34)  results  in  the  relationship 


S3  S  II-(35) 


as  the  range  of  values  0  for  which  waveguiding  (without  cutoff)  can 
occur,  as  determined  by  Eq .  (31). 

TM  Modes 
Appendix  B  also  contains  a  derivation  of  E  and  H  for  TM  modes 
propagating  in  a  homogeneous  asymmetric  waveguide.   The  solutions  are 


7*1        * 


II-(36) 


m       U      *  II-C37) 


+f-=  H  ^ 


y.;^=  Ho „(x)eU   I  ' 

!i  J  II-(40) 


II-(38) 


II-C39) 
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Here  $  has  the  value  given  in  Eq .  (22)  and  is  a  constant.   The  x 
dependence  is  given  by 

~^^e'L    *-3    doo    (  Qo^Je**       W    x±0       II- (41) 


M*^    {    "^a^      ^9**    C)      4r-W±X± 


J  II- (42) 


i€ 


-a^e'    ^(C+m )e^ 


4,^   *=-i\/      n-(43) 


From  Eqs.  (7),  (9),  and  (11),  E  =  (i/k  n2)3H  /9z  and  E  = 

x       o     y         z 

2 
(-i/k  n  )9H  /9x.   The  exact  form  of  E   and  E   is  given  in  Eqs.  B-(157) 
o     y  ox      oz 

through  B-(162).   The  coefficient  A_  is  arbitrary,  as  was  the 

coefficient  F„  in  the  TE  case.   The  parameters  a,  g,  and  s  are  given 

TM      TM 
by  Eqs.  (26)  -  (28).   The  quantities  6        and  <f>   are  given  by 
j      -i  .  a.g  gs 

a5  II- (44) 


C-  W[(^K'J(^)] 


II-(45) 


It  may  be  noted  from  Eqs.  (29)-(30)  and  (44)-(45)  that 


w  C  -  KV"/)  w  ft1 


TF 
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II- (46) 


II- (47) 


In  the  course  of  obtaining  the  solutions  for  E   and  H   in 
Appendix  B,  the  characteristic  waveguide  and  cutoff  equations  for  the 


11 


TM  case  were  obtained: 

°  °  II-(48) 
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2   2 
Eq .  (49)  is  simi].ar  to  Eq .  (32)  except  for  the  factor  n  /n   in  the 

argument  of  the  arctangent.   From  Eqs.  (48)  and  (49),  it  may  be  seen 

that  the  value  of  W  for  the  TM  case  must  be  larger  than  W  in  the  TE 

TM    TE 
case  from  Eqs.  (31)  and  (32),  owing  to  the  fact  that  <j>   >  <b   and 

ag    ag 

TM    TE 
d>   >  <j>  J.   Finally,  Eq .  (35)  applies  equally  well  to  TM  modes  for 
gs    gs 

the  allowable  values  of  0  for  waveguiding  without  cutoff,  as 

determined  by  Eq.  (48). 

TE  +  TM  Modes 

As  mentioned  above  and  shown  in  more  detail  in  Appendix  A,  a 

field  of  any  polarization  in  a  homogeneous  asymmetric  planar  guide  can 

be  expressed  as  the  sum  (i.e.,  the  superposition)  of  a  TE  field  and  a 

TM  field.   Further,  the  components  obtained  from  the  TE  derivation 

[E  ,  H  ,  and  H  ]  obey  the  waveguide  and  cutoff  equations  given  by 

Eqs.  (31)  and  (32),  whereas  the  components  obtained  from  the  TM 

derivation  [E  ,  E  ,  and  H  ]  obey  Eqs.  (48)  and  (49).   For  propagation 
x   z       y 

of  a  TE  +  TM  wave  in  a  homogeneous  guide  of  thickness  W,  the  angle 

0  =  0   in  Eq .  (31)  is  not  equal  to  the  angle  0  =  0   in  Eq .  (48)  due 

TE   TM    TE       tm 
to  the  different  phases  d>   ,  0    ,  d>   ,  and  <f>     As  a  consequence, 

ag'   ag    gs'      gs 

the  propagation  constant  3  must  be  different  for  the  TE  and  TM  cases, 
i.e., 


$r  =  K»^%e  II_(50) 
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frx"^0"^  ^M  «>C51) 


Hence,  as  the  TE  +  TM  wave  propagates  down  the  guide  according  to 

-if3z 
e     ,  the  phase  difference  a3  between  the  TE  and  the  TM  parts  of  the 

TE  +  TM  wave  is 


Ap  ((3T£-(Ol  n-(52) 

=>6  n^  (i^v^re.-  W  9-m)  L        ii- (53) 

where  L  is  the  length  of  the  guide.   Thus  a  TE  +  TM  wave,  initially 

linearly  polarized,  will  become  more  or  less  elliptically  polarized, 

depending  on  0   and  0    [which  are  determined  by  Eqs.  (31)  and  (48)] 

and  L,  for  a  given  k  and  n  . 
o      g 

Concerning  the  minimum  thickness  (W  .  )__  ,  _,  necessary  to 

mm  IE  +  TM 

propagate  a  given  mode,  the  larger  of  (W  .  )__  or  (W  .  )_.,  [given  by 

mm  TE      mm  TM  °  J 

Eqs.  (32)  and  (49)]  must  be  the  minimum  thickness.   Owing  to  the 

2   2 
factor  n  /n   in  the  argument  of  the  arctangent  of  Eq.  (49),  (W  .  )_, 

g  a  °  min  TM 

is  the  larger  (since  n  >  n  ).   Hence  the  cutoff  thickness  for  TE  +  TM 
g    a 

modes  is 


00™=    Cw^,^ 


II-(54) 


Quasi-Geometrical  Ray-Optics  Approach 

It  is  possible  to  derive  the  waveguide  equations  using  a  quasi- 

3 
geometrical  approach  similar  to  one  given  by  Tien  and  Ulrich.    Referring 

to  Figure  II-2,  the  plane  wave  propagating  down  the  guide  may  be  thought 

of  as  a  zig-zagging  within  the  guide,  plus  an  evanescent  wave 
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®k— vJ-t°~e-~ ^ 


♦•2 


Figure  II-2.   Geometry  of  a  homogeneous  optical  waveguide, 
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propagating  in  the  superstrate  and  in  the  substrate  at  the  same 
velocity  in  the  z-direction  as  in  the  guide.   This  is  consistent  with 
the  statement  made  earlier  in  this  chapter  that  the  wavevector  component 

A 

along  the  z-axis,  3,  has  a  constant  value  k  n   sinQ  in  all  three 

0  g 

regions  (guide,  substrate  and  superstrate)  for  a  TE  or  for  a  TM  wave. 
Physically,  if  this  were  not  so,  eventually  a  point  would  be  reached 
along  the  guide  for  which  a  wave  within  the  guide  would  be 
unaccompanied  by  an  evanescent  wave  outside  the  guide  —  clearly  an 
impossibility  for  a  totally  internally  reflected  guided  wave. 

The  optical  rays  in  Figure  II-2  represent  normals  to  the  pro- 
pagating, zig-zagging  plane  wavefronts  within  the  guide.   It  may  be 
seen  that  the  time  it  takes  for  the  light  to  travel  along  the  path 
PQR  must  be  equal  to  the  time  it  takes  for  the  light  (the  evanescent 
wave)  in  the  n   -  region  (or,  equivalently ,  in  the  n   -  region)  to 
travel  along  the  path  PR.   Another  way  of  expressing  this  is  to  state 
that  the  phase  delay  experienced  by  the  light  wave  in  the  guide  must 
be  equal  (to  within  an  even  multiple  of  tt)  to  the  phase  delay 
experienced  by  the  evanescent  wave,  i.e., 


(akdsc    A-tl^)^    "t^mir  -    (^J^e  <id°y) 


n-(55) 


the  evanescent  wave's  phase  delay  is: 


(jpkue  fel^-  (P*)^  UWW9-)|3 


II-C56) 


The  phase  delay  in  the  guide  has  three  contributions,  one  due 
to  the  optical  path  along  PQR,  another  due  to  the  phase  change 
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~2^ac,   resulting  from  total  internal  reflection  from  the  n  -  n 
ag  g    a 

boundary,  and  the  third  change  -2<j)   due  to  total  internal  reflection 

gs 

from  the  n   -  n  boundary : 

g    s        } 

uw  Jd^3  =  CFm)>.«3  -s^a  -J&  n_(57) 

Substituting  Eqs.  (56)  and  (58)  into  Eq .  (55)  produces: 

/i3^9-=aii/3wM^^-3f3fc-a«^         n_(59) 


Since  3  =  kQn   sinO,  Eq .  (59)  can  be  written  as 


=  (svoJ?0is)  [  (i-4^3e)/e«a]  n-(6i) 

=  (3M^j)  c^e-  IM62) 


awi6n,  ci^>&-  a^p   -A«p   =a 


3s     o^-rr  n-(63) 


Eq .  (63)  is  identical  with  the  characteristic  waveguide  equations 

(31)  and  (48),  providing  that  the  correct  values  of  6        and  6       are  used 

Tag     Tgs 
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4 
According  to  Born  and  Wolf,   the  phase  changes  <f>   and  d>   for  the  TE 

ag     rgs 

and  TM  cases  of  total  internal  reflection  are  identical  to  the  values 

given  by  Eqs.  (29),  (30),  (44)  and  (45).   Further,  once  the  waveguide 

equations  are  known,  it  is  then  possible  to  obtain  the  waveguide 

cutoff  equations  for  (W  .  )    and  (W  .  )_,.   Finally,  the  extension  to 
mm  TE       mm  TM 

TE  +  TM  waves  is  the  same  as  given  above  —  by  the  superposition  of 

TE  and  TM  waves. 

Physically,  the  waveguide  equation  (63)  [i.e.,  Eqs.  (31)  or  (48)] 

represents  the  condition  that,  after  two  reflections  (one  at  each 

boundary),  the  wavefront  is  in  phase  with  the  original  wavefront.   If 

this  were  not  so,  many  wavefronts,  each  differing  in  phase  by  some  amount, 

might  eventually  add  up  somewhere  in  the  guide  such  that  they  would 

destructively  interfere.   By  conservation  of  energy,  the  light  would 

then  have  to  pass  out  into  the  n  or  n  regions,  i.e.,  the  light 

as 

could  no  longer  propagate  in  the  n   -  medium. 

It  may  be  seen  that  a  quasi-geometrical  ray-optics  approach  can 
be  a  useful  way  of  treating  wave  propagation  in  an  optical  waveguide. 
It  provides  a  means  for  understanding  the  physics  of  the  problem  at 
hand;  at  the  same  time  it  results  in  an  accurate  representation  of 
the  homogeneous-waveguide  equation. 


CHAPTER  III 
BASIC  PHYSICS  OF  INHOMOGENEOUS  ASYMMETRIC  PLANAR  OPTICAL  WAVEGUIDES 

Introduction 
An  inhomogeneous  optical  waveguide  is  one  that  has  a  variable 

index  of  refraction  n  within  the  guide.   The  refractive  index  profile 

g 

treated  in  this  dissertation  is  assumed  to  lie  in  the  x-direction 

only,  i.e.,  transverse  to  the  net  direction  of  propagation.   Further, 

the  index  profile  is  assumed  to  be  a  monotonic  function  of  x,  such 

that  the  value  of  the  index  is  higher  near  the  superstrate.   That 

is,  as  indicated  in  Figure  III-l,  the  refractive  index  of  the  guide  has 

its  maximum  at  the  n   -  n  boundary  and  then  continuously  decreases 
a    g 

in  the  negative  x-direction.   The  index  can  be  written  as 


fl^n30)  tti-cl) 


where 


n3W 


III-(2) 


Quasi-Geometrical  Ray-Optics  Approach 
It  is  convenient  to  begin  with  a  quasi-geometrical  approach. 
The  justification  of  the  results  will  he  discussed  afterwards. 
Letting  3  still  be  the  wavevector  component  along  the  z-axis,  then 
can  be  written  as 

17 
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x*0 


x-W  -I 


na0c)>  n,>tta     -kr  oJ/ 


3 


x*0 


>     n-  Cx) 


x=-vJ 


Figure  III-l.   Cross-sectional  view  of  light  propagating  in  an 

inhomogeneous  asymmetric  planar  optical  waveguide. 
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(3-  ho  a  Od  **B(x) 


III-(3) 

.   v  ;  ^am  rr:  x  i 


In  Appendix  C,  solutions  were  obtained  to  the  wave  equation  assuming 
that  6  was  a  constant  independent  of  x  for  both  TE  and  TM  waves,  i.e., 
as  in  the  homogeneous  case 


f 


III-C4) 


III-C5) 


where 


yos    »1flCx)&<^  B{k)~  Q>ufW' 


III-(6) 


Eq.  (4)  implies  that  the  wave  travels  at  the  same  velocity  v   in  the 

z  direction  in  all  three  regions:   guide,  superstrate,  and  substrate. 

Fq.  (6)  is  simply  another  way  of  writing  Snell's  Law. 

Referring  to  Figure  III-l,  it  may  be  seen  that  as  the  ray  travels 

from  point  P  to  point  Q,  the  ray  is  constantly  bent  away  from  the 

normal  to  the  n   -  n  boundary,  just  as  predicted  by  Eq .  (6).   Upon 
a    g 

reflection  at  the  n  -  n  boundary,  the  ray  travels  from  point  Q  to 
g    s 

point  R,  this  time  being  constantly  refracted  toward  the  normal.   At 

the  n   -  n  boundary  it  is  totally  internally  reflected  once  again, 
a    g 

In  order  to  obtain  the  characteristic  waveguide  equation  for  the 
inhomogeneous  case,  the  same  argument  given  in  the  homogeneous  case 
is  useful.   That  is,  the  time  it  takes  for  the  light  to  travel  along 
the  path  PQR  must  be  equal  to  the  time  it  takes  for  the  light 


20 


(an  evanescent  wave)  in  the  n  -  region  to  travel  along  PR.   As 
before,  this  can  also  be  expressed  in  terms  of  phase  delays: 

The  evanescent  wave's  phase  delay  is 


III-C7) 


III-(8) 


The  phase  delay  in  the  guide  has  three  contributions:   one  due  to  the 
changing  optical  path  caused  by  the  index  gradient,  /k  •  dr,  where 

A      A      A. A    A         A  . 

r  =  xx  +  yy  +  zz  (x,  y,  and  z  are  unit  normals):  the  other  two  are 

due  to  the  reflections  at  the  n  -  n  and  n  -  n  boundaries:  -2d> 

a    g      g    s  Yag 

and  -26      .   The  values  of  6        and  6        can  be  obtained  using  Eqs. 
gs  Tag     Ygs 

II-(29)  and  II-(30)  for  the  TE  case  and  Eqs.  II-(44)  and  II-(45)  for 
the  TM  case.   Defining 

%(x)  =    n3(>tut&(x)  m_(9) 


then 


III-(IO) 


III-(ll) 


III-(12) 


III-C13) 
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where 


o.-K  (y<)U 


2   .  =\'/a 


^=K  (n2'<) 


III-C14) 


III-U5) 


The  phase  delay  in  the  guide  can  thus  be  written  as 


(^*'t))3=  JJ*-d*  -af^-  s<p3. 


Further, 


III-(16) 


(PR) 


PaR  °  (FR)  ~" 


III-(17) 


The  first  integral  on  the  right-hand  side  of  Eq .  (17)  is  due  to  the  x- 
component  of  the  optical  ray  path  along  PQ.   [The  product  k  •  r  is 
negative,  hence  the  factor  (-dx) . ]   The  second  integral  is  due  to  the 
x-component  of  the  optical  ray  path  along  QR.   [The  product  k  '  r  is 
positive.]   The  third  integral  is  due  to  the  z-component  of  the  path  along 
PR.   [The  product  k  •  r  is  positive  also.]   Since,  by  Eq .  (6), 
n  (x)sinO(x)  is  a  constant,  Eq  .  (17)  becomes 

0 

jXa-r  =   2J?o  j  ,na(K)<uo06OdLx    +    (AQrj)  C^)  m.(18) 

Thus,  with  the  aid  of  Eqs.  (8),  (16),  and  (18),  Eq .  (7)  becomes 
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0 

j 


^(Fi?)  =  [^„  j  "3«  <*>«**)<*•<  +  (A^)(pr) 


III-(19) 


0 
PJfc      f  „   GO  tO0&6c)d*  -  «P<P      -  5><Pas-  5-^n- 

°  J   3  '  J       °  III-(20) 

— W 

Eq .  (20)  is  thus  the  waveguide  equation  for  asymmetric  inhomogeneous 
waveguides,  as  determined  by  the  quasi-geometrical  ray-optics  approach. 
It  may  be  noted  that  when  n  (x)  and  cosO(x)  are  constants   [as  in  the 
homogeneous  case],  Eq.  (20)  reduces  to  Eqs.  II-(31)  or  II-(48),  the 
waveguide  equations  for  homogeneous  waveguides. 
Writing 


cjk&(x)=    [/-  *^a960]6 


III-(21) 


then  the  integral  in  Eq.  (20)  can  be  transformed  to 

°y^  o  °y^  j       J  m-(22) 

UJr<o    J  L   3      /  J    A  III-(23) 


where  the  last  step  follows  from  Eq .  (6).   Thus  Eq .  (20)  becomes 
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which  makes  it  easy  to  directly  incorporate  the  mathematical  form  of 

the  index  gradient  n  (x) .   Sometimes  the  form  of  n  (x)  allows  the 
g  g 

integral  to  be  directly  integrated  analytically.   Otherwise  (such  as 
with  an  error  function  profile)  it  can  be  integrated  numerically. 
From  Eqs.  (10)- (15),  it  may  be  seen  that  Eq .  (24)  is  transcendental  in 
both  n  and  W.   Writing 


0  III-(25) 


then,  for  a  given  value  of  n  (0),  perhaps  the  easiest  way  to  solve 
Eq .  (24)  is  to  pick  a  value  for  0(0),  then  try  values  of  W  until  one 
is  found  which  satisfies  Eq.  (24).   By  using  different  values  of  0(0), 

a  set  of  curves  can  be  generated  for  the  different  modes  m.   Also, 

TM    TE      TM    TE 
since  is        >   A        and  d>    >  d>   ,  then  the  value  of  W  obtained  from  Eq . 
ag    ag      gs    gs' 

(24)  for  the  TM  case  is  greater  than  W  obtained  for  the  TE  case,  as 
in  the  homogeneous  case. 

For  the  cutoff  condition,  waves  will  begin  to  escape  when  either 


5^i9(c>)^  n^/m^a 


Swi9(-vo)*  ns/xi. 


III-(26) 


III-C27) 


where 


na o.-  nj(°)  ni-(28) 

A^~  "/-w)  ln.(29) 


24 


Using  F,q .  (6),  Eqs.  (26)  and  (27)  are  equivalent  to  the  following 
equations 


n=  f]  (x)-&nQ6c)~  na 


III-C30) 


n  -  f\  (x)  £w\  6)6*)  -  n  < 


III- (31) 


respectively.   It  was  assumed  for  homogenous  waveguides  that  n   >  n  . 

Hence  as  0  decreases  from  a  maximum  value  of  tt/2  (i.e.,  propagation 

straight  down  the  guide),  then,  for  a  given  functional  form  of  n  ,  the 

g 

cutoff  condition  first  to  be  met  is  given  by  Eq .  (31).   That  is,  waves 

will  start  to  leak  out  at  the  n  -  n  boundary  when  sinQ(x)  =  n  /n  (x) 

g    s        J  s'  gv 

[i.e.,  when  sinQ(-W)  =  n  /n  (-W) ,  or  equivalently ,  when  sinG(O)  =  ^  /n  (0) 
s   g  t»   g 

The  waves  will  still  be  totally  internally  reflected  at  the  n  -  n 
boundary.   This  result  is  identical  with  the  result  obtained  for  the 
homogeneous  case.   Substituting  Eq .  (31)  into  Eqs.  (10)  -  (15)  produces 


C=<=0 


7F 

5 


jA  '<? 


ay* 


=  h*Cl  [(«*-*+!  /(^  ~^J 


III-(32) 


III-(33) 


£  ^-'l^n^i'kin^fi  ni_m) 


Then  at  cutoff,  Eq .  (24)  becomes  for  the  TE  case 

o 


J  [<w-<]'/sa<  •   4irimlf  h 


buv\ 


-ly^) 


0        3\'6 


(«/' 


|>*S 


IIII- (35) 


TF 


and,    for   the   TM  case, 

o  , 


rmfT  + 


W\ 


3 
l3 


3\'/a 


".?    (nj-n?f  J^II-(36) 


As  mentioned  above,  Eqs .  (35)  and  (36)  hold  when  n   >  n  .   However, 

for  the  inhomogeneous  case,  the  situation  n   <  n   is  different  than 
b  s    a 

n   >  n   [since  n  (0)  4   n  (-W) ] .   Hence  if  n   <  n  ,  then  the  cutoff 
s    a         g       g  s    a 

TE       TM 
condition  is  given  by  Eq .  (30),  in  which  case  $  '  4   0,  (j)   f  0, 


i>    J  =  4>   =0,  and  Eq .  (24)  becomes 
ag    gs 


kiaa! 


OTiS-ncTjj     in-(37) 


and 


n,f    K1-/)" 


aV/a 


III-C38) 


It  may  be  seen  that,  indeed,  the  substitution  of  n   for  n  and  vice 
versa  in  Eqs.  (37)  and  (38)  does  not  result  in  the  same  relations  as 
given  by  Eqs.  (35)  and  (36)  because  n   4   n      ,    unlike  the  homogeneous 
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Since  sin0(0)  cannot  exceed  a  value  of  unity,  then  from  Eq .  (6) 

[or  Eq.  (25)],  n  can  never  exceed  the  value  n  (0).   In  fact,  n  cannot 

equal  n  (0)  either,  as  may  be  seen  from  the  following  argument.   If 
g 

sin0(0)  =  1,  then  0(0)  =  tt/2  ,  and  from  Snell's  Law  [Eq.  (6)] 
0(x)  =  tt/2   for  all  x.   Hence  the  integrand  in  Eq .  (20)  would  be  zero 
for  all  x,  and  in  order  to  pronagate  a  light  wave,  only  an  infinitely 
thick  (V  ->  °°)  guide  could  satisfy  Eq .  (20),  which  is  physically 
unrealistic.   Thus,  using  Eq .  (25),  the  range  of  values  of  0(0)  for 
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waveguiding  in  the  case  n   >  n   is  given  by  the  relationship 


n  >f  ^s 


III-(39) 


and  for  the  n   >  n   case, 
a    s 


ki  <>  ^>n, 


III-(40) 


Rqs.  (39)  and  (40)  applv  to  both  TE  and  TM  modes  as  determined  by  Eq . 
(24).   It  may  be  seen  that  Eq .  (39)  is  identical  to  Eq .  II-(35)  for 
the  homogeneous  case. 

Solutions  for  E  and  H 
As  in  Appendix  B,  in  Appendix  C  trial  solutions  are  chosen  for 
the  vector  potentials  F  and  A,  only  now  a  provision  is  made  in  the  form 
of  the  solution  to  allow  for  an  index  gradient  within  the  guide. 
Specifically,  instead  of  the  product  gx  in  Eqs.  II-(24)  and  II-(42), 
where  g  is  a  constant  given  by  Eq .  II-(27),  the  choice  is  g  =  g(x).   The 
function  g(x)/k   is  often  called  the  eikonal.    Using  the  boundary 
conditions   on  E  and  H  together  with  the  assumption  that  the  index 
in  the  guide  is  a  slowly  varying  function  of  x,  the  following  solutions 


ire  obtained  in  Appendix  C  for  TE  modes: 


III-(41) 


III-(42) 


{  a  f  '%»)  Fae";^  ^  (fj  +mrr)e: 


>fy+*J 


Irsr   X^-W  itt_(43 


(43) 
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where  Eqs.  II- (17)  through  II- (21)  apply  equally  well  to  the  TE 
inhomogeneous  case.   The  exact  form  of  Hqx  and  Hqz  is  given  in  Eqs. 
C-(135)  through  C-(138).   It  may  be  seen  from  Eq .  (42)  that  both  the 
amplitude  and  phase  of  E   (x)  are  functions  of  x.   The  function  g(x) 


is  given  by 


3Gc)=\  j  ?W* 


IIT-(44) 


where  q (x)  is  given  by  Eq.  III-(9).   The  method  used  in  Appendix  C 
for  solving  the  wave  equation  for  an  inhomogeneous  medium  is  called  the 
WKBJ  method  "    (after  Wentzel,  Kramers,  Brillouin,  and  Jeffries).   The 
assumption  that  the  guide  has  a  slowly  varying  index  profile  leads  to  the 
follwing  requirement: 


i 


^•^  '  III-(45) 


where  q'  and  q';  are  the  first  and  second  derivatives  of  q  with  respect 
to  x.   Eq.  (45)  is  the  WKBJ  approximation  for  the  TE  case.     It  may 
be  seen  that  for  large  wavenumbers  k   (i.e.,  small  X  ) ,  the  WKBJ 


app 


roximation  is  more  easily  satisfied.   However,  as  q  =  n  cosO 


approaches  zero  (as  n  decreases  and  0  approaches  tt/2),  Eq .  (45)  cannot 
be  satisfied,  regardless  of  how  small  q'  and  q"  are.  This  situation 
will  be  discussed  in  Chapter  IV.  When  Eq .  (45)  is  satisfied,  all  of 
the  TE  results  obtained  from  the  quasi-geometric  ray-optics  approach 
are  valid,  at  least  to  the  level  of  approximation  given  by  Eq .  (45). 
In  the  homogeneous  case  q  is  a  constant.  Hence,  the  function  g(x)  in 
Eq.  (44)  can  be  easily  evaluated  to  give  (kQq)x  =  (k  n  cos0)x,  which 


28 


is  the  same  as  the  product  gx  originally  obtained  in  Chapter  II,  as 
mentioned  above. 

For  the  TM  case,  the  following  solutions  are  obtained  in  Appendix 

/        D™ 

III-(46) 
SW  =f-a^h  W  V  V   «[gW+^  .^  ni_(47) 

^Wi'HV*1"3  wif«?^)wx^w  ITI.(48) 


where  Fqs.  IT-(36)  through  II-(40)  apply  to  the  TM  inhomogeneous  case. 

The  exact  form  of  E   and  E   is  given  in  Eqs.  C-(224)  through  C-(229). 
ox      oz 

It  mav  be  noted  that  the  amplitude  dependence  of  H   on  x  has  two 

r  oy 

-1/2 
contributions,  n  (x)  and  q     (x) ,  whereas  for  the  TE  amplitude 
g 

-1/2 
dependence  only  one  factor,  q     (x) ,  is  present.   The  WKBJ  approximation 

for  the  TM  case   is  given  by 


Jl 


s/iy-^l.-vW-^ 


.3  /   o-  3 


zfi     'V  i       r  )  n3   \r^ 


«        IIT-(49) 


which  has  the  additional  terms  (given  in  the  second  brace)  due  to  the 

first  and  second  derivatives  of  n  .   The  situation  when  Eq .  (49)  is 

g 

not  satisfied  will  be  discussed  in  Chapter  IV.   When  Eq .  (49)  is_ 
satisfied,  all  of  the  TM  results  obtained  from  the  quasi-geometric 
method  are  valid,  at  least  to  the  level  of  approximation  given  by  Eq . 
(49). 


CHAPTER  IV 
BASIC  PHYSICS  OF  TOTAL  INTERNAL  REFRACTION 
IN  INHOMOGENEOUS  PLANAR  OPTICAL  WAVEGUIDES 

Introduction 
It  was  stated  in  Chapter  III  (and  shown  in  Appendix  C)  that  the 
WKBJ  solutions  for  E  and  H  are  valid  as  long  as  Fqs.  Ill- (45)  and 
III- (49)  for  TE  and  TM  modes  respectively   are  satisfied.   When  either 

Eq.  III-(45)  or  III-(49)  is  not  satisfied,  it  may  be  seen  from  Eqs. 

-1/2 
III-(42),  C-(135),  III~(47),  and  C-(225)  that,  due  to  the  factor  q 

at  least  two  field  components  (E   and  H   for  the  TE  case,  H^  and  E^ 

for  the  TM  case)  become  infinite  as  q (x)  =  n  (x)  cosO(x)  approaches  zero 

as  0(x)  -*■  it/2.   Since  this  is  physically  unrealistic,  some  other 

approach  must  be  used  to  solve  the  wave  equation  near  the  level  q  =  0. 

11         12 
In  the  field  of  radio  ionospheric  physics,  Budden   and  Wait   have 

developed  a  useful  approach  for  describing  the  propagation  of  light 
near  q  =  0.   This  method  will  be  adapted  to  optical  waveguides.   Before 
doing  this,  however,  it  is  convenient  to  use  quasi-geometrical  arguments 
to  form  a  physically  based  explanation  of  the  results  obtained  in  the 
more  rigorous  derivation.   This  permits  a  practical  understanding  of 
the  phenomenon  of  total  internal  ref raction-and-ref lection  in  general- 
index-gradient  planar  optical  waveguides. 

Quasi-Geometrical  Ray-Optics  Approach 
First  of  all,  it  should  be  mentioned  that  -W  is  now  chosen  to  be 
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the  value  of  x  for  which  0  =  tt/2  (and  q  =  0)  for  a  given  mode  m  and 
incident  angle  0(0),  as  indicated  in  Figure  IV-1 .   That  is,  at  the 
level  x  =  -W,  0(-W)  =  tt/2,  q(-W)  =  0,  and  total  internal  refraction 
occurs.   It  is  convenient  to  let 

'3  iv-  (i) 


where 


n^Co)-  r\a+t«  IV_(2) 


at  the  top  of  the  guide,  and 


Y\^L~b)=r}0  iv-o) 


where  -D  is  the  penetration  depth  of  the  index  profile  such  that  beyond 
x  =  -D  there  is  no  longer  an  index  gradient.   For  instance,  nQ  might 
be  the  index  of  the  film  before  ion  implanation,  and  An  the  maximum 
index  change  (where  An  <<  n  )  which  is  generally  a  measurable  quantity. 
The  function  b(x)  thus  gives  the  index  profile,  where 


b(o)=l  IV_(4) 

b(~»)-0  IV- (5) 
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K=0 


*=-H/ 


*■  2 


Figure  IV-1.   Cross-sectional  view  of  light  propagating  in  an  inhomo- 

geneous  planar  optical  waveguide  by  means  of  total  internal 
refraction  and  reflection. 
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As  shown  in  Figure  IV-1,  after  reflecting  from  the  n   -  n  (0) 
boundary,  the  wave  is  continually  refracted  away  from  the  normal  to 
the  boundary  until,  at  the  level  x  =  -W  (and  6  =  tt/2)  ,  the  wave  is 
refracted  back  toward  the  normal  as  it  returns  toward  the  n&  -  n„(0) 
boundary.   The  phase  change  suffered  by  this  wave  has  three  contributions: 

a)  the  integral  of  the  changing  optical  path  k  ■  dr  over  the  path  POR; 

b)  the  phase  change  -2d)   upon  total  internal  reflection  at  the  n  -  n  (0) 

)  e       3.%  a      g 

boundary:  and  c)  the  phase  change  upon  total  internal  refraction.   It  is 
the  third  contribution  which  is  the  more  elusive  quantity. 

From  one  point  of  view,  since  the  total  internal  refraction 
phenomenon  is  purely  refraction,  no  phase  change  might  be  expected. 
On  the  other  hand,  since  the  light  does  in  fact  reverse  its  direction, 
then  a  phase  change  of  it  might  seem  reasonable,  since  in  the  limit  as 
0  ->  tt/2,  the  phase  change  2(f)   due  to  total  internal  reflection  approaches 
tt.   However,  neither  of  these  points  of  view  is  correct.   In  fact, 

the  phase  change  is  tt/2,  not  0  or  it. 

1  3 
It  is  known   that  a  phase  change  of  tt  occurs  when  light  passes 

through  the  focus  of  a  spherical  lens:  for  a  cylindrical  lens,  there  is 

a  phase  change  of  tt/2.   In  Figure  IV-2a,  it  may  be  seen  that  as  a  beam 

of  light  is  refracted  at  a  boundary  with  an  abrupt  refractive  index  change, 

the  beam  width  decreases  when  n   <  n  .   For  a  continuously  varying 

index,  a  similar  phenomenon  occurs,  as  shown  in  Figure  IV-3.   Here  the 

path  of  the  ray  is  shown  surrounded  by  a  small  beam  width.   Because  the 
waveguide  is  a  planar  slab  (as  opposed  to  a  fiber),  it  may  be  seen  that 
in  the  region  x  =  -W  and  (q  =  0) ,  the  beam  undergoes  a  cylindrical  focus. 
Hence,  it  is  reasonable  to  expect  a  phase  change  of  tt/2  upon  "reflection", 
i.e.,  total  internal  refraction,  at  this  level. 
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Figure  IV-2a 


Figure  IV-2b 


Figure  IV-2.  Refraction  of  light  at  a  boundary  with  an  abrupt  change 
of  refractive  index  (Figure  IV-2a)  and  in  a  medium  with 
a  continuously  varying  index  (Figure  IV-2b) . 
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Figure  IV-3.   Total  internal  refraction  in  a  planar  optical  waveguide 
showing  the  cylindrical  focussing  effect. 
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From  these  various  phase  contributions  the  characteristic  wave- 
guide equation  is  thus  obtained: 

o 

which  applies  to  both  TE  and  TM  modes.   [A  second-order  correction  for 

TM  modes  will  be  discussed  later.]   Comparison  of  Eq.  (6)  with  Eq.  III- 

(20)    shows   that   the   only   change   is    to   replace   the   term   2<j>        with   tt/2. 

8s 

Eq.  (6)  can  also  be  written  as 

o 


s\  { [  t^fo-fr** -  a<JL,  -  f  --£„-*  IV-(7) 


similar  to  Eq.  Ill- (24).  Eqs.  (6)  and  (7)  are  transcendental  in  n, 
which  is  equivalent  to  saying  that  it  is  transcendental  in  0(0), 
according  to  Eq .  Ill- (25).  It  may  or  may  not  be  transcendental  in  W 
too,  depending  on  the  functional  form  of  n  (x)  and  the  integral  in 
Eq.  (7).  Eqs.  (6)  and  (7)  are  useful  for  describing  the  propagation 
of  light  in  such  optical  waveguides  and  has  general  applicability  to 
any  slowly  varying  monotonic  index  profile  (including  the  error  function 

profile);  it  also  allows  computation  of  numerical  answers.   As  before, 

TM    TE 
since  i>        >  <j>    ,  the  value  of  W  for  the  TM  case  will  exceed  W  for  the 
ag    ag 

TE  case.   It  is  interesting  to  note  that,  in  the  limit  as  0(0)  ■*  tt/2, 

14 
Eq .  (7)  approaches  an  equation  for  TE  modes  derived  by  Marcuse. 

Following  arguments  given  in  Chapter  III,  the  cutoff  condition 

occurs  when  either  [from  Eq .  Ill- (30)] 
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^=>v  IV-(8) 

or    [similar   to   Eq.    Ill- (31)] 

fj~r\0  IV-  (9) 

For  n   >n   (which  is  usually  the  case) ,  the  cutoff  condition  is  given 
o   a 

by  Eq.  (9),  i.e.,  when  n  =  n  .   When  n  <  n  ,  the  light  can  no  longer 
be  totally  internally  refracted,  and  the  light  escapes  below  the  guide. 

For  this  cutoff  condition,  there  is  thus  no  tt/2  phase  change.   For  the 

TF        TM 
exact  form  of  the  phase  changes  2*  J  and  2<j>   ,  Eq .  (9)  can  be  substituted 

ag       a¥, 

into  Eqs.  III-(IO),  III-(12),  and  III-(U).   Further,  from  Eqs.  III-(6) 
and  III-(9),  q (x)  can  be  written  as 

'/a 

IV- (10) 


qlx)-    h^*)'*}9] 


Using  Eqs.  (2),  (9),  and  (10),  then 

'/a 


^(o)=  (Sn0  A*) 


IV-(ll) 


where  the  (An)   term  has  been  neglected  compared  to  2nQAn.   Thus,  at 


cutoff  the  phases  are 

.Tt~ 


ft 


-  it 


0L*\ 


<-* 


~a*\ 


*0+  gn>7   fop  ~"*J 


IV- (12) 


IV-(13) 


Eq.  IV- (7)  then  becomes  an  integral  equation  for  Wmin: 


IV-(14) 
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where  for  (W  .  )__.  Eq.  (12)  is  to  be  used  in  Eq.  (14),  and  for 
mm  TE 

(W  .  )™  Eq.  (13)  is  to  be  used. 

min  TM 

In  the  event  that  n  <  n  ,  a  much  less  common  situation  [where 
o    a 

n  <  n  +  An,  in  order  to  insure  total  internal  reflection  at  the 
a    o 

n  -  n  (0)  boundary],  then  the  cutoff  condition  is  given  when  n  =  n 
a    g 
[Eq.  (8)],  and  the  light  leaks  out  into  the  superstrate.   In  this  case 

i,        =0.   Thus  Eq.  (7)  becomes 

j  LnaM-n^]  ^=   Aft(Wl)/?  iv-(i5) 

which  is  the  same  for  both  TE  and  TM  modes  (except  for  the  second- 
order  correction  for  the  TM  case,  to  be  discussed). 

As  explained  in  Chapter  III,  n  can  never  exceed  or  equal  the 

value  n  (0).   Thus  the  range  of  values  of  0(0)  for  waveguiding  without 

g 
cutoff  in  the  case  n   >  n  is  given  by  the  relationship 

yi0-|-An  >  Kl  >  n0  IV-(16) 

and  for  n   <  n 

no  t&„  ;>  y  >  na  iv-u7) 

Since  0(-W)  =  it/2,  then  from  Eqs.  (9)  and  III-(6),  n  =  n  (-W)  =  nQ. 
Referring  to  Eq.  (1),  this  can  be  satisfied  only  when  b(-W)  =  0,  which 
from  Eq.  (5)  shows  that,  at  cutoff,  W  equals  the  maximum  index  gradient 
penetration  depth  D. 

It  was  mentioned  above  that  there  is  a  second-order  correction  to 
the  characteristic  waveguide  equation  for  TM  modes.   Physically,  the 
reason  for  this  is  as  follows.   In  the  TE  case,  the  only  component 

of  "E   is  E  ,  which  is  transverse  to  the  index  gradient.   Hence,  as 
TE     y' 

a  TE  wave  propagates  within  the  guide,  the  E-vector  at  any  one  point, 
say  x  ,  within  the  guide  will  "see"  only  one  index  —  ng(x0)  "~  at  a 
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time.   For  the  TM  case,  however,  since  there  is  a  component  of  E 

along  the  index  gradient  (the  E  component),  the  situation  becomes 

more  complicated.   Thus  when  deriving  the  tt/2  phase  change  for  the 

characteristic  waveguide  equation  for  the  TM  case,  it  was  implicitly 

assumed  that  the  TM  wave  "sees"  some  sort  of  average  index  n  (x  )  at 

go 

at  any  given  position  x  on  the  ray  path.  As  shown  below,  this 

averaging  takes  place  over  some  fraction  of  a  wavelength  of  light, 

say  C.A(x  ),  where  C  is  a  constant  to  be  determined  and  where  A(x)  = 

X  /n  (x) .   Any  shift  5W  in  the  level  of  total  internal  refraction  at 
o  g 

x  =  -W  can  then  be  considered  as  due  to  the  change  of  wavelength  AA 

caused  by  a  different  refractive  index  above  and  below  the  ray  path. 

That  is,  the  index  above  the  ray  path,  say  n  (-W),  is  slightly 

greater  than  n  (-W) ,  and  the  index  below  the  path,  say  n  (-W) , 

is  slightly  less  than  n  (-W) .   Thus  the  wavelength  A   slightly  above 

the  ray  path  is  smaller  than  A(-W),  and  the  wavelength  A  slightly 

below  the  ray  path  is  slightly  larger  than  A(-W).   Hence  the  median 

position  of  A  (over  which  E  is  averaged)  occurs  at  a  level  slightly 

below  x  =  -W.   The  net  effect  is  that  the  E  component  sees  the  index 

n  (-W)  before  the  ray  actually  reaches  the  level  x  =  -W.   As  a  result, 

the  TM  ray  is  internally  refracted  "sooner"  than  in  the  TE  case,  i.e., 

at  a  value  of  |x|  which  is  slightly  smaller  than  in  the  TE  case. 

This  has  a  slight  tendency  to  offset  the  fact  that  the  value  of  W  for 

the  TM  case  (from  the  waveguide  equation)  is  larger  than  W  for  the  TE  case. 

It  is  possible  to  derive  the  functional  form  of  this  shift  6W. 
From  the  above  argument,  it  is  reasonable  to  expect  that  6W  is  pro- 
portional to  the  change  in  wavelength  AA  due  to  the  slightly  different 
+ 


re 


fractive  indexes  n  and  n  ,  i.e., 
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sw=eaw 


IV- (18) 


where  C  is  another  constant.   Since  A  =  A  /n  ,  then 


hX  = 


-A, 


3 


M 


IV-(19) 


Further,  for  an  index  gradient  dn  /dx  near  the  level  x  =  -W,  the 
change  in  index  seen  over  a  distance  C  A  is 


An    =    fl.^l 


IV-(20) 


ii. 


=  C'A»Vn3W    <WUW 


IV-(21) 


Combining  Eqs.  (18)  -  (21)  produces 


x=->V 


IV-(22) 


Using  Eq.  (1) 

Ax" 


=  An  fc'(-w) 


*=-W 


IV- (2 3) 


n^"w) =  ^o +  5rlo  ^"  +  S*oM  -h^J 


IV-(24) 


Since  An  <<  n  ,  then  Eq.  (22)  becomes  approximately 


^W 


$  _q  n    \Q    Ak?  bY-wj 


IV-(25) 


Writing   A' 


2        2 
Att   /k      ,    then 
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.   ,  IV- (26) 


'o 


The  product  C  C.  is  evaluated  in  the  more  rigorous  derivation  given 

2 
in  Appendix  D  and  has  the  value  3/8tt   [from  Eq.  D-(122)].   Thus 

°        5     tf  3^  5  IV-(27) 


In  order  to  get  a  rough  estimate  of  how  large  6W  might  be,  if 

An  =  0.1,  D  =  lum,  and  W  =  D  (i.e.,  the  wave  is  propagating  nearly  at 

cutoff),  then, for  a  linear  index  profile,  b'(-W)  =  O.lum   ,  which  is 

a  fairly  steep  profile.   Further,  letting  A  =  0.63ijm  and  n  =  1.5, 

then  6W  -  4.5  x  10  um  =  0.45  A,  which  is  very  small.   For  less  steep 

profiles,  6W  will  be  smaller  still.   Thus  to  a  very  good  approximation 

6W  can  be  neglected  in  calculations  concerning  the  TM  characteristic 

waveguide  equation  as  well  as  the  cutoff  thickness  (W  .  )__,. 

mm  TM 

Solutions  for  E  and  H 

1.   The  Linear  Approximation 

As  mentioned  in  the  introduction  of  this  chapter,  the  WKBJ 

solutions  for  the  field  components  E   and  H   (for  the  TE  case)  and 

r        oy      ox 

H   and  E    (for  the  TM  case)  become  infinite  as  q(x)  approaches 
oy      ox 

zero.   In  order  to  get  around  this  problem,  theoretical  techniques 

11         12 
developed  by  Budden   and  Wait   for  radio  ionospheric  physics  can  be 

adapted  to  optical  waveguides. 

The  assumption  is  now  made  that,  for  a  slowly  varying  index  pro- 

2 
file,  q  is  approximately  linearly  proportional  to  x  near  the  level  q  =  0, 

i.e.,  near  x  =  -W.   This  is  nearly  the  same  thing  as  saying  that  n  (x) 

is  linear  in  this  region,  as  may  be  seen  from  the  following  argument. 
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Expanding  n  (x)  in  a  Taylor  series  expansion  around  x  »  -W,  produces 

3  "  a 

IV-(28) 
which,  from  Eq.  (1),  becomes 

0 


which  is  linear  in  x.   Further,  for  small  An 

fl/6c)s  nf  +  2*0  An  [i^j+bY-^x+w)]      iv_(30) 

which  is  also  linear  in  x.   At  the  level  q  =  0  (and  x  =  -W) ,  0(-W) 
has  the  value  tt/2,so  from  Eq.  III-(6) 


a=,  ki  7-W)  s^s96w) 


1  '"a 


IV-(31) 


-nV-w) 

"^  L    '  IV-(32) 


Thus,  from  Eq.  (30), 


3         ~"       -       4*b(-ltf) 


IV-(33) 


hence  [from  Eq.  (10)] 


<,*(*)-     »*(*)-?  =  3o.&»  b'H(x+w)  IV.(34) 


which  is  linear  in  x  too.   Eq.  (34)  is  equivalent  to  a  Taylor  expansion 


2 
of  q  (x)  in  the  region  near  x  =  -W. 
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2.   TE  Modes 

Defining  the  following  quantity: 


IV- (35) 


2  2 

where  (q  ) '  is  the  derivative  of  q   [from  Eq.  (34)]  with  respect  to  x, 

then,  as  shown  in  Appendix  D,  the  wave  equation  for  F   [Eq.  A-(92),  or 
equivalently,  Eq.  D-(20)]  can  be  written  as 


IV-(36) 


which  is  called  the  Stokes  differential  equation.   Its  solution  as 
discussed  by  Budden   is  the  following: 


F  =  K     e  AC&; 


X  'Vlr  "^/  IV-(37) 


where  Km   is  a  constant  and  where  A±(0    is  the  Airy  integral  function 
TE  

defined  as  follows: 


>w<?)=  ir\  ^fe^  ^^)d^ 


IV-(38) 


Tables  of  Ai(c)  are  given  by  Miller.    A  plot  of  Ai(c)  versus  X,   is 
given  in  Figure  IV-4.   It  may  be  seen  that  beyond  the  region  £  =  0 
(q  =  0,  x  =  -W)  the  solution  Ai(c)  is  an  exponentially  damped  wave  — 
an  evanescent  wave. 
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Osc.'ILt 


+  *(?) 


Figure  IV-4.   Plot  of  the  Airy  integral  function  Ai(c)  as  a  function 

2  11,12 
or  c,  or,  equivalently,  as  a  function  of  q  . 
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There  is  another  solution  to  Eq.  (26)  since  it  is  a  second-order 
differential  equation.   It  is  another  Airy  integral  function  Bi(^). 
However,  beyond  the  level  x  =  -W  at  q  =  0,  the  function  Bi(?)  becomes 
indefinitely  large   as  ?  increases,  i.e.,  as  x  ■>  -°°.   Since  it  is 
physically  unrealistic  for  F   (and  thus  E)  to  get  larger  and  larger  at 
distances  farther  and  farther  away  from  the  waveguide,  then  the  solution 

to  Eq.  (26)  cannot  contain  a  multiple  of  Bi(^). 

2 
One  of  the  interesting  properties  of  Ai(c),  or  equivalently  Ai(q  ), 

is  that  for  large  q  (corresponding  to  large  cosQ,  far  away  from  the 

level  q  =  0) ,  Ai(c)  asymptotically  approaches  the  WKBJ  solution  for  E, 

as  shown  by  Budden.    Referring  to  Figure  IV-5  and  following  the 

11  12  2 

arguments  of  Budden   and  Wait,    it  is  assumed  that  q   is  linear  over 

the  entire  region  be.   There  is  a  certain  range  cd_  where  |q|  is  so 

small  that  it  violates  the  WKBJ  approximation,  so  the  WKBJ  solutions 

cannot  be  used.   Outside  this  range   the  WKBJ  solutions  are  good 

approximations.   In  the  regions  be  and  cle  the  asymptotic  approximations 

of  Ai(<;)  can  be  used,  and  they  must  be  matched  to  the  WKBJ  solutions. 

In  other  words,  outside  the  region  near  q  =  0,  the  WKBJ  solutions 

should  be  valid;  inside  the  region,  Ai(?)  must  be  used.   At  the 

"boundaries"  of  this  region,  at  say  x  =  -W  ±  A  (where  A  is  to  be 

determined) ,  the  two  solutions  must  merge  smoothly,  thus  providing 

boundary  conditions  to  be  satisfied.   When  this  is  done,  as  shown  in 

Appendix  D,  the  phase  change  6   suffered  by  the  wave  upon  propagating 

2 
into  the  linear  n   (i.e.,  linear  q  )  region,  reflecting  (i.e., 

g 
internally  refracting)  at  the  level  x  =  -W,  and  then  propagating  out  of 

the  linear  n  region,  is  given  by 
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Figure  IV-5.   Plot  of  q  (x)  as  a  function  of  x  showing  the  linear 
region  be. 
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IV- (39) 


e  *  =  i   e     '""* 


Aside  from  the  factor  i  =  e17r   ,  this  phase  change  6R  is  nothing  more 
than  what  would  be  predicted  by  the  WKBJ  method  even  though  the  WKBJ 
solution  itself  breaks  down  in  the  region  q  =  0.  When  6  is  combined 
with  the  phase  change  suffered  by  the  wave  in  the  rest  of  the  guide, 
the  characteristic  waveguide  equation  [Eq.  (6)]  results. 

The  effect  of  curvature  of  the  index  profile  must  be  investigated 
in  order  to  determine  the  validity  of  Eq.  (34)  and  the  resultant 
conversion  of  the  wave  equation  to  the  Stokes  equation.   If  an 
additional  term  is  kept  in  Eq.  (29),  then 

2 


IV-(40) 


and 


IV-(41) 


Then,  from  Eq.  (10) , 


where 


lV-(42) 


IV-(43) 


IV- (44) 
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When  Eq.  (42)  is  substituted  into  the  wave  equation  for  ?x  [Eq.  A-(92)], 
then  as  shown  in  Appendix  D  the  wave  equation  can  be  converted  to  the 
Stokes  equation  for  which  Ai(?)  is  solution  only  if 

bi-i5?^B=r/3Bi| 

1  V  l  IV-(45) 


and  this  must  hold  for  values  of  t,   large  enough  for  the  asymptotic 

approximation  of  Ai(c)  (and  thus  the  WKBJ  solutions)  to  be  valid. 

According  to  Budden,    this  will  happen  when  \t,\    "»   1.  For  |c|  =  1, 
Eq.  (44)  requires  the  following  linearity  criterion: 


which  from  Eqs.  (1),  (43),  and  (44),  requires  that 

«:'(-*)«*  fa»*K9f  [»*&/* 


3 


IV- (47) 


The  condition  that  |?|  =1  (at  x  =  -W  ±  A,  the  boundaries  of  the  linear 

n  region,  as  shown  in  Figures  IV-4  and  IV-5),  corresponds  to 
g 

1 


la*?n0  n£(-»)l 


It  is  now  possible  to  write  down  the  solutions  for  ETE  and  HTE. 
The  solutions  for  x  =  0  are  given  by  Eqs.  III-(41),  C-(135) ,  and  C-(136) . 
For  the  region  (-W  +  A)  =  x  =  0,  the  solutions  are  the  WKBJ  solutions 
given  by  Eqs.  III-(42),  C-(135) ,  and  C-(137).   For  the  region 
(_W  -A)  =  x  =  (-W  +  A) ,  the  solutions  involve  the  Airy  integral  function: 
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y*)  =  f  ^  Ai® 


IV-(49) 


where  1C,  is  defined  in  Eq.  D-(58)  [using  Eq„  D-(48)].   The  exact 

form  of  H   and  H   is  given  in  Eqs.  D-(94)  and  D-(96) .   Finally,  for 
ox      oz 

the  region  x  =  (-W  -  A) ,  the  WKBJ  exponentially   decaying  solution  is 


% 


i    ~  IV- (50) 


where  F,  is  defined  in  Eq.  D-(62).   The  exact  form  of  Hqx  and  Ho 
in  this  region  is  given  in  Eqs.  D-(99)  and  D-(100). 

3.   TM  Modes 

7,11 
Defining  the  following  quantities: 


IV-(51) 


^)l  \     Ml 


±x    J    ai9n?    ^xs   IV-(52) 


'3     L    ^*     J      oUl.rjj 

(where  Q  is  an  "effective  value  of  q") ,    then,  as  shown  in  Appendix 
D,  the  wave  equation  for  A   [Eq.  A- (93)]  can  be  simplified  to 


IV-(53) 


similar  to  Eq.  D-20.   By  making  the  assumption  that  n  (x)  varies 

linearly  near  the  level  x  =  -W,  then  it  is  shown  in  Appendix  D  that 

the  difference  between  Q2  and  q2  from  Eq.  (52)  is  proportional  to 

[An  b'(-W)/k  ]2,  which  is  a  very  small  quantity  for  optical  waveguides. 

Hence,  near  the  level  x  =  -W,  the  wave  equation  Eq.  (53)  becomes 

approximately 
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IJJ   +  \S  %%)  V=0  IV-<54) 

and  all  the  techniques  described  for  the  TE  case  can  be  used  to  convert 
this  form  of  the  wave  equation  to  the  Stokes  equation,  and  to  then 

find  the  Airy  integral  solutions. 

2      2 
Neglecting  the  difference  between  Q  and  q  is  equivalent  to  the 

assumption  that  they  have  a  zero  at  nearly  the  same  point.   As  discussed 

in  the  quasi-geometrical  approach,  it  is  possible  to  calculate  this  (very) 

slight  difference  SW  between  the  TE  and  TM  reflection  points  (i.e., 

the  points  of  total  internal  refraction) .   Since  5W  is  so  small, 

Eq.  (54)  is  valid  for  TM  modes  to  a  very  good  approximation. 

2 
The  quantity  Q  (x)  is  subject  to  two  linearity  criterions,  as 

2  2 

compared  to  just  one  [Eq.  (46)]  for  q  (x) .   Assuming  Q  obeys  the 

criterion  given  by  Eq.  (46)  concerning  the  second  derivative  n  "(-W), 

2 
then  a  Taylor  series  expansion  of  Q  from  Eq.  (52)  can  still  result  in 

a  quadratic  term  which  has  the  first  derivative  n  ' (-W)  only  [Eq.  D-(123) ] . 

g 
2 
As  a  result   [from  Eq.  D-(128)],  Q  may  be  written  as 


GL*W=  (y±Yi)Gl.  ~(x+rtfQ 


a  IV- (55) 


where  6W  has  been  neglected  and  where  the  constants  Q»  and  Q„  are 
defined  in  Eqs.  D-(120)  and  D-(126)  respectively.   Eq.  (55)  is  of  the 
same  form  as  Eq.  (42);  thus  the  quadratic  term  in  Eq.  (55)  is  also  subject 
to  the  linearity  criterion  given  by  Eq.  (46),  i.e., 


=  la/3 


A  I  °  IV-(56) 
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Using  the  definitions  of  (^  and  Q2  together  with  the  fact  that  the 
minimum  value  of  n  is  n  before  cutoff,  Eq.  (56)  requires  that 

^ »  <~*    ^  h'('w^ 


Budden   defines  the  quantity  B  as  follows: 


BH 


or  k/K(x)J 


U— W  J  IV-(58) 


[where  optical  waveguide  notation  has  been  used  for  the  right-hand  side 
of  Eq.  (58)  instead  of  Budden's  radio  ionosphere  notation].   Using  Eq. 
(58),  Eq.  (57)  becomes 

8»/-/7  IV-(59) 

For  optical  waveguides,  Eq.  (58)  almost  always  satisfies  Eq.  (59). 
That  is,  if  the  linearity  criterion  given  by  Eq.  (46)  is  satisfied  for 
TE  modes,  then  it  is  satisfied  for  TM  modes  also,  to  an  extremely  good 
approximation.   For  instance,  if  the  index  profile  is  exponential  such 
that 

Y)(x)=    n    +  An[|-  e  \    (l-e    ) 

'§'  ^  J  '  IV-(60) 

then  Eq.  (58)  satisfies  Eqs.  (1)  -  (3).   If  the  steep-index-profile 
values  used  with  Eq.  (27)  are  chosen,  and  the  e~   point  of  the  index 
profile  occurs  when  x  =  -D/2  =  -  0.5ym,  then  d  =  2um  .   Also,  letting 
n  =  n  =1.5  and  W  =  D,  (i.e.,  the  wave  is  propagating  nearly  at  cutoff), 
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then  b'(-W)  =  b'(-D)  =  -  And/(1  -  edD)  =  0.313.   Thus,  from  Eq.  (56), 
B  =  51,  which  certainly  satisfies  Eq.  (57).   For  waveguide  index 
profiles  with  An  =  0.01  over  a  distance  D  =  lOum,  B  has  the  value  of 

1089. 

2 
It  might  be  mentioned  that  in  the  radio  ionospheric  case,  (An) 

9 
may  not  be  small  but  may  approach  nQ.   In  such  a  situation,  the  problem 

.,  11   . 
is  more  complicated,  and  Budden  discusses  it  in  some  detail.    As 

he  explains  it,  the  field  component  Ex  imparts  vertical  motions 
to  the  electrons  and  if  E  is  large,  then  within  one  cycle  any  one 
electron  sees  a  different  E  at  different  levels.   In  other  words, 
the  electron  does  not  see  a  constant  force.  Hence  its  motion  does  not 
vary  harmonically  with  time  and  is  therefore  anharmonic.   Because 
of  this,  some  energy  goes  into  harmonics,  and  the  original  wave  is 
attenuated.  As  a  result,  the  phase  factor  encountered  upon  reflec- 
tion (or  total  internal  refraction)  must  have  an  imaginary  component 
to  account  for  this  loss  in  the  TM  case.   Further,  the  phase  change 
is  less  than  it/2.  Maximum  attenuation  (i.e.,  conversion  to  harmonics) 
occurs  when  B  ^  0.45,  in  which  case  only  about  72%  of  the  incident 
energy  is  reflected  by  the  ionosphere. 

For  the  optical  waveguide  case,  B  is  very  large;  hence  a 
negligible  loss  should  occur  per  reflection  (i.e.,  per  total  internal 
refraction).   However,  unlike  radio  waves  which  reflect  once  from 
the  ionosphere,  light   in  a  waveguide  might  be  reflected  many  times 
during  its  passage  within  the  guide.   Perhaps,  for  certain  index 
profiles,  enough  energy  can  be  converted  to  higher  harmonics  to 
appreciably  attenuate  the  TM  wave.  And  perhaps,  with  the  proper 
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phase  matching,  non-linear  effects  such  as  second-harmonic  generation 
or  parametric  upconversion  might  take  place.   As  Budden  points  out, 
the  problem  of  TM  attenuation  is  very  difficult,  requiring  numerical 
computations.   Since  from  Eq.  (16)  or  (17)  [and  Eq.  III-(6)],  0(x) 
is  very  small  for  small  An,  then,  for  a  propagation  length  of  a  few 

millimeters,  the  number  of  total  internal  refractions  will  be  quite 

2 
small.   Further,  since  B  is  so  large  (on  the  order  of  10  or  more) 

for  such  profiles,  a  negligible  loss  per  refraction  should  occur. 

Thus,  no  significant  attenuation  should  occur  for  a  TM  wave  in  a  slowly 

varying  index  optical  waveguide  which  is  only  a  few  millimeters  in 

length  (such  as  the  zinc  sulfide  films  used  in  this  research  project). 

It  may  be  that,  waveguides  with  steeper  profiles  and  longer  lengths  will 

produce  significant  attenuation  of  a  TM  wave. 

It  is  now  possible  to  write  down  the  solutions  for  E   and  H^. 

For  the  region  x  =  0,  the  solutions  are  given  by  Eqs.  III-(46),  C-(224), 

and  C-(227).   For  the  region  (-W  +  A)  =  x  =  0,  the  solutions  are  the 

WKBJ  solutions  given  by  Eqs.  III-(47),  C-(225),  and  C-(228).   For  the 

region  (-W  -  A)  =  x  =  (-W  +  A),  the  solutions  involve  the  Airy  integral 

function: 

+fo  (*)  =    a  i.  (3  K^  At  (5)  iv- (6i) 


The  exact  form  of  E   and  E   is  given  in  Eqs.  D-(151)  and  D-(152). 

OX        OZ  -in' 

Finally,  for  the  region  x  =  (-W  -  A),  the  WKBJ  exponentially  decaying 
solution  is  y 

S  ~^\V       ^ie  IV_(62) 
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where  A.  is  defined  in  Eq .  D-(144).   The  exact  form  of  F   and  F   in 
4  ox      oz 

this  region  is  given  in  Eqs.  D-(154)  and  D-(155). 

From  the  analysis  for  TE  and  TM  modes  given  above  and  in  Appendix 
D,  all  the  results  obtained  from  the  quasi-geometric  ray-optics  approach 
can  be  verified  from  considerably  more  rigorous  mathematics.   Thus  it  can 
be  concluded  that  a  practical  and  accurate  understanding  of  the  phenomenon 
of  total  internal  refraction  in  optical  waveguides  can  be  obtained  by 
using  physical  insight  arguments  without  the  necessity  of  finding  the 
actual  solutions  for  E  and  H.   Of  course,  the  conditions  under  which  such 
arguments  are  valid  must  be  determined  from  the  more  rigorous  approach. 


CHAPTER  V 
GUTDING  OF  LIGHT  IN  EPITAXIAL  ZINC  SULFIDE  FILMS 

Introduction 

Cubic  zinc  sulfide  (also  called  B-ZnS,  zincblende,  or,  for  naturally 

1  6 

occurring  crystals,  zincblende  or  sphalerite)  has  a  bulk  refractive  index 

n  =  2.35  at  X   =  0.6328pm.   The  epitaxial  zinc  sulfide  films  fabricated 

by  the  Central  Research  Laboratories  of  Itek  Corporation  are  grown  on 

17  18 
gallium  arsenide  (GaAs)  substrates.   '     The  bulk  refractive  index  of 

19  ~ 
GaAs  is  complex:     n_  =  n_  +  ±Kn   =  3.87  +  0.087i  at  X   =  0.6328pm. 

(j        Lr  L>  O 

[Although  the  imaginary  part  k  of  the  index  is  only  about  2% 

G 

20 
of  the  real  part  n  ,  the  absorption  coefficient    a   is  approximately 

4   -1 
1.7  x  10  cm      ,    so,  for  X   =  0. 6328pm,  no  light  propagation  can  take 

place  in  GaAs  crystals  of  length  L  greater  than  a  few  micrometers.] 

I—  i  21 

Since  k   <<  n  ,  then  |n  |  ~  n  ,  and  the  reflection  coefficient   R  has 

little  dependence  on  k  .   Further,  since  n   >  n  ,  then  light  propagation 

in  an  epitaxial  zinc  sulfide  film  on  a  gallium  arsenide  substrate  cannot 

take  place  solely  by  means  of  total  internal  reflection.   Instead,  either 

partial  internal  reflection  or,  possibly  (as  discussed  below),  total 

internal  refraction  must  be  the  mechanism  for  guiding  light  at  or  near 

the  ZnS-GaAs  boundary. 

Total  Internal  Refraction 

Cubic  GaAs  has  a  lattice  spacing  of  5.64  A;   cubic  ZnS  has  a  natural 

spacing  of  5.41  A  (Ref.  17).   If  there  is  a  perfect  lattice  match  at 
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the  ZnS-f ilm/GaAs-substrate  boundary,  then  the  spacing  of  the  ZnS  film 

o  o 

lattice  will  change  from  5.64  A  to  a  more  closely  packed  5.41  A, 
a  difference  of  4.25%.   It  seems  possible  that,  because  of  the  greater 
ZnS  lattice  spacing  at  the  ZnS/GaAs  boundary,  a  lower  ZnS  density  as 
well  as  a  lower  ZnS  refractive  index  should  result.   Then,  as  the 
epitaxial  ZnS  lattice  changes  from  5.64  A  to  5.41  A,  an  index  gradient 
could  exist  within  the  guide,  with  a  higher  value  of  the  index  n  at 
distances  increasingly  farther  from  the  ZnS/GaAs  boundary.   Thus,  light 
could  propagate  in  the  ZnS  film  by  means  of  total  internal  refraction  at 
some  level  above  the  ZnS/GaAs  boundary. 

There  is  a  second  possible  mechanism  which  might  give  rise  to  a 
positive  index  gradient  for  n   in  addition  to  the  lattice  mismatch.   The 
coefficient  of  linear  thermal  expansion  (at  200°C)  of  zinc  sulfide  is 
6.3  x  10   /°C:  for  gallium  arsenide  the  value  is  6.5  x  10   /°C  (Ref.  22). 
After  crystal  film  growth  at  typically  400°C,  the  ZnS  and  GaAs  must  cool 
down  simultaneously.   In  doing  so,  the  smaller  thermal  expansion  of  the 
ZnS  film  causes  a  small  bending  of  the  ZnS/GaAs  combination,  such  that  the 
ZnS  film  is  on  the  inner  diameter.   As  a  result,  the  density  of  the  ZnS 
film  is  perhaps  slightly  higher  at  the  ZnS/superstrate  boundary.   Such 
a  bending  was  observed  previously  by  scientists  at  Itek  (unpublished), 
and  a  Kfocal  length"  f  of  approximately  one  meter  was  observed  for  a 
collimated  beam  of  laser  light.   (To  the  eye,  the  surface  appears  perfectly 
flat.)   Thus  the  radius  of  curvature  of  the  bending  is  2f.   For  a  ZnS 
film  thickness  of  W  =  20um  or  so  (a  typical  value  of  the  Itek-grown  films 
used  for  the  light  guiding  experiments),  then  only  a  small  change  in  the 
radius  of  curvature  2f  is  present  between  the  top  and  bottom  of  the  film. 
For  a  crystal  length  dimension  L  at  the  top  surface,  then  the  angle 
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subtended  by  the  top  surface  of  the  crystal  is  L  /4irf:  for  the  bottom 

surface,  the  same  angle  is  subtended  and  can  be  written  as  L  /2ir(2f  +  W) , 

B 

where  L   is  the  length  of  the  film  along  the  bottom  of  the  ZnS  film. 

Equating  the  above  two  expressions  leads  to  L  =  L  (1  +  W/2f).   Since 

b    r 

W/2f  ^10   ,  it  is  reasonable  that  any  index  gradient  due  to  a 
differential  thermal  expansion  should  be  negligible  compared  to  an  index 
gradient  caused  by  the  lattice  mismatch. 

The  question,  then,  is  how  far  into  the  ZnS  film  does  the  lattice 

23  24 
mismatch  have  any  effects.   Experiments  by  other  authors   '    with  zinc 

sulfide  films  deposited  on  glass  and  fused  silica  substrates  have 

confirmed  the  presence  of  index  inhomogeneities  extending  anywhere  from 

o  o 

100  A  to  800  A  (i.e.,  from  20  to  150  lattice  spacings)  away  from  the 
substrate/film  boundary.   Of  course,  since  the  substrates  were  amorphous, 
the  zinc  and  sulfur  atoms  had  complete  freedom  to  align  themselves:  in 
fact,  the  ZnS  films  were  nearly  amorphous  themselves,  showing  a 
minimum  of  crystalline  structure.   On  the  other  hand,  the  epitaxial 
growth  of  zinc  sulfide  films  on  gallium  arsenide  substrates  requires 
that  the  ZnS  molecules  align  themselves  according  to  the  cubic- 
structured  GaAs  lattice.   Hence,  the  possibility  exists  that  any 
inhomogeneity  due  to  a  lattice  mismatch  could  extend  more  than  a  1000  A 
into  the  ZnS  film. 

Partial  Internal  Reflection 
Partial  internal  reflection  at  the  ZnS/GaAs  boundary  is  another 

possible  mechanism  for  guiding  of  light  in  a  zinc  sulfide  thin  film. 

25 
Such  a  light  guide  is  a  leaky  one   since  some  energy  is  lost  upon 

25 
reflection  by  waves  at  the  ZnS/GaAs  boundary.   Ulrich  and  Prettl 
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discuss  leaky  light  guides  in  detail  and  arrive  at  the  following 
relationships  for  a  given  mode  number  m: 

^m  '-       7    a      aw \3  (l*    «A,j  V-(3) 


TH  =    nfiraAft (gUii!  /  \ 


V-(4) 


where 


Jt.W.j-    *»W+     c    =.^y/a       (-TH-mafe)      ? 


3  "*^  '  V-(5) 


V-(6) 


Eq .  (2)  can  be  considered  as  a  characteristic  waveguide  equation  for 
leaky  light  guides.   The  quantities  n  ,  n  ,  n  and  W  are  defined  in 
Figure  II-l  except  that,  for  leaky  guides, 

>  n 


m5>  n3  >  7W>  nfc  v_(7) 


A 

The  quantity  B   is  similar  to  the  wavevector  component  along  the  z-axis 

m 

[from  Eq .  II-(22)],  only  now  a  provision  is  made  for  attenuation  of  the 
light  by  including  the  term  -lot  . 
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25 
The  quantity  W   is  the  "equivalent  thickness"   of  the  guide  and 

is  the  result  of  the  fact  that  the  totally  internally  reflected  wave  at 

the  n  -  n  boundary  penetrates  slightly  into  the  n  region.   For  A  = 

0.6328um  and  W  =  20um,  then  k  W  =  200:  for  n  =  n,  =  2.35  and  n  =1.63 

o  g    Z  a 

(the  largest  index  of  the  superstrate  materials  used  in  waveguiding 

2     2  -1/2 
experiments),  then  (n   -  n   )     =  0.6  <<  k  W.   Fence,  for  all  practical 
g     a  o 

purposes,  W  can  be  used  instead  of  W   in  Fqs.  (1)  -  (6)  for  both  TE  and 

t  t-  eq 

TM  modes  in  Itek-grown  epitaxial  zinc  sulfide  light  guides. 

25 
Eqs.  (1)  -  (4)  were  derived   under  the  assumption  that  the  imaginary 

part  k   of  the  refractive  index  n  is  much  less  than  the  real  part  n  ,  i.e., 

K   <<  n  and  In  I  ^  n  .   As  explained  in  the  introduction  to  this  chapter, 

s     s      '  s '     s 

the  values  of  n  and  k   satisfy  this  assumption.   Using  the  value 

n  =  n„  =  3.87  together  with  A   =  0.6328um  and  W  =  20pm,  then,  from 
s    G  o 

Eqs.  (3)  and  (4) 

7T        ,a  ,  -3 


V-(8) 


-C    ^Cw  +  ir^x/o^m^1)  v"(9) 

Since  the  samples  used  for  light  guiding  experiments  had  lengths  of  a 
few  millimeters,  negligible  attenuation  should  occur  for  the  lower- 
order  modes  (m  <  5). 

Ignoring  the  term  lot   in  Eq.  (1),  then  from  Eq .  II- (22),  r\      can  be 
written  as 


1 


V-(10) 


m~    "J  ^"   "  ™ 
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where  0   is  the  value  of  0  in  Figure  II-l  for  a  given  mode  m.   Combining 

m  no  o 

Eqs.    (2)    and    (10), 


n3(i-^emy  JUL   &± 
an3     $tVy 


S^9^=     I    -        -A (^-hJ 


Using  the  values  of  X    ,    n  =11,,,  and  W  used  above, 
o'   g    Z' 


V-(ll) 


n'W3"  v-(i2) 


Se*&„  =  I"   O^*^  Xm-rlj         v_(13) 

It  may  be  seen  that,  to  a  very  good  approximation,  0   is  very  nearly 

m 

equal  to  (but  always  slightly  less  than)  tt/2  —  propagation  straight 
down  the  guide.   This  is  true  for  all  but  the  highest  order  modes, 
which,  from  Eqs.  (8)  and  (9),  are  strongly  attenuated  anyway. 


CHAPTER  VI 
THE  LINEAR  ELECTRO-OPTIC  EFFECT  OF  EPITAXIAL  ZINC  SULFIDE  FILMS 

The  Principal  Axes 
The  optical  properties  of  a  crystal  can  be  described  in  terms  of 
an  optical  indicatrix,  a  refractive  index  ellipsoidal  surface.   Since 

zinc  sulfide  is  isotropic,  its  indicatrix  can  be  represented  by  the 

26 
following  equation: 

a    -  I  vi-(i) 

which  is  the  equation  of  a  sphere.   Here  n  =  n „,  the  bulk  refractive 
index.   Upon  application  of  an  external  electric  f ield  "e  (not  to  be  confused 
with  the  electromagnetic  field  E  of  a  light  wave)  for  the  linear  electro- 

•  u  26 

optic  effect,  then  the  indicatrix  becomes 

a       s      2 


H1L   ^  3r/excj2fe^  +  e2^)--| 

a  ^ '  v      y  J  vi-(2) 


n 


Here  r,   is  the  electro-optic  coefficient  of  zinc  sulfide. 
41 

The  epitaxial  films  grown  by  Itek  Corporation's  Central  Research 

A 

Laboratories  are  oriented  in  the  (100)  direction,  defined  as  the  x 
direction  in  Figure  VI-1.   On  top  of  this  film  an  electrode  is  deposited. 

A 

Thus,  the  applied  field  across  the  film  is  in  the  x  direction  only,  i.e., 
e  =  e  =  0.   Eq .  (2)  then  reduces  to 

y   z 
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Figure  VI-2.   Crystal  orientation  of  epitaxial  zinc  sulfide  films 
grown  by  the  Central  Research  Laboratories  of  Itek 

A 

Corporation.   The  principal  axes  are  given  by  the  x' , 
y',  and  z'  directions  which  are  also  the  normals  to 
the  natural  cleavage  planes. 
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XV+2    +  ^,e<a2=l  vi-o) 


By  a  proper  rotation  of  the  index  ellipsoid,  a  new  set  of  principal 


axes  (x1,  y',  z')  can  be  obtained  such  that 

,2  Mi3  2/a 


x; J 

n-a  +   ,a  *   ^    '  vi-(4) 

j 


Comparing  Eqs.  (3)  and  (4)  it  may  be  seen  that  the  x  axis  remains 
unchanged,  i.e., 


A  /   A 

X  =  X  VI-(5) 


and 


Y\      s.  y\  VI-(6) 


The  problem  thus  reduces  to  calculating  the  angle  0,   through  which  the 

A         A  A  A 

y  and  z  axes  must  be  rotated  to  y'  and  z'.   As  shown  in  Appendix  E, 

A.     *  A 

ft   =  tt/4.   The  new  axes  x' ,  y',  and  z'  are  shown  in  Figure  VI-1.   They 
coincide  with  the  normals  to  the  natural  cleavage  planes  of  cubic  zinc 
sulfide. 

In  the  course  of  obtaining  the  principal  axes,  the  dependence  of 
the  principal  refractive  indexes  n' ,  n' ,  and  n'  on  the  applied  field  e^ 
are  also  obtained  in  Appendix  E.   The  results  to  first  order  are 

^"a  V!-(7) 


n,'=  n  -  n  3  n  e«£ 
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VI-(8) 


l^=  ^o+  ^o%(^  vi-(9) 


3 
where  n   rne  /2  <<  n  .   It  may  be  seen  that  n'  and  n'  have  a  linear 
o   41  x       o  y      z 

dependence  on  the  applied  field  e  . 

The  Transverse  Electro-Optic  Effect 

For  propagation  in  either  the  z  or  z'  direction  [the  (001)  and  (Oil) 

directions,  respectively],  the  light  has  a  k  -  vector  component 

6  =  k  n   sinO  which  is  perpendicular  (i.e.,  transverse)  to  the  applied 
o  o 

electric  field  e.   This  results  in  the  transverse  electro-optic  effect. 
[This  is  not  to  be  confused  with  the  transverse  electric  (TE)  mode  of 
propagation. ]   The  light  also  has  a  k-vector  component  k  n  cos©  which 
is  parallel  to  e,  resulting  in  the  longitudinal  electro-optic  effect. 

In  order  to  get  intensity  modulation  (as  opposed  to  pure  phase 
modulation)  when  the  ZnS  guide  is  placed  between  crossed  polarizers,  it 
is  necessary  that  the  polarization  of  the  incident  light  be  altered.   In 
order  to  determine  if  changes  in  polarization  will  occur,  all  incident 
polarization  components  must  first  be  resolved  along  the  principal  axes 
(so  long  as  the  axes  are  perpendicular  to  the  net  direction  of  light 
propagation) .   If  there  are  two  polarization  components  each  parallel  to 
a  different  principal  axis  which  is  perpendicular  to  the  net  direction  of 
travel,  then  each  polarization  component  will  "see"  a  different  refractive 
index  if  an  electric  field  is  applied  to  the  crystal.   As  a  result,  each 
component  will  have  a  different  phase  retardation,  and  the  light  will 
become  elliptically  polarized  as  it  propagates  down  the  guide.   The  net 
effect  is  to  rotate  the  plane  of  polarization  of  the  incident  light. 
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A  combination  transverse  and  longitudinal  electro-optic  effect  can 
in  general  be  very  complicated.  However,  for  optical  waveguides,  there 
are  several  factors  which  tend  to  reduce  the  complexity. 

First  of  all,  the  angle  0  will  be  very  nearly  90°  for  epitaxial  zinc 
sulfide  films  which  guide  light  at  or  near  the  ZnS/GaAs  boundary  either 
by  partial  internal  reflection  or  by  total  internal  refraction.   In  the 
former  case,  this  is  obvious  from  Eq.  V-(13).   Since  the  longitudinal 
component  depends  on  the  cosine  of  0,  it  will  be  negligible  for  0  ^  tt/2. 
In  the  latter  case,  the  change  in  index  due  to  a  lattice  mismatch  is 
probably  no  more  than  4.25%,  i.e.,  An  =  0.0425n  .   From  Eqs.  IV-(16) 
and  D-(73),  sin0(O)  cannot  be  any  less  than  n^l (nz  +  An)  =  0.96,  i.e., 
0(0)  cannot  be  any  less  than  74°.   Of  course,  0(x)  changes  from  0(0)  to 
0(-W)  =  90°.   Hence, very  roughly,  the  minimum  "average''  value  of  0(x) 
would  be  82°.   Since  sin  82°=  0.99  and  cos  82°  =  0.14,  it  may  be  seen 
that  the  longitudinal  component  (from  cos  82°)  is,  at  the  very  most, 
only  14%  of  the  transverse  component. 

Even  more  significantly,  the  voltage  required  to  produce  a 
longitudinal  electro-optic  effect  in  a  thin-film  waveguide  is  orders  of 
magnitude  greater  than  the  voltage  necessary  to  produce  a  transverse 
effect.   For  a  pure  longitudinal  effect,  the  half-wave  voltage  (defined 

as  the  applied  voltage  required  to  change  the  plane  of  polarization 

,26 
of  the  incident  light  by  90  )  is  given  by 

^o 

'wa'S     £f),?  rr  vi-(io) 


■    26 
For  a  pure  transverse  effect,  the  half-wave  voltage  is 
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V,/a) 


W 


3 


r««      n0V    L  vi- (ii) 


-  ^_W_  (Ya/3l)l^  vi-(i2) 


Using  A   =  0.6328um,  n  =  n_  =  2.35,  and  r.  n  =   2.07  x  10    cm/volt 
°     o  o    Z  41 

(Ref.  26),  then  (V,  ._),     =  1.2  x  10  volts.   For  W  =  20um  and  L  =  1cm, 
'        X/2   Long 

then  (V,  ,„)„     =  10~3  (V, /0)T     =  47  volts.   Thus,  with  an  applied 
X/2   Trans  X/2   Long 

voltage  of  50  volts  or  so,  a  large  rotation  of  the  polarization  will  occur 
as  a  result  of  the  transverse  electro-optic  effect:  a  negligible  effect 
will  be  observed  for  the  longitudinal  effect. 

Further,  since  optical  waveguides  are  quite  thin,  dielectric  break- 
down will  occur  before  a  very  large  voltage  can  be  applied  to  the  thin- 
film  material.   For  zinc  sulfide,  the  breakdown  voltage  for  practical 
purposes  is  approximately  10  volt/cm  =  10  volt/um  (Ref.  27).   Thus,  a 
20um-thick  zinc  sulfide  film  would  break  down  long  before  any  significant 
longitudinal  effect  occurred. 

Thus  it  can  be  concluded  that  the  transverse  electro-optic  effect 
is  the  primary  mechanism  for  electro-optic  modulation  in  zinc  sulfide 
optical  waveguides.   It  might  be  pointed  out  that,  since  6  =  ^Qnz   sinQ 
is  a  constant  for  homogeneous  as  well  as  inhomogeneous  waveguides  (as 
explained  in  Chapters  II-IV) ,  the  transverse  component  of  the  k-vector 
is  thus  a  constant  for  either  case. 

Intensity  Modulation  for  Propagation  in  the  (011)  Direction 

1.   TE  Modes 

For  a  pure  TE  mode,  the  electric  vector  E  is  polarized  parallel  to 


the  y   axis ,  i.e., 
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M     c   I 


E        -         L4   ^ u  VI- (13) 


A  /  ,- 


~  ij  VI-(14) 


where  E  =  E  (x)  el(wt    3z) .   The  wave  will  "see"  an  index  n^  =  nQ  =  nz 

3 
with  zero  applied  field  e  ,  and  an  index  n'  =  nz  ~  nz  r4lex/'2  with  an 

applied  field.   Thus,  phase  modulation  can  occur  for  the  TE  case.   However 

since  the  light  is  polarized  parallel  to  only  one  principal  axis,  then 

for  any  value  of  e   the  E-field  sees  only  one  index  of  refraction  at  any 

one  time.   Hence  the  polarization  will  not  be  altered,  and  no  intensity 

modulation  is  possible. 

2.   TM  Modes 

For  a  pure  TM  mode,  the  electric  vector  E  can  be  written  as 

£   a  X   C   +-  2  t^ 
TM        X  VI- (15) 

=  X      Ez^&^-i     £c<jS&-tM  VI- (16) 

Since  the  propagation  is  in  the  z'  direction,  it  may  be  seen  that  there 
is  only  one  polarization  component  along  a  principal  axis  perpendicular 
to  the  net  direction  of  travel.   Hence,  the  polarization  will  not  be 
altered,  so  no  intensity  modulation  is  possible;  phase  modulation 
also  is  not  possible. 
3.   TE  +  TM  Modes 

For  TE  +  TM  incident  light,  the  E-vector  can  be  written  as 

Ere,7M  s     Etf  +    Etm  vi_(i?) 
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-  U  £*   +  ^  Bu    +  *  fczJ  VI'(18) 


In  this  situation,  with  an  applied  field  e^  the  speed  of  propagation 

for  E'  is  different  from  the  speed  of  propagation  for  E*.  Thus  after 
x                                             y 

travelling  a  length  L  through  the  crystal,  unequal  phase  shifts  6^  and 
6'  are  introduced: 

E±=    fe^^TA,  )  ^  VI_(20) 

E„     ~         ^6    e  VI-(21) 


where 


VI-(22) 


and 


%l  i  ft   L 

dM   ^Te  vi- (23) 


PTW      °   x       r/^  VI- (24) 


?Te-Ana^&Te- 


The  values  of  0^  and  9_,  are  to  be  determined  from  the  appropriate 

TE       J  M. 

characteristic  waveguide  equation,  depending  on  the  mechanism  for  light 
guiding. 
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It  is  now  assumed  that  the  zinc  sulfide  light  guide  is  placed 
between  crossed  polarizers  (where  the  input  polarizer  is  oriented  so  as 
to  pass  a  TE  +  TM  wave).   Upon  application  of  a  field  e^,    the  E^  and  E^ 
components,  each  at  angle  45°  to  the  analyzer  axis,  have  components 
-E'//2  and  +E'//2  along  the  analyzer  axis.   The  amplitude  EA  of  the 
light  passed  by  the  analyzer  is  thus 


£A   =  (-£„'  4  F'j/iJ 


VI-(26) 


where  the  z'  component  has  been  neglected  since  0^  *   ir/2,  as  explained 
above.   From  Eqs.  (20)  and  (21),  Eq .  (26)  becomes 

o  x  l  o  u    \   t  to  L 


VI-(27) 


The  mean  (i.e.,  time-averaged)  intensity  I  transmitted  by  the  analyzer 


■£    ~~    ^A  A   '   **  VI- (28) 


From  Eq.  (27)  this  results  in 


I  m  m  ^       ,  M   ^        VI-(29) 


-t-e 

Defining 


6       '  VI-(30) 
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then 


1=  (iji)l  s^\M  +  h3^&m  ^o  &-lj)] 


VI-(31) 


Using  Eqs.  (22)  -  (25)  together  with  Eqs.  (7)  and  (8),  then  the  phase 
shifts  can  be  written  as 


3 


Thus, 


where 


%  ~  K^o1  C^Qm-^^Tt 


Further,  since 


^   t^^VL^ 


VI-(32) 


VI-(33) 


sx'-s;>  s.+W  ^.Le./a 


VI-(35) 


e    =  v/vu 

^       /  VI-(36) 


where  V  is  the  voltage  across  the  crystal,  then 

*      Jj      °  tl  I  VI-(37) 


3 

f/     /  VI-C38) 
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If   the  approximation  is  made   that 


<9'rK  ~    V'Te°Z    "A  VI-(39) 


which  is  valid  for  epitaxial  zinc  sulfide  light  guides  (as  explained 
above}  then  Eqs.  (31)  and  (38)  become 


1=   X0    w^  (hx'-l') 


VI-(AO) 


and 


3 


^-K)-v<r«nrM 


This  is  the  same  result  obtained  for  the  bulk  transverse  electro-optic 
effect,    and  intensity  modulation  will  thus  occur  as  the  applied 
voltage  is  varied  from  zero  to  (Vi  ^Trans' 

Tt  should  be  mentioned  that  if  there  is  a  dielectric  layer  isolating 
the  metal  electrode  from  the  zinc  sulfide  film,  then  the  voltage  V  in 
Eq.  (41)  will  not  be  the  total  voltage  V  ,  but  only  the  voltage  across 
the  ZnS  film.   The  relationship  between  V  and  V   can  be  obtained  from 
elementary  physics  since  the  layer  and  ZnS  film  are  effectively  two 
parallel-plate  capacitors  in  series.   The  result  is 

V-        W€i  V 

V  ~   v~  VI-(42) 

where  e,  and  W,  are  the  low-frequency  dielectric  constant  and  thickness, 
d      d 

respectively, of  the  layer;  e      is  the  low-frequency  dielectric  constant 
of  zinc  sulfide  and  has  the  value  8.37  +  0.8  picofarad/meter  (Ref.  28). 
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Intensity  Modulation  for  Propagation  in  the  (001)  Direction 
1.   TE  Modes 

For  a  pure  TE  mode,  the  electric  vector  E  is  polarized  to  the  y 
axis,  i.e., 

15       J  VI-(43) 

°  VI- (44) 


=(<J'->/^ 


VI-(45) 


The  polarization  components  of  the  incident  light  must  always  be  at  right 

A  A 

angles  to  the  direction  of  propagation.   Since  the  y   and  z   axes  are 

A 

not  perpendicular  to  the  net  direction  of  propagation  for  the  z-propagation 

case,  then  the  method  of  resolving  the  incident  polarization  along  principal 

axes  is  not  a  proper  way  of  examining  the  problem.   Instead,  the  effective 

contributions  of  the  axes  can  be  evaluated  by  resolving  these  axes  along 

the  axis  formed  by  the  incident  polarization.   In  this  situation,  any 

change  An'  due  to  the  y'  axis  is  identically  compensated  by  any  change 

An'  =  -An'  [from  Eqs.  E-(17)  and  E-(18)]  due  to  the  z'  axis.   Hence  the  TE 
z      y 

wave  "sees"  a  constant  index  n  =  n„.   The  net  effect  is  no  change  at  all 

o    L 

(not  even  phase  modulation  when  a  field  is  applied  to  the  crystal)  .  In 

A 

other  words,  for  propagation  in  the  z  direction,  at  least  for  TE  modes,  the 
zinc  sulfide  appears  isotropic,  with  or  without  an  applied  electric  field. 
2.   TM  Modes 

For  a  pure  TM  mode,  the  electric  vector  E  can  be  written  as 

VI-(46) 
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A    /  /  A/  /  A/         / 


TH 


VI-(47) 


VI-(48) 


Since   0        ^  tt/2,    then   Eq .    (48)    reduces   to 

^V/H  *  VI- (49) 


which  is  the  same  as  Eq .  (16).   Thus  no  intensity  modulation  is  possible 

(and   since  n'  =  n   is  independent  of  the  applied  field  e  ,  no  phase 
'        x    o  x 

modulation  is  possible  either).   Even  if  0TM  is  not  quite  equal  to  tt/2, 
there  will  still  be  no  net  phase  change  from  the  y   and  z   terms  for 
the  same  reasons  given  for  TE  modes  propagating  in  the  z  direction. 

3.   TE  +  TM  Modes 

For  TE  +  TM  incident  light,  the  E-vector  can  be  written  as 


£rz  +  in~     STir  +    BtH 


VI-(50) 


=   x   tx   rye    +■  -z  c2 


^  VI-(51) 


-  K 


VI-(52) 


Again,    since   0        ^  tt/2,    then 


TM 


^TFrTTH  '  VO  /  VI- (53) 
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Again,  the  same  arguments  used  for  TE  and  for  TM  modes  lead  to  the 

A  A 

conclusion  that  the  y'  and  z'  term  will  produce  no  net  effect.   Thus, 
even  though  from  the  appearance  of  Eq .  (52)  it  might  seem  that  there 
should  be  some  modulation  (either  intensity  or  phase), such  is  not 
the  case.   Hence  the  z  direction  is  an  optic  axis  of  cubic  zinc  sulfide 
since  the  refractive  index  is  the  same  for  TE,  TM,  or  TE  +  TM  polar- 
izations . 


CHAPTER  VII 
EXPERIMENTAL  PROCEDURE 

Sample  Preparation 

Most  of  the  Itek-grown  samples  have  a  circularly  shaped  zinc 
sulfide  film  which  does  not  extend  all  the  way  to  the  edge  of  the 
gallium  arsenide  substrates.   Thus,  in  order  to  propagate  light  in  the 
z'  direction  [the  (Oil)  direction  shown  in  Figure  VI-1] ,  it  is  necessary 
to  cleave  each  sample  twice  (indicated  in  Figure  VII-1) ,  thereby 
exposing  one  edge  for  input  coupling  of  the  light  and  another  edge 
at  the  other  end  for  the  output.   This  is  accomplished  by  placing  the 
sharp  edge  of  a  one-sided  razor  blade  in  the  proper  orientation  on  the 
GaAs  side  of  the  ZnS-on-GaAs  sample  and  then  tapping  the  blade  slightly. 
In  this  way  a  microscopically  smooth  edge  of  the  ZnS  film  is  obtained. 

For  propagation  in  the  z  direction  [the  (001)  direction],  it  is 
necessary  to  slice  the  input  and  output  ends  of  the  ZnS-on-GaAs  samples 
(using  a  diamond  saw)  at  a  45   angle  to  the  normals  to  the  cleavage 
planes.   The  sample  is  then  sandwiched  between  two  glass  plates  and 
subsequently  fine  polished.   The  glass-plate  sandwich  is  used  because 
the  first  polishing  attempts  without  the  plates  resulted  in  a  bevelled 
ZnS  edge.   As  a  result,  light  incident  on  the  film  edge  was  refracted 
into  the  GaAs  substrate,  and  propagation  down  the  length  of  the  thin- 
film  guide  could  not  be  achieved.   Schematic  diagrams  of  the  resulting 
samples  for  z  propagation  are  shown  in  Figure  VII-2. 
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Figure  VII-1.   Top  view  of  cleaved-edge  sample  for  light  propagation 
in  the  z'  direction  [the  (Oil)  direction]. 
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Figure  VII-2.   Top  view  of  sliced-edge  sample  for  light  propagation 
in  the  z  direction  [the  (001)  direction]. 
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For  electrode  preparation,  the  top  surface  of  the  sample  is 
masked  off  with  Scotch  tape,  leaving  exposed  the  desired  area  for  the 
electrode.   Then  electrically  conductive  silver  paint  (Epoxy  Products  Co., 
New  Haven,  Conn.)  is  applied  to  this  area  and  allowed  to  dry,  after 
which  the  tape  is  removed.   The  sample  is  then  mounted  as  shown  in 
Figure  VII-3.   A  glass  substrate  is  used  as  an  electrically  insulating 
mechanical  base.   Half  of  the  top  surface  of  the  plate  is  covered  with 
electrically  conducting  silver  paint,  and  the  GaAs  side  of  the  ZnS-on- 
GaAs  sample  is  placed  on  one  portion  of  this  paint  before  it  dries. 
(Gallium  arsenide  has  a  much  lower  resistivity  than  zinc  sulfide,  and 
therefore  it  acts  as  the  bottom  electrode  to  the  ZnS  light  modulator.) 
To  the  other  portion  of  this  paint  an  electrical  lead  is  affixed  with 
Dupont  Duco  cement  (since  the  silver  paint  does  not  provide  much 
mechanical  strength) . 

On  the  other  half  of  the  glass  substrate  more  silver  paint  is 
deposited  (electrically  isolated  from  the  first  half).   Then  an  electrical 
lead  is  affixed  to  this  paint  and  to  the  top  electrode  on  the  ZnS 
film  using  some  more  paint.   Finally,  another  electrical  lead  is 
bonded  with  Duco  cement  to  the  silver  paint  on  the  glass  base.   The 
sample  is  thus  mechanically  mounted  for  waveguiding  experiments  and 
electrically  connected  for  electro-optic  measurements. 

Experimental  Set-up 
A  helium-neon  laser  is  used  as  the  primary  light  source  for  the 
experiments.   The  output  beam  of  light  (A   =  0.6328pm)  is  linearly 
polarized  so  that  use  of  a  half-wave  plate  and  an  angularly  calibrated 
linear  polarizer  allows  easy  selection  of  TE,  TM,  or  TF  +  TM  light.   The 
light  is  coupled  into  the  guide  by  means  of  a  20X  microscope  objective. 
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Figure  VIT-3.   End  view  of  mechanical  and  electrical  mounting  for 
waveguiding  and  electro-optic  investigation  of  zinc 
sulfide  thin-film  samples. 
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At  the  output  end,  a  microscope,  consisting  of  a  20X  or  44X  objective 
and  a  10X  eyepiece,  is  used.   This  experimental  set-up  is  shown 
schematically  in  Figure  VII-4.   The  input  objective,  the  ZnS-on-GaAs 
sample  mount,  and  the  output  microscope  are  mounted  independently  on 
translational  stages  each  capable  of  precision  motion  in  three  dimensions. 
The  sample  mount  has  an  angular  degree  of  freedom  also.   For  the  electro- 
optic  modulation  experiments  an  analyzer  (with  transmission  axis  90°  to 
the  input  polarizer  axis)  is  placed  between  the  ZnS  sample  and  the  output 


microscope. 
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Figure  VII-4.   Side  view  of  experimental  set-up  for  waveguiding  experi- 
ments (not  drawn  to  scale) . 


CHAPTER  VIII 
EXPERIMENTAL  RESULTS  AND  OBSERVATIONS 

Light  Guiding 
Experiments  with  light  guiding  in  epitaxial  zinc  sulfide  films  were 
performed  to  first  determine  whether  partial  internal  reflection 
or  total  internal  refraction  was  the  primarv  mechanism  for  light  guiding 
at  or  near  the  ZnS/GaAs  boundary.   Preliminary  experiments  had  resulted 
in  an  output  intensity  pattern  with  a  maximum  nearer  the  ZnS/superstrate 
boundary,  suggesting,  perhaps,  total  internal  refraction.   This  is 
shown  in  Figure  VIII-2:  Figure  VIII-1  is  a  photograph  of  the  entire  guide 
illuminated  in  transmission  by  white  light.   Subsequent  experiments  with 
several  different  samples  showed  that  the  maximum  of  the  distribution 
could  be  shifted  to  the  center  of  the  guide  (Figures  VIII-3  and  VIII-4) , 
and  to  the  bottom  of  the  guide  (Figure  VIII-5),  all  for  the  same  mode. 
(This  is  essentially  an  m  =  0  mode,  although  the  beginning  of  a  conversion 
to  an  m  =  1  mode  can  be  seen,  especially  in  Figure  VIII-2).   This  result 
was  obtained  by  translating  the  ZnS  sample,  located  at  the  focal  point 
of  the  input  20X  objective,  in  the  x  direction  (i.e.,  perpendicular  to 
the  direction  of  light  propagation).  In  this  way  the  effective  angle  of 
incidence  of  the  light  (defined  by  the  central  ray  of  the  cone  of  light 
from  the  objective)  could  be  changed  slightly,  but  not  so  much  as  to 
convert  the  waveguided  light  to  a  higher  order  mode.   Since  the  entire 
puide  (and  not  just  the  top  of  the  guide)  could  be  selectively 


81 


82 


Figure  VIII-1.   End-view  photograph  of  the  cleaved  edge  of  an  18um-thick 
zinc  sulfide  light  guide  illuminated  in  transmission  by 
white  light.   Magnification:   440X.   Some  microcrystalline 
defects  are  evident .   The  gallium  arsenide  substrate  is 
below  the  guide  (i.e.,  below  x  =  -W) ,  the  air  superstrate 
above  the  guide  (i.e.,  above  x  =  0) . 
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Figure  VIII-2.   End-view  photograph  of  the  output  intensity  distribution 

of  the  cleaved-edge  18um-thick  ZnS  guide  in  Figure  VTII-1, 
but  with  a  higher  magnification  (780X) ;  wavelength 
\     =   0.6328pm.   Figure  VIII-2  is  a  mirror  image  of  Figure 
VIII-1,  as  may  be  seen  by  comparing  the  asymmetry  of  the 
crystal  defect  in  the  center  of  Figure  VIII-1  with  the 
(same)  defect  in  the  upper  right  of  Figure  VIII-1.   It 
may  be  noted  that  the  intensity  distribution  has  a 
maximum  near  the  top  of  the  guide.   The  mode  is  m  =  0 
but  approaching  m  =  1. 
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Figure  VIII-3.   Same  as  in  Figure  VIII-2,  only  now  the  intensity  distribu- 
tion has  a  maximum  near  the  center  of  the  guide. 
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Figure  VIII- 4.   Plot  of  the  intensity  distribution  in  Figure  VIII-3 
showing  the  symmetry. 
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Figure  VIII-5.   Same  as  in  Figure  VIII-2,  only  now  the  intensity  distribu- 
tion has  a  maximum  at  the  bottom  of  the  guide. 
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"illuminated"  in  this  manner,  this  suggested  that  the  ZnS  films  were 
fairly  homogeneous  throughout;  if  some  inhomogeneity  is  present, 
then  the  index  profile  is  an  extremely  slowly  varying  one.   Hence, 
continual  internal  refraction  is  probably  not  the  primary  mechanism  for 
light  guiding  within  a  ZnS  film. 

The  question  remains:   what  about  the  effects  of  the  lattice  mis- 
match between  the  ZnS  and  the  GaAs?   During  further  experiments  with 
various  samples  illuminated  in  transmission  with  white  light,  an 
interesting  phenomenon  was  discovered.   Inadvertently,  the  white  light 
bulb  used  for  illumination  overfilled  the  input  20X  objective  with  the 
result  that  the  ZnS-on-GaAs  sample  was  illuminated  in  reflection  also, 
i.e.,  light  reflected  off  the  face  of  the  output  microscope  objective. 
As  a  result,  it  was  observed  that  these  samples  showed  a  thin  black  line 
at  the  ZnS/GaAs  boundary.   This  line  was  not  visible  (or,  at  most,  was 
barely  discernible)  when  the  sample  was  viewed  in  reflected  light  only. 
This  black  line  is  shown  in  Figures  VIII-6  and  VIII-7  and  is  seen  to  be 
about  lym  thick.   The  thickness  of  this  line  did  not  increase  when  a 
voltage  of  two  hundred  volts  or  more  was  applied  to  the  crystal,  as  might 
have  occurred  if  the  line  was  due  to  some  semiconductor  junction  effect. 
Thus,  it  seems  reasonable  to  conclude  that  this  line  is  a  permanent 
"built-inr  part  of  the  ZnS-GaAs  structure. 

This  black  line  must  be  either  zinc  sulfide  or  gallium  arsenide 
since  no  other  material  was  present  when  the  ZnS  film  was  epitaxially 
grown.   (If  there  were  some  other  material,  such  as  oxygen,  then  a  zinc 
oxide  layer  would  have  prevented  growth  of  a  good,  cubic-structured  ZnS 
film.)   If  the  line  were  GaAs,  then  light  would  strongly  reflect  from 


Figure  VIII-6.   End-view  photograph  of  the  cleaved  edge  of  an  llpm-thick 
zinc  sulfide  guide  on  a  gallium  arsenide  substrate 
illuminated  simultaneously  in  transmission  and  in 
reflection  with  white  light.   Magnification:   4A0X.   A 
black  line  of  about  lum  thickness  is  evident  at  the 
ZnS/GaAs  boundary.   The  somewhat  thicker  black  line  at 
the  top  of  the  ZnS  film  is  due  to  a  3um-thick  layer  of 
dielectric  parylene  used  to  isolate  the  metal  electrode 
from  the  zinc  sulfide.   Vertical  microcrystalline 
cracks  in  the  sample  are  also  evident. 
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Figure  VIII-7.   Same  as  Figure  VIII-6  but  with  a  shorter  exposure  to 

enhance  the  contrast  of  the  black  line  at  the  ZnS/GaAs 
boundary. 
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it  (as  seen  from  the  bottom  of  Figures  VIII-6  and  VIII-7).   This  is  due 
to  the  high  index  n  =  3.87,  resulting  in  a  high  normal-incidence  GaAs/air 
reflectivity  of  35%  or  so.   Therefore, the  black  line  is  of  zinc  sulfide 
material.   Further,  since  this  black  line  is  prominent  only  when  the 
white  light  has  travelled  the  length  of  the  thin  film  crystal,  it  seems 
plausible  to  also  conclude  that  this  line  is  a  result  of  the  lattice 
mismatch.   A  lum  layer  corresponds  to  approximately  2000  lattice  spacings . 
As  discussed  in  Chapter  V,  it  is  not  unreasonable  to  expect  that  the 
lattice  mismatch  could  have  effects  extending  beyond  1000  A.   Thus,  it 
appears  that  although  an  inhomogeneity  may  exist  within  an  epitaxial 
zinc  sulfide  film,  it  does  not  extend  beyond  lym  from  the  ZnS/GaAs 
boundary,  a  distance  of  only  5  -  10%  of  the  total  film  thickness.   Thus, 
90  -  95%  of  the  film  is  quite  homogeneous. 

As  stated  in  the  caption  to  Figure  VIII-6,  the  thicker  black  line 
at  the  top  of  the  ZnS  film  is  due  to  a  3pm-thick  superstrate  layer  of 
parylene  (a  transparent  dielectric  material  obtainable  from  Union  Carbide 

Corn..  Bound  Brook,  M.  J.;  refractive  index  n  =  1.63).   Since  n   <  n„, 

a  a    Z 

the  parylene  must  be  a  leaky  light  guide  with  losses  proportional   to 

the  cube  of  the  reciprocal  of  the  guide  thickness  (cf.  Chapter  V) . 

Using  Eq .  V-(3),  the  attenuation  coefficient  a    is  approximately  1.3mm   ; 

m 

since  the  length  of  the  ZnS  crystal  and  narylene  is  5mm,  no  light  can 
propagate  within  the  parylene  but  must,  instead,  eventually  leak  into  the 
ZnS  film.   Further,  since  the  parylene  is  transparent  and  since  only  a 
small  amount  of  light  (^  4%)  will  reflect  from  the  parylene/air  boundary, 
the  illuminating  light   reflecting  from  the  output  microscope  objective 
must  be  transmitted  in  the  -z  direction  into  the  parylene.   This  is  why 
the  parylene  appears  black. 
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The  same  line  of  reasoning  can  be  applied  to  the  zinc  sulfide 
black  line.   In  addition,  the  black  appearance  cannot  be  attributed 
to  scattering,  since  forward  and  backward  scattering  would  tend  to 
illuminate  the  region  as  seen  through  the  microscope.   Hence,  depending 
on  whether  or  not  the  inhomogeneous  region  has  a  smooth  transition 
to  the  more  homogeneous  region,  it  must  be  concluded  that  the  inhomo- 
geneity  is  great  enough  that  either: 

1.  A  very  thin  and  very  leaky  light  guide  has  been  formed  (since  both 
the  substrate  index  n  and  the  superstrate  index  n^     of  the  rest  of  the 
ZnS  film  exceed  the  index  of  the  inhomogeneous  region),  in  which  case 
light  from  the  main  ZnS  guide  is  guided  by  total  internal  reflection 
from  this  region:  or, 

2.  Total  internal  refraction  of  light  from  the  main  ZnS  guide  into  and 
out  of  this  region  is  taking  place. 

For  either  situation,  there  is  thus  a  good  physical  reason,  backed 
up  by  strong  experimental  evidence,  for  believing  that  a  dielectric, 
inhomogeneous, ZnS   isolation  layer  exists  between  the  main  zinc  sulfide 
guide  and  the  gallium  arsenide  substrate.   Thus  the  main  ZnS  guide  is 

probably  not  a  leaky  light  guide,  but  instead  is  a  true  dielectric 

29 
waveguide.   This  conclusion  contradicts  the  statements  by  Martin 

30 
and  by  Hall  and  Yeh   based  on  their  results  and  analyses  with  epitaxial 

ZnS  samples  grown  at  Ttek  Corporation  and  supplied  to  them  for  examination. 
This  is  not  surprising  in  view  of  the  fact  that  the  inhomogeneous  region 
is  on  the  order  of  lym  thick  and  only  the  combination  of  transmission 
and  reflection  illumination  revealed  its  existence. 

Photographs  of  mode  structures  (other  than  m  =  0)  for  light  pro- 
pagating in  the  z'  direction  [the  (Oil) direction]  are  given  in  Figures 

A. 

VIII-8  through  VIII-10,  and,  for  light  propagating  in  the  z  direction 
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[the  (001)  direction],  in  Figures  VIII-11  through  VIII-13.   The  various 

A 

modes  were  obtained  by  translating  the  sample  in  the  x  direction.   It 
may  be  seen  in  Figure  VIII-14  that  the  guide  output  edge  no  longer 
appears  rectangular.   This  is  due  to  the  great  difficulty  in  obtaining 
a  well  polished  edge  down  to  the  dimensions  of  micrometers  (which  is  an 
entirely  different  problem  from  fine  polishing  a  broad  surface  to  such 
dimensions).   As  a  result,  the  light  is  either  scattered  or  refracted 
out  of  the  field  of  view  of  the  microscope  by  the  sliced  edge.   White- 
light  photographs  of  the  sliced  edge  are  not  given  here  since  very  little 
of  the  edge  was  in  focus  for  a  given  microscope  position.   As  mentioned 
in  the  caption  to  Figure  VIII-11,  a  pure  m  =  0  mode  could  not  be  obtained 

A 

for  light  propagation  in  the  z  direction. 

Finally,  it  should  be  mentioned  that  the  mode  patterns  shown  in 
the  figures  are  representative  of  TE,  TM,  or  TE  +  TM  modes.   It  was 
observed  that  when  the  input  polarization  was  changed  (e.g.,  TF  to  TM) , 
it  was  necessary  to  slightly  re-adjust  the  position  of  the  sample  with 
respect  to  the  focal  point  of  the  input  20X  objective  in  order  to  re- 
obtain  the  same  mode  structure. 

Flectro-Optic  Modulation 

After  becoming  familiar  with  the  light  guiding  properties  of 
epitaxial  zinc  sulfide  films,  electro-optic  experiments  were  begun. 
Initially  a  large  amount  (a  few  milliwatts)  of  laser  light  was  coupled 
into  the  guides.   With  the  samples  placed  between  crossed  polarizers  and 
a  voltage  applied,  a  transverse  electro-optic  effect  was  observed. 
However,  for  those  samples  with  parylene  between  the  electrode  and  the 
ZnS  film,  the  light  transmitted  upon  initial  application  of  the  voltage 
gradually  decreased  in  intensity  with  time  (over  a  few  seconds).   Upon 
re-application  of  the  voltage  (i.e.,  off  then  on  again),  the  transmitted 
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Figure  VIII-8.   End-view  photograph  of  the  cleaved  edge  of  an  18ym-thick 

zinc  sulfide  film  showing  an  m  =  1  mode  pattern  (i.e., 

one  horizontal  node);  wavelength  A  =  0.6328pm. 
Magnification:   440X. 
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Figure  VIII-9.   Same  as  Figure  VIII-8  but  for  an  m  =  2  mode  pattern 

(i.e.,  tvjo  horizontal  nodes).   On  the  left  an  m  =  3  mode 
pattern  can  be  seen. 
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Figure  VIII-10.   Same  as  Figure  VIII-8  but  for  an  m  =  5  mode  pattern. 
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Figure  VIII-11.   End-view  photograph  of  the  sliced  edge  of  a  20ym-thick 
zinc  sulfide  film  showing  an  m  =  1  mode  pattern;  wave- 
length X     =  0.6328pm.   Magnification:   440X.   (A  pure 
m  =  0  mode  pattern  was  not  obtainable.) 
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Figure  VIII-12.   Same  as  Figure  VIII-11  but  for  an  m=2  mode  pattern. 
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Figure  VIII-13.   Same  as  Figure  VIII-11  but  showing  how  the  edge  of  a 
sliced-edge  ZnS  sample  does  not  appear  rectangular 
because  of  the  difficulty  in  polishing  the  edge,  a 
problem  which  is  not  present  with  the  cleaved-edge 
samples. 
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light  reached  its  previous  peak,  but  again  gradually  decreased.   This 

27  31 
is  the  result  of  a  small  photoconductive  sensitivity   '    in  the  red 

region  of  the  visible  spectrum  due  to  absorption  from  small  amounts  of 

impurities  in  the  zinc  sulfide  film.   [Since  parylene  (deposited  on  top 

of  the  ZnS)  has  an  extremely  high  resistivity,  the  photoproduced  carriers 

built  up  at  the  ZnS/parylene  interface  reduced  the  voltage  in  the  ZnS  film. ] 

Subsequently,  the  intensity  of  the  input  light  was  reduced  to  a  few 

microwatts,  after  which  no  further  significant  photoconductive  effect 

was  present  to  interfere  with  the  electro-optic  measurements.   (Although 

the  voltages  applied  were  dc  only,  at  modulation  frequencies  above  1 

Hz  or  so  the  photoconductive  effect  should  not  be  significant  even  for 

milliwatt  or  greater  intensity  levels.   This  is  because  the  charge  carriers 

generated  by  the  light  can  no  longer  follow  the  changes  in  the  applied 

31 
electric  field.) 

For  quantitative  measurements,  a  highly  light-sensitive  photomultiplier 
tube  was  used  as  a  detector  since  the  photoconductive  relaxation  effect 
took  place  for  high  intensities, precluding  the  use  of  less  sensitive 
detectors.   For  qualitative  visual  observations,  a  greater  light  intensity 
was  required,  but  since  a  switch  was  used  manually  to  apply  and  re-apply 
the  voltages,  this  did  not  present  any  serious  difficulties. 

Three  samples  were  investigated,  as  indicated  in  Table  VIII-1. 
With  one  sample  (Sample  #1),  modes  m  =  0,2,  and  4  were  obtainable  in 
pure  form.   (Apparently,  proper  phase  conditions  for  propagating  a 
pure  m  =  1  or  3  mode  were  not  met.)   This  sample  also  resulted  in  the 
greatest  amount  (95%)  of  light  extinction.   For  the  other  two  samples, 
only  a  pure  m  =  0  mode  was  obtainable:  all  other  mode  patterns  were 
highly  multi-moded.   Typical  results  of  the  observed  electro-optic 
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modulation  are  presented  in  Figures  VIII-14  through  VIII-18  and  are 
summarized  in  Table  VIII-1.   Since  Sample  #2  had  a  3ym-thick  layer  of 
parylene  between  the  electrode  and  the  ZnS,  Eq .  VI- (42)  was  used  to 
calculate  the  voltage  across  just  the  ZnS  crystal  (using  a  value  of 
E,  =  3.1  for  the  parylene). 

From  the  graphs  it  may  be  seen  that   with  the  ZnS  sample  placed 
between  crossed  polarizers  the  minimum  is  shifted  away  from  the  V  =  0 
situation.   This  can  be  attributed  to  some  residual  birefringence  within 
the  ZnS  crystal,  as  well  as  to  the  fact  that,  as  stated  in  Chapter  II 
(d.  12),  a  TE  +  TM  wave  initially  linearly  polarized  will  become 
elliptically  polarized  as  it  propagates  down  the  guide. 

For  the  modulation  curves  of  the  m  -  2  and  4  mode  patterns  (Figures 
VIII-15  and  VIII-17) ,  it  is  evident  that  some  extinction  was  obtainable. 
It  should  be  noted,  however,  that  only  extinction  occurred,  i.e.,  no 
increased  transmission  occurred  for  the  m  /  0  modes.   Visually,  it 
was  observed  that  the  mode  pattern  would  degenerate  into  a  highly 
multi-moded  pattern  as  the  voltage  was  increased;  at  the  same  time,  the 
over-all  output  intensity  first  decreased  then  increased  back  again  to 
the  intensity  present  for  V  =  0.   Similar  measurements  for  light  pro- 
pagating in  the  (001)  direction  gave  no  modulation  of  the  light  for  any 
polarization  or  any  mode,  as  predicted  by  theory.   However,  in  the  (011) 
case,  visual  observations  revealed  that  for  a  given  mode  pattern  at 
V  =  0,  the  pattern  was  re-established  for  a  voltage  approximately  halfway 
between  the  minimum  and  the  second  maximum.   This  can  probably  be 
explained  by  the  fact  that  the  proper  phase  conditions  for  propagating  a 
given  mode  were  not  met  when  the  voltage  corresponding  to  the  maximum 
output  signal  (which  was  multi-mode)  was  applied.   Another  way  of  stating 
the  above  is  the  following.   Although  some  amplitude  modulation  was 
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Figure  VIII-14.   Plot  of  the  transmitted  intensity  I/I  versus  applied 
voltage  from  V  =  0  to  250  volts  for  an  m  =  0  mode  in 
Sample  #1. 
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Figure  VIII-15.   Plot  of  I/I   from  V  =  0  to  500  volts  for  an  m  =  2  mode 
in  Sample  #1. 
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Figure  VIII-16.   Same  as  Figure  VIII-15  but  for  an  m  =  k   mode. 
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Figure  VIII-17.   Plot  of  i/l  from  V  =  0  to  200  volts  for  an  m  =  0  mode 
in  Sample  #2. 
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Figure  VIII-18.   Plot  of  I/I   from  V  =  0  to  200  volts  for  an  m  =  0  mode 
in  Sample  #3. 


Table  VIII-1.   Summary  of  Characteristics  of  Samples  Used  for  Electro- 
Optic  Experiments 
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Sample  # 

Thickness 

Interaction 

Width 

Experimentally- 

Predicted 

W 

Length  L 

h 

Determined 
Value  of 

(VA/2)Trans 

* 

Value   of 

^VX/2^ Trans 

1 

15ym 
±  1.5pm 

4  .0mm 
±  0.25mm 

3mm 

94  volts 
±  5  volts 

103  volts 
+  5  volts 

2 

25um 
±  1.5ym 

8.0mm 
±  0.25mm 

2  mm 

80  volts 
+  5  volts 

86  volts 
±  5  volts 

3 

17.5pm 
±  1.5pm 

4 . 25mm 
+  0.25mm 

4 .  5mm 

125  volts 
±  5  volts 

113  volts 
±  5  volts 

Obtained  by  substituting  W  and  L  into  Eq .  VI-(ll),  where  r 
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,-10 


1.82  x  10    cm/volt  and  n  =  2.33  (Ref.  18),  and  X   =  0.6328ym. 


The  measured  voltage  V  was  105  volts  ±  5  volts,  but  since  Sample 
#2  had  a  3um-thick  layer  of  parylene,  the  value  given  above  was 
determined  from  Eq .  VI-(42)  for  the  voltage  V  across  the  zinc  sulfide 
only. 
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observed  for  modes  m  >  0,  the  primary  effect  was  phase  modulation,  resulting 
in  mode  conversion.   Finally,  continued  experimentation  with  the  high 

voltages  indicated  on  the  graphs  resulted  in  destructive  dielectric  break- 
down of  some  of  the  samples. 

In  Table  VIII-1  it  may  be  noticed  that  the  measured  values  of  the 
half-wave  voltage  correspond  closely  with  the  values  predicted  for  the 
bulk-crystal  case  [from  Eq .  VI- (11)],  which  was  to  be  expected  as  discussed 
in  Chapter  VI.   It  was  found  that,  even  though  the  measured  values  of 
^VA/2^Trans'  pluS  or  minus  the  experimental  error,  overlapped  the 
predicted  values,  nevertheless  the  predicted  values  were  consistently 
lower  than  the  measured  values.   Thus,  in  calculating  the  expected 

values,  the  bulk-crystal  value   of  the  electro-optic  coefficient  r,„  was 

41 

1  R 
not  used  in  Eq .  VI- (11).   Instead,  a  measured  value   of  r.,  for  Ttek- 

41 

grown  epitaxial  zinc  sulfide  films  was  used.   This  value  was  obtained 
from  longitudinal  electro-optic  experiments  for  light  incident  upon  the 
ZnS  from  the   x  direction.   The  result  is:   r,_  =  1.82  x  10    cm/volt, 
which  is  12%  lower  than  the  bulk-crystal  value.   It  may  be  seen  that  the 
results  of  the  transverse  electro-optic  waveguiding  experiments  presented 

I  Q 

here  tend  to  corroborate  the  measurements  of  r,   obtained  previously 
in  separate  experiments  in  which  the  light  was  propagated  along  an 
entirely  different  direction  in  the  crystal. 


SUMMARY 

This  dissertation  presents  the  results  of  a  theoretical  and  an 
experimental  research  effort  on  general  index  gradient  waveguides  and 
their  relationship  to  guiding  and  electro-optic  modulation  of  light  in 
thin-film  zinc  sulfide  crystals.   The  basic  physics  of  homogeneous 
asymmetric  waveguides  has  been  reviewed,  allowing  ready  comparison  of 
similar  concepts  (such  as  the  characteristic  waveguide  equation)  between 
homogeneous  and  inhomogeneous  optical  waveguides.   It  has  been  shown  that 
a  quasi-geometrical  ray-optics  approach  leads  to  the  same  results 
obtained  from  more  rigorous  mathematics,  permitting  a  practical  understand- 
ing of  the  phenomenon  of  total  internal  refraction  and/or  reflection  in 
general-index-gradient  nlanar  optical  waveguides.   As  a  result,  the 
characteristic  waveguide  equation  for  describing  the  propagation  of 
light  in  such  guides  can  be  easily  derived.   This  equation  has  general 
applicability  to  any  slowly  varying  monotonic  index  nrofile  and  allows 
computation  of  various  waveguide  propagation  parameters.   Of  particular 
note,  the  phase  change  of  tt/2  (not  0  or  it)  occurring  upon  total  internal 
refraction  has  been  explained  as  being  due  to  a  cylindrical  focussing  of 
the  light  within  the  guide. 

Waveguiding  and  the  electro-optic  modulation  of  light  in  epitaxial 
zinc  sulfide  films  on  gallium  arsenide  substrates  have  also  been  discussed. 
Two  mechanisms  for  light  guiding  at  the  ZnS/GaAs  boundary  have  been 
described:  total  internal  refraction  (resulting  from  a  4%  lattice 
mismatch)  or  partial  internal  reflection  (since  the  refractive  index  of 
GaAs  exceeds  that  of  ZnS) .   In  the  course  of  waveguiding  experiments, 
the  discovery  was  made  of  a  one-micrometer-thick  region  of  inhomogeneous 
zinc  sulfide  which  is  transparent  in  visible  light  and  which  has  a  lower 
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index  than  the  rest  of  the  ZnS.   As  a  result,  it  is  concluded  that 
1)  90-95%  of  the  zinc  sulfide  thin  film  is  homogeneous,  and  2)  due  to  the 
thin  inhomogeneous  layer,  the  zinc  sulfide  guide  is  isolated  from  the 
gallium  arsenide  substrate;  therefore  the  guide  is  not  a  leaky  one  but 
instead  is  a  true  dielectric  waveguide.   Finally,  electro-optic 
measurements  have  been  presented  which  agree  with  those  predicted  from 
theory.   In  addition,  results  from  these  experiments  have  been  shown  to 
corroborate  previously  obtained  measurements  by  other  authors  on  the 
electro-optic  coefficient  of  Itek-grown  epitaxial  zinc  sulfide  films. 


APPENDIX  A 

APPLICATION  OF  MAXWELL'S  EQUATIONS  AND  VECTOR  POTENTIALS 

TO  HOMOGENEOUS  AND  INHOMOGENEOUS  PLANAR  OPTICAL  WAVEGUIDES 

Maxwell's  Equations 
Maxwell's  equations  are: 

V-$=<iirp  a-(d 


V-B  ^0 


A-(2) 


V         C  IT  A-(3) 


X  *  ft   __  T  c  ^ 


A- (4) 


where  D  =  eE,  B  =  uH,  and  c  is  the  soeed  of  light  in  free  space.   For 

2 
a  non-permeable  medium  (u=l)  ,  £=n  .   For  an  isotropric  dielectric  with 

no  charge  density  (p=0)  and  no  currents  (J=0) ,  Eqs.  (l)-(4)  can  be  re- 
written as: 


A-(5) 


A- (6) 


v,r=-xM 
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Vxftr  X  ±-CrTe)  A_(8) 


For  a  homogeneous  medium,  Eq .  (5)  may  be  re-written  as: 

nfl  V-f*o  A_(9) 

Assuming  that  n  is  also  not  a  function  of  time,  Eq .  (8)  may  be 
expressed  as: 


C    ^t  A- (10) 


Taking  the  curl  of  Eq .  (7)  and  interchanging  the  curl  and  time  deriva- 
tive operations  gives: 

VX  C   ^tV  A-(ll) 

Using  a  trigonometric  identity  for  the  left-hand  side  (L.  H.  S.)  of 
Eq.  (11),  and  substituting  Eq .  (10)  into  the  right-hand  side  (R.  H.  S.) 
yields 

.2  .a=* 


*(7-*)-(v.r)F--£0 


By  Eq.  (9),  the  first  term  of  the  L.  H.  S.  of  Eq .  (12)  is  zero.  Writing 

2 
V  '  V  as  V  ,  Eq.  (12)  becomes  the  wave  equation: 
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Similarly,  a  wave  equation  for  H  may  be  obtained  by  taking  the  curl 
of  Eq.  (10): 


Ca  ±t  ^  >  A-(14) 


or,  by  Fq.  (7) , 


Thus 


A  plane-wave   solution   to   Eq .    (13)    is   of    the    form 
00 


A-(15) 


A-(16) 


A-(17) 


A- (18) 


A- (19) 


.A     A       A       AAA  A 

Here  w  is  the  angular  frequency  and  r  =  xx  +  yy  +  zz  (x,  y,  and  z 

are  unit  normals).   Also,  k  is  the  propagation  vector  which  has 

a  direction  normal  to  the  wavefront  and  a  magnitude  of  nk  ,  where 

o 

n   is    the   index  of   refraction   and   k     =   co/c   =   2tt/A      (and    \      is    the 

o  oo 

free-space  wavelength) ;  6  is  the  wavevector  component  along  the 

A 

z-axis  [as  shown  in  Figure  II-l]  and  has  the  value  k  n   sin6. 

°   g 


It  may  be  noted  that  any  x  or  y  dependence  has  been 

.a 
absorbed  into  the  vector  E  (x,y) ,  where 


eo(m)=  *X +  $*•« -^5 


A-(20) 


A  similar  solution  exists  for  H: 


A- (21) 


where 


K^3>^*.^3  V*H^ 


A-(22) 


From  the  geometry  of  Figure  II-l  for  a  planar  (slab)  waveguide, 

«*     -* 

E  and  H  have  no  y-dependence.   Hence, 

oo         J 


*.=  *.M  ■  e  =  e(*,z) 


A-(23) 


and 


lf,«  K(0  ;    H=  i+M 


A- (24) 


and  all  8/8y  terms  in  Eqs.  (13)  and  (16)  are  zero. 


Vector  Potentials 


icot 


Returning  now  to  Eqs.  (5)  -  (8),  if  the  time  dependence  e 
in  Eqs.  (19)  and  (21)  is  used,  then  the  following  set  of  equations 


result  (assuming  n  is  not  a  function  of  time) : 

V-  iff)  -  o 


2-»       -» 
In  view  of  the  divergenceless  character  of  n  E  and  of  H,  it 

is  possible  to  express  the  fields  in  terms  of  a  magnetic 

vector  potential  A  or  in  terms  of  an  electric  vector  potential 

-»  1  2  ■* 

F.   '   Any  vector  field  H  whose  divergence  vanishes,  as  in  Eq. 

(26) ,  can  be  written  as 


Substituting  Eq .  (29)  into  Eq.  (27)  produces 

Any  curl-free  vector  field  can  be  expressed  as  the  gradient 
of  a  scalar  potential.   Thus 

E   +  CA0A  =  V  f 


A- (25) 


A-(26) 


A-(27) 


A-(28) 


A-(29) 


A- (30) 


A  -(31) 
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A  -» 

where  i>     is  an  electric  scalar  potential.   Re-writing  Eq.  (31),  E 

may  be  written  as 

^  ~  -  iM  i  -+-  v  ^A 

<-^on  T     V   T  A- (32) 


2-* 
In  an  exactly  analogous  fashion,  Eq.  (25)  allows  n  E  to  be 

written  as 


Substituting  Eq.  (33)  into  Eq.  (28)  produces 
which  will  be  satisfied  if 


A-(33) 


A- (34) 


A- (35) 


■p 
Here,  ty     is  a  magnetic  scalar  potential.   Eqs.  (32)  and  (35)  are 

dual  equations. 

To  obtain  an  equation  for  A,  Eqs.  (29)  and  (32)  are  substituted 

into  Eq.  (28) : 

b  A- (36) 


To  obtain  an  equation  for  F,  Eq.  (33)  and  (35)  are  substituted  into 
Eq.  (27): 
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**[^V*^f)]=  -A.Vf  -c\V^ 


A-(37) 


Returning  now  to  Eqs.  (27),  (28),  (29),  and  (33),  the  following 
may  be  said.   Solutions  for  E  are  given  by  Eq.  (28)  and  also  by  Eq . 
(33) ;  hence  the  superposition  of  these  two  partial  solutions  is  also 
a  solution,   i.e., 

which,  by  Eq .  (29),  becomes 

f»  ~^    Wf)  -  WKrf)VxVa  A-(39, 

Similarly,  a  complete  solution  for  H  is  given  by  the  superposition 
of  Eqs.  (27)  and  (29): 

H~  YKA  -r-  (t/K)  Vxi  A- (40) 

which,  by  Eq.  (33),  becomes 


A-(41) 


Homogeneous  Waveguides 
For  homogeneous  waveguides,  the  refractive  index  is  not  a 
function  of  position;  hence  Eq.  (37)  reduces  to 


A- (42) 


11/ 


Using  the  vector  identity 

Eqs.  (36)  and  (42)  can  be  written  as 


A- (43) 


A- (44) 


A-(45) 


There  is  a  good  deal  of  arbitrariness  in  the  choice  of  the  vector 
potentials  A  and  F,  and  it  is  possible  to  specify  the  divergence  of  them. 


Thus  letting  i 

-  l.   3  ,|,A 


A-(46) 


A- (47) 


then  Eqs.  (44)  and  (45)  reduce  to  the  wave  equations 


VaA-t->„VA  =  o  A_(48, 


?=F  +-4.VP-0  A-(49) 

For  homogeneous  waveguides,  Eqs.  (39)  and  (41)  reduce  to 

E=  ~  VxF  -(c/*.tta)v*Vxfi  A_(5o) 
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Ff=  Vx7~(c/h)VxV*F 


A-(51) 


which,  by  Eq .  (43),  become 


A- (52) 


A- (53) 


At  this  point  it  is  convenient  to  consider  some  particular  choice 
of  potentials.   Letting 


A-0 


A- (54) 


P^*P 


A- (55) 


then  the  wave  equations  [Eqs.  (48)  and  (49)]  reduce  to 


4r 


^^;a-o 


A-(56) 


where  the  3/9y  term  has  been  set  equal  to  zero  for  the  slab  waveguide 
case.   Further  the  equations  for  E  and  H  [Eqs.  (52)  and  (53)]  become: 


£  =  -Vx  f=  - 


^x       u     ±z 

Px      0      0 


A-(57) 


f=  ~t/K)iv(m-  7af] 


A-(58) 
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Expanding  Eqs.  (57)  and  (58)  into  rectangular  coordinates  results  in 
the  following: 


£   =0 

*  A-(59) 


13      ^*  A- (60) 


E  =    O 

Z    w  A-(61) 


*„=^M^ 


A- (62) 


^~  0  A-(63) 


^=-C^) 


^x 


^X  A"(64> 


For  this  choice  of  A  and  F  it  may  be  seen  that  there  is  no  z-component 
of  E;  for  this  reason  this  is  called  the  transverse  electric  (TE)  mode 
of  propagation.  Letting 

-*    A 

A  ~  X  Ax  A-(65) 

r*  o 

A-(66) 

then  the  wave  equations  reduce  to 


*1*      d-x"  *  A-C67) 


From  Eqs.  (52)  and  (53),  the  following  equations  for  E  and  H  result 


--f^)^fe-)-v^j 


A- (68) 


tt=     Vk(\* 


A     a.     A 
K         14    Z 


3 


*>*      0     0 


A-(69) 


or,  in  rectangular  coordinates, 


A-(70) 


V° 


e =  4M 


A-(71) 


A- (72) 


*x  =  0 


A-(73) 


*..= 


_    iA, 


^       *z 


A- (74) 


na=o 


A-(75) 


For  this  choice  of  A  and  F  there  is  no  z-component  of  H;  hence  this 
is  called  the  transverse  magnetic  (TM)  mode  of  propagation. 


Finally,  for  a  situation  where  neither  a  pure  TE  nor  a  pure  TM 
exists,  the  following  choice  of  A  and  F  is  useful: 


A-xA, 


Then,    from  Eqs.     (52)    and    (53) 


X  W      oy'S    ^2.a 


£  =  KM 


v-(^)- 


^F. 


±Z: 


14 - 


-    ^ 


A- (76) 


F=  x  F 

*  A-(77) 


A-(78) 


A-(79) 


*  4*dX  A-(80) 


A- (81) 


^  ^  A- (82) 


y  ^zax  A"(83) 


It  may  be  seen  that  Eqs.  (76)  -  (83)  are  nothing  more  than 
the  superposition  of  Eqs.  (60)  -  (65)  [the  TE  case]  and  Eqs. 
(70)  -  (75)  [the  TM  case].   Thus  a  field  of  any  polarization 
in  a  homogeneous  planar  guide  can  be  expressed  as  the  sum  of 
a  TE  field  and  a  TM  field. 


Inhomogeneous  Waveguides 
Returning  now  to  Eq .  (37),  use  of  the  vector  identity 


A- (84) 


allows  Eq.  (37)  to  be  written  as 


v*[£  (?«*k?  +  n* Vx f)] =-^0Vf- ;a„ Vt F 


A~(85) 


The  assumption  is  now  made  that  if  a  refractive  index  gradient  exists, 
then  it  lies  in  the  x  direction  only,  i.e., 


Y[-V\(x) 


A- (86) 


Further,  F  and  A  are  still  chosen  to  lie  in  the  x  direction  [i.e., 

F  =  xF  ,  A  =  xF  ]  as  in  the  homogeneous  case.   As  a  consequence,  the 
xx 

term  Vn   x  F  in  Eq.  (85)  has  the  value 


A 

A      A 

X 

3      * 

dj£ 

0  C) 

*>X 

Fx 

0  0 

-  0 


A- (87) 


Hence,  Eq.  (85)  reduces  to 


A-(88) 


which  is  identical  to  Eq.  (42)  for  the  homogeneous  case.   Eq.  (36) 

a 
for  A  is  still  valid  for  the  inhomogeneous  case. 

Using  Eq.  (43),  Eqs.  (36)  and  (88)  can  he  written  in  the  form  of 

Eqs.  (44)  and  (45).   Then  Eqs.  (46)  and  (47)  are  again  used  as  the  choice 

for  the  divergence  of  the  vector  potentials  A  and  F.   Substitution  of 

Eqs.  (46)  and  (47)  into  Eqs.  (44)  and  (45)  yield  the  following  equations. 


A-(89) 


A- (90) 


Since  'i:A  =  (-i/k  n2)V  *  A.  Eq .  (89)  becomes 


n  -0/^)(V'*)CVr?)±K^O 


A-(91) 


It  may  be  noted  that  Eq .  (91)  for  A  for  the  inhomogeneous  case  is 

identical  to  Eq .  (48)  for  the  homogeneous  case,  except  for  the  term 

•*    2  -* 

-  (V'A) (Vn  )  which  is  zero  for  n  =  constant.   Eq.  (90)  for  F  is 

•a 
identical  to  Eq.  (49),  i.e.,  F  obeys  the  same  wave  equation  for  both  the 

homogeneous  and  inhomogeneous  cases.   Using  Eqs.  (84)  and  (85),  then 

Eqs.  (90)  and  (91)  become  the  wave  equations 


0  +£5  *4.V&0F,-o 

d2       ^Ka  A-(92) 
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The  equations  for  E  and  H  are  given  by  Eqs.  (39)  and  (41). 
Employing  the  vector  identies  given  in  Eqs.  (43)  and  (84),  Eqs.  (39) 
and  (41)  can  be  written  as 


E=  ~^  (*»'+  f  W  Fx  F)-(i/>.ii')[7f7.j?)-7^] 


fB  VKA~(^)7x[^(^%p+nVxf)] 


A- (94) 


A-(95) 


2   -* 
Since  Vn   x  F  is  zero,  Eqs.  (86)  and  (87)  become 


£  -   -VxF  -(^0na)[^(r.A)-^A] 


ff*     Vx/t    ~(C/K)VxVxF 


-   7*A  -(c/Ao)[7^f)-^p| 


A-(96) 


A-(97) 


A-(98) 


Again,  using  the  choices  A  =  xA  and  F  =  xF  ,  E  and  H  in  Eqs.  (96) 
and  (98)  then  become,  in  rectangular  coordinates, 


**-  LtfAMOl  % 


A- (99) 


e  =  -*S 


** 


A- (100) 


LV  - 


A-(101) 


ik- 


A- (102) 


4.= 


-    ^A* 


&        ** 


A-(103) 


n««  ~(;/k) 


A-(104) 


It  may  be  noted  that  Eqs.  (99)  and  (104)  are  identical  to  Eqs.  (86)  - 

(91)  if  the  refractive  index  n  is  a  constant.   Further,  it  was  the  choice 

A  =  xA   and  F  =  xF   (i.e.,  that  A  and  F  lie  in  the  direction  of  the 
x  x 

index  gradient)  that  resulted  in  considerable  simplification  of  the 
inhomogeneous  waveguide  equations,  especially  Eqs.  (85),  (94)  and  (95). 
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~*  "*"   A 

In  the  homogeneous  case,  the  choice  A  =  0,  F  =  xF   led  to  Eqs. 

x         n 

(60)  -  (65),  the  TE  mode  of  propagation.   So  too  for  the  inhomogeneous 

•A 

case.   If  A  =  0,  then  E  ,  E  ,  and  H   in  Eqs.  (93)  -  (98)  are  identically 

«k 
zero,  producing  a  pure  TE  mode.   If  F  is  set  equal  to  zero,  then  E  , 

y 

Hx,  and  Hz  in  Eqs.  (101)  -  (104)  are  identically  zero,  producing  a 
pure  TM  mode.   Further,  as  with  the  homogeneous  case,  for  the  situation 
where  neither  a  pure  TE  nor  a  pure  TM  exists,  by  choosing  A  and  F  to 
be  both  non-zero,  a  field  of  any  polarization  in  an  inhomogeneous 
planar  guide  can  be  expressed  as  the  sum  of  a  TE  field  and  a  TM  field. 


APPENDIX  B 
HOMOGENEOUS  ASYMMETRIC  WAVEGUIDES 


TE  Modes 

In  order  to  obtain  solutions  for  E  and  H  for  the  TF.  case,  the 

-*   A 

following  choice  is  made  for  E  [F  =  xE  ] : 

x       x 


B-(l) 


^C#>e(ut^%^5X+F3e->] 


3-     j        B-(2) 

■for  -Whx^O 


^sK^^)^^^^         W    X*-*       B"(3) 


where  F   F0 ,  F~,  and  F,  are  constants  (possibly  complex).   The  choice 

for  F  given  in  Eqs.  (1)  -  (3),  when  substituted  into  Eqs.  A-(59) 
through  A- (64),  insures  that  the  wave  is  confined  to  the  n   region 
(the  dielectric  waveguide)  and  that  the  field  goes  to  zero  as  x— ►  +  °°. 
It  is  also  consistent  with  the  form  of  E  and  ~F  given  by  Eqs.  A- (19) 
and  A-(21). 
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Substituting  Eqs.  (1)  -  (3)  into  Eq.  A- (56),  the  wave  equation 
for  F,  yields 

O  B-(4) 


(*a-^%*:V)Fxl= 


:^K'^i9a  -fo\  #$  fe  •*"  ♦ ^  -  o     B"<5) 


sa-p/+^ns3)F„=0 


B-(6) 


isx      — i£x 
Since  e    and  e     are  independent  functions,  Eq.  (5)  can  be  re- 
written as  two  equations: 


i<txV_„3  .aJ  ,  ia^\_ 


Re-expressing  Eqs.  (4)  -  (8)  leads  to  the  following: 


a  -    (I  a  ..  a    .a 


B-(7) 


B-(8) 


(C  =  A.V  + 


B-(9) 
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a   2^  n  2  „,  2    ; 


B-(10) 


(V*~  A>a"5a  +  5' 


B-(ll) 


It  is  possible  to  determine  the  relationships  among  B  .  6  . 

and  3  by  imposing  the  standard  boundary  conditions  at  the  dielectric 

_*     fc 

interfaces,  namely  that  the  components  of  E  and  H  (E  and  H  )  which 

y     z 

are  tangential  to  the  n   -  n  and  n   -  n  boundaries  be  continuous 
a    g      g    s 

across  the  boundaries. 


isv  m 


B-CL2) 


M.  -  M 


5 


ot    X-0 


B-C13) 


fc,].=     DM 


w,  -  w 


at  X*-W 


B-(14) 


B-(15) 


However,  Eqs.  A- (60)  and  A- (64)  imply  that  the  boundary  conditions 

on  Hz  are  equivalent  to  the  requirement  that  (i/k  )9E  /3x  be  continuous 

across  the  boundaries.   Thus  Eqs.  (13)  and  (15)  can  be  re-written  as 
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$e< 


L  ^x 


at  y-0 


B-(16) 


A3 


flutjC'-Vl/  B-(17) 


From  Eqs.  (1),  (2),  and  A- (60),  Eq.  (12)  implies  that,  at  the  n  -  n 


boundary, 


1  a.  xx 


X^O 


-  *h  % 


vx=Q 


B-(18) 


or, 


B-(19) 


hence. 


^^^(FtFje^^ 


B-(20) 


Since  F  ,  F„,  and  F   are  constants,  for  arbitrary  z  the  only  possible 


relationship  between  (3   and  3   is 
a      g 
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?-'?<. 


B-(21) 


Similarly,  at  the  n   -  n  boundary,  Eqs.  (2),  (3),  (14)  and  A-(60) 

6       B 


imply  that 


^•F«Lw-ftF' 


X=~W      '  O   3  'x  =  -W  B-(22) 


or, 


F  e-^jZ -5M/     .{A-zTp.    -<3W         .  w7 


23) 


-sW     -ig''T         ieW 
Since  F^e    ,  F?e   n  ,  and  F  e  M     are  constants,  the  following  relationship 


must  exist 


^ 


B-(24) 


Recalling  that  £  is  the  wavevector  component  along  the  z-axis  and  that, 

for  the  n    region  (the  waveguide),  g  =  k  n  sin6,  then  it  is  convenient 
8  o  g 

to  let 
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(3=  jkQn^^&    y^^=&s  B-(25) 


where  0  is  defined  in  Figure  II-l.   Eq .  (25)  thus  implies  that  the 

wave  propagates  in  the  z-direction  at  the  same  velocity  v  in  the  n  , 

n  ,  and  n  regions,  since  v  =  cB. 
g'      s  z 

Eqs.  (9)  -  (11)  can  now  be  transformed  with  the  aid  of  Eq.  (25): 


^  &An$sCr?9--n?) 


a=  k(.»{"-  ^5^5&y 


B-(26) 


B-(27) 


°   3  B-(28) 


s>  =  A  (^  s^^G--  y\s a) 


a\7A 


B-(29) 


Use  of  Eqs.  A- (60),  and  A-(64)  allows  Eqs.  (12),  (16),  (14),  and 
(17)  to  be  written  respectively  as 


lF*4^0^     L^J^q  8-<30) 


x-O 
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B-(31) 


i^l^~~    I 


F 


3  K---*J 


B-(32) 


ur. 


<L 


=.  uf 


^x 


x  — W 


x=-W 


B-(33) 


From  Eqs.    (1)    -    (3)    and    (25),    Eqs.    (30)    -    (33)   become 


B-(34) 


-aF>=  faF»-caFa 


B-(35) 


F,<f*".   F.«-<«W+F  e*W 


B-(36) 


*V    S   ^F^   8   ^3^      B-C37) 


iaW 


Multiplying  Eq.  (34)  by  +a   and  adding  it  to  Eq.  (35)  produces 


£(«-+^  =  -F3(*-fc3) 


B-(38) 
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Multiplying  Eq .  (36)  by  -s  and  adding  it  to  Eq.  (37)  produces 


F„  (-«♦'• 


<:qW 


:3)e-'9     =F3Cs+^)e,'3 


,\»W 


B-(39) 


Dividing  Eq.  (38)  by  Eq.  (39)  results  in 


^  e^    ^      ,aW 


S-c< 


?+cv 


& 


B-(40) 


a^lN_     (o^)   (s^) 


B-(41) 


From  complex  algebra 


Z=  X,  l  I=ti(,I#u^e1^ 


B-(42) 


z?=?;-iy=-i[+lx+1)=-Ux*+x2  e 


=r^  4lW(2/x) 


B-(43) 


and 


Z*      X-X         ^W(x/t) 
-  -e 


Z        X+lT 


B-(44) 
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Thus,  Eq.  (41)  can  be  re-written  as 


<~      ~  €  e  B-(45) 


2igW     ,     .        2igW  -(2innr)  . 
With  no  loss  of  generality,  e     can  be  written  as  e    e        where 

m  =   0,    1,    2,    ...     .      Thus 


e*p[c  (agW-a^Tr)J~e^|t[^tityrV(V^)+At 


'-"'■B-(46) 


aaW-a^-*^ -amir  B_(47) 


where 


1*"3  d  B-C48) 


?7  =  w  h/%) 


B-(49) 


Thus,    from  Eqs.     (26),    (28),    and    (29), 
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€[--  ^"'KV^'e-  -<tls  o-o-l 


B-(50) 


(f£»    W'[(KijaA^ft9--iisa/a/naCod^]  b-(51) 


Using  Eq .  (28),  Eq .  (47)  can  be  written  as 


a  A.  w  *  «we-  -  9  <f^"  -  a  <>37  =  amir    ,.02) 


It  may  be  noted  that  Eq .  (52),  the  waveguide  equation,  is  a  transcen- 
dental equation  in  0. 

Eq.  (52)  is  useful  for  calculating  the  minimum  thickness  (W   ) 

min'TE 

necessary  for  propagating  a  given  mode  of  order  m.   Solving  for  W  gives: 


W«  L<£  +  <$£  +mr)/A*j ^^  B-(53) 


When  sin  6  =  n  /n   =  sin  G   ,  then  6  is  at  the  critical  angle  8    for 
s      g        gs'  gs 

total  internal  reflection  at  the  n   -  n  boundary.   Since  n   >  n  ,  total 

g    s  s    a 

reflection  still  occurs  at  the  n   -  n  boundary.   However,  if  6  <  6   , 

a    g  gs' 

plane  waves  escape  or  leak  out  into  the  n  medium  and  propagation  ceases, 

i.e.,  the  mode  is  cut  off.   The  cutoff  angle  8    is  therefore 

'  °    cc 
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&>V&CC   5    S£\&as=   "J"} 


B-(54) 


By  Eqs.  (50)  and  (54),  at  cutoff, 

^  ^[W-^ti-^t]   ,(55) 


^  V'a   /       =/  ~\'lx 


**>©-*  ((-^\)  «(l-ns7<) 


3s      v'   ^  ^s;     ^    '*'  V         b-(56) 


Thus,  for  TE  modes  at  cutoff, 


Further,  by   Eqs.  (51)  and  (54),  for  TE  modes  at  cutoff, 


T3S 


B-(57) 


B-(58) 
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2     2  1/2 
Since  1/k  n  cos  9   =  X  /2ir(n   -  n   )    ,  Eq.  (53)  may  be  written  for 
o  g      gs    o     g     s 

TE  modes  as 


In  order  to  find  the  final  solutions  for  E  and  H,  it  is  necessary 
to  calculate  the  coefficients  F  ,  F„ ,  F~,  and  F,  ,  one  of  which  is 

arbitrary,  say  F  .   Using  Eq.  (38),  it  may  be  seen  that  F„  and  F  differ 
only  by  a  phase  factor.   Thus  letting 


Y   -  e     '      '- 
2> 


B-(60) 


then  Eqs.  (38)  and  (39)  become 


af 


B-(61) 


V  B-(62) 
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Taking  the  complex  conjugate  of  Eq.  (62)  and  dividing  it  into  Eq. 
(61)  produces 

(S+Ca)         ^  ^s^2\)  B_(63) 


or, 


C</cf>  ^    ^       aC3u/ 


(<t-^      (stt-.)  B-(") 


Using  Eqs.  (44),  (48),  and  (49),  Eq.  (64)  may  be  written  as 


e  e  o    e,        ai    ^      d    e  B_(65) 


Using  Eq.  (47),  the  exponentials  in  Eq .  (65)  yield  the  following  relationship: 


re  TF 

1     ~       Ti,  B-(66) 


Turning  now  to  F   and  F,  ,  Eqs.  (34)  and  (36)  can  be  solved  in  terms 
of  F  .   Using  Eqs.  (60)  and  (66),  Eq .  (34)  becomes 
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a.  v  / 


=aFae-£<^  ^(qq 


B-(67) 


B-(68) 


From  Eq.    (36) 


F.-t^F.^  +  F^] 


2 


B-(69) 


5  Lc  +e  °      °  B-(7( 


or,    from  Eq .    (47), 


F^V^^^e^-)] 


-  5Fa  e^  e-'^a   <Uw  ((p™  4-mtt 


Finally,    from  Eqs.     (2),    (60),    and    (66) 


B-C71) 


B-(72) 
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Te  -rrr 


B-(73) 


^  e       ^e^oCgx-ff^)        B-(?4) 


It  is  now  possible  to  write  down  the  TE  solutions  for  E  and  H 
by  using  Eqs.  A- (59)  through  A- (64)  together  with  Eqs.  (1)  -  (3), 
(68),  (72),  and  (74). 


t=  k  H,  +  z  U 


B-(75) 


B-(76) 


Af^.^^^e^^ 


B-(78) 


B-(79) 
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where 


E^.aF.e-A^^e—    *  x*o      ,(80) 


From  Eqs.  A-(60)  and  A- (62),  H  =  (-i/k  )3E  /3z  for  all  x;  hence  from 
Eqs.  (1)  -  (3) 


A  J-  ~  fo/M)  E  Wol^ 


3  B-(83) 


Last  of  all,  from  Eqs.  A-(60)  and  A-(64),  Hz  =  (i/kQ) 3Ey/9x;  hence 


\z  ^-l^/^jWg^  ^^3     Jv-V^x  fcQ 


B-(84) 


B-(85) 
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hJ$*  ^^s/jge  (x)        w  x^-w    b- 


(86) 


TM  Modes 


Since  the  wave  equation   [Eq.  A- (67)]  for  A  has  the  same  form  as  Eq, 
A- (56)  for  F,  the  following  choice  is  made  for  A   [A  =  x  A  ]: 


x         x 


a^(v^kmaV<* 


fjDC^O    B-(87) 


v  om^-^h^\^-^j 


B-(88) 


^=W.MtW"^e5'x    V    x^-W 


<L     W   B-(89) 


where  A  ,  A„ ,  A~  and  A,  are  constants  (possibly  complex).   Substituting 
Eqs.  (87)  -  (89)  into  Eq.  A-(67)  leads  to  the  following  equations: 


?" 


.a=  AVh,'  6_(90) 


q      or)      i 
f  ^  "     °  n^  ""  9<  B-(91) 
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fr 


n  2      a 


•«>  '5 


B-(92) 


Again  it  is  oossible  to  determine  the  relationships  among  3,3,  and  3 

by  imposing  the  boundary  conditions  on  F.  and  H.   For  TM  modes,  continuity 

— *•     _^. 
of  the  tangential  components  of  E  and  H  (i.e.,  F   and  F  )  are  given  by 


[^  1    =   I 


e. 


W.-M 


cd    K  =0 


B-(93) 


B-(94) 


fcV    [E.1 


LMr  W, 


of  K--IV 


B-(95) 


B-(96) 


However,  Eqs.  A-(72)  and  A-(74)  imply  that  the  boundary  conditions  on  Er 

2 
are  equivalent  to  the  requirement  that  (-i/k  n  )3H  /3x  be  continuous 

o     y 

across  the  boundaries.   Thus,  Eqs.  (93)  and  (95)  can  be  written  as 


B-(97) 
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_1_ 


rat 


$     L 


H*      d« 


V  (^J3  *  "-* 


B-(98) 


Using  Eqs.  (87)  and  (88),  Eq .  (94)  implies  that  at  the  n  -  n  boundary 


-^  A., 


x^o 


a  k. 


o 


B-(99) 


or,  from  Eqs.  (87)  and  (88), 


A,e42:^+Aj)rft 


R-(IOO) 


Since  A^,  A^,    and  A^  are  constants,  for  arbitrary  z  the  following  relation- 
ship must  exist: 


?-% 


B-(101) 


Similarly,  at  the  n   -  n   boundary,  the  boundary  conditions  require  that 

6       S 


if. 


-up5  Ays 


=  -c|3.A 


v*-W 


^  •jCs-W' 


B-(102) 


146 


5  ~~B-(103) 


Eq.    (103)    can  be   satisfied  only  if 


ftrP' 


B-(104) 


Again,  since  6  =  k  n  sinO  is  the  wavevector  component  along  the  z-axis 
within  the  guide,  it  is  convenient  to  let 


=  Kn^V-^-ffp. 


B-(105) 


Substituting  Eq .  (105)  into  Eqs.  (90)  -  (92),  and  solving  for  a.,, 
'  ,  and  s. ,  leads  to  the  following: 


V  \(n'^9-ntf 


U 

B-(106) 


0l  °  B-(107) 
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*-<y 


B-(108) 


Comparing  Eqs.  (106)  -  (108)  with  Eqs.  (26)  -  (29)  for  the  TE  case,  it  may 
be  seen  that 


1=0. 


>     ^r%    )    S,  =  S 


B-(109) 


Also,  Eqs.  (105),  A-(72),  and  A-(74)  allow  Eqs.  (97),  (94),  (98), 
and  (96)  to  be  written  respectively  as 


X  liiwl     _  _L  [h< 


K*l   **  lx^      *?l**Ls 


*  =  o 


B-(llO) 


M_  -  KL 


v^o 


x~o 


B-(lll) 


_L 


khsJ 


[He 


s  i  **j    ,    ^ i  sn        B_(112) 


J  r  =  -VU  J  v  =  - 


K  =  -W 


X  =  -Y\/ 


B-(113) 


Substituting  Eqs.  (87)  -  (89)  and  (109)  into  Eqs.  (110)  -  (113)  produces 
the  following  boundary  relations: 
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\    %^  ~C^\A3  B-(114) 


K     ^A=  c  f),      a  A     -oV'qA. 


A   -    A^  +  4 


~3  „  -sW 


*  ^  B-(115) 


",»V-   -V^Ase-a^v^B-a 


<0 


_2 
Multiplying  Eq .  (115)  by  +an    and  adding  it  to  Eq.  (114)  produces 


~  fca («■  "I*  + 1 *~Q  j)  -  A3(*  n:a-  m'*  g)    B. 


(118) 


_2 
Similarly,  multiplying  Eq.  (117)  by  -sn    and  adding  it  to  Eq.  (116) 

produces 


Kir^uyf^l^-s^n^-^ 


Dividing  Eq .  (118)  by  Eq .  (119)  results  in 


(-vVu^)    ^  (n+c+j^V 


(-Shf    tc*j%)  O^-Cn-4,) 


^     U^g       0       B-(120) 
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B-(121) 


Using   Eq.    (44),    Eq.    (121)    can  be   re-written  as 


!J|Wo  eaiWl(a<9/n-s)    aiW'^O^^4) 


B-(122) 


Again,  with  no  loss  of  generality,  e  1S  "  can  be  written  as  e   g  m 
where  m=0,  1,  2,  ...   .   Thus 


e*p[i  Ca8W"M-  «^>£ai(C  ^,T)1 


j    Yji  J  J        B-(123> 
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where 


^W-#CJ  -A<f£  ^^-r 


C=  ^-'[W?)KV<)] 


«C  »W'[t^  («8V»/fl 


B-(125) 


B-Q26) 


Then,    from  Eqs.    (106)    -    (109), 


fij  -  w  {(^)K^a^-n.a)y,a^i]B_(i27: 


C-  ^-'{(V^I^*8-  -^//S0^} 


B-(128) 


TM       TM  TE       TE 

It  mav  be  noted  that  <j>    and  i>        are  identical  to  i>    J   and  <p         in  Eqs. 
ag       gs  Tag       gs 

2   2       2   2 
(50)  and  (51)  except  for  the  factors  n   /n   and  n   /n   in  the  argument 

gags 

of  the  arctangent.   Using  Eqs.  (107)  and  (109),  Eq.  (124)  becomes 


™  ,TM 


a^W^o»o-^  -sup™  =  amir 


B-(129) 
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which  is  the  waveguide  equation  for  the  TM  case. 

Solving  for  W  and  following  the  TE  arguments  (which  apply  equally 
well  to  TM  modes)  preceding  Eq.  (54)  for  the  cutoff  condition,  the 
following  equations  result: 


€'-  ^  {in^l^-^^i 


B-(130) 


ft'   O 

'$S  B-(131) 


so  that  the  following  equation  is  obtained  for  (V  .  )_„: 
ft   n  mm  TM 


V  »^)7M     5Tr(^3-0'a(_         [  >V  CVlJ'-M  •)"aJ  B-C132) 


2    2 
Eq.  (132)  is  similar  to  Eq.  (59)  except  for  the  factor  n   /n 

g   a 

In  order  to  find  the  final  solutions  for  E  and  H  for  the  TM  case, 
it  is  necessary  to  solve  for  A  ,  A„,  and  A,  in  terms  of  the  arbitrary 
amplitude  A  ,  similar  to  what  was  done  for  the  TE  case.   Letting 
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^3  ^a,  B-(133) 


then  Eqs.    (118)    and    (119)   become 


-(*<  +U^  =  e6   f  a^^-iVy 


-<*^V^^ 


Taking  the  complex  conjugate  of  Eq.  (135)  and  dividing  it  into  Eq.  (134) 
produces 


;«* ^V)  «e«r  on^ctf 


(sn"fCni  k^r?) 


B-C136) 


B-(137) 
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Using  Eqs.  (44),  (125),  and  (126),  Eq.  (137)  becomes 


B-(138) 


or,  with  no  loss  of  generality, 

e 


'<-e^K+(C^^  -*8* 


B-(139) 


Using  Eq.  (124),  the  following  relationship  exists  for  the  phase  factor; 


0  B-(140) 


Turning  now  to  A±   and  A^ ,  Eqs.  (115)  and  (117)  can  be  solved  in  terms 
of  A„. 


B-(141) 


*-%)  B-(142) 
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Also, 


V-e-fv'^tA,^] 


W 

B-(143) 


^•^e-^fe-'^^/^O 


B-(144) 


or,    from  Eq .    (124), 

B-(145) 


5ft/        -£<P7W 


5rv         _  L  (J/  /        -T)U  \ 


B-(146) 


Finally,  from  Eqs.  (88)  and  (133), 


^\  tf*  -_  A   e4  <  [e!^  flA^ 


^C  'B-(147) 


-a 


-3Aae-^   Cao(8*+$ 


B-(148) 
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It  is  now  possible  to  write  down  the  TM  solutions  for  E  and  H  by 
using  Eqs.  A-(70)  through  A-(75)  together  with  Eqs.  (87)  -  (89), 
(142),  (146),  and  (148) 


Emz   K  p*  +  z  ^ 


W  3  »3 


e.M-^S1'*?** 


where 


.i<e 


-cyx 


B-(149) 


B-(150) 


B-(151) 


B-(152) 


B-O-53) 


B-(154) 


^■-^•'^c-df*-,^^  B.(156) 

for    V4-\\j 
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From  Eqs.  A-(70)  and  A-(74),  E  =  (i/k  n2)3H  /3z  for  all  x;  hence  from 
Eqs.  (87)  -  (89) 

B-U57) 


F„  W  =  (pA.Mj*)  Ctx)    ^-WiiXhO 


ok 


»«5~v       ^-yv*x±o  B_(15g) 


^cx  ^)  *  ifK*?)  H-^iy.)         hr    x^  ~W  b-cis9) 


Last  of  all,  from  Eqs.  A- (72)  and  A- (74),  E  =  (-i/k  n  )3H  /3x;  hence 


V  "  B-(160) 


for     -Wi^fefl 


^^•-i^JM     •*«-  **-h/ 


B-(162) 


APPENDIX  C 
INHOMOGENEOUS  ASYMMETRIC  WAVEGUIDES 

TE  Modes 

The  following  choice  of  trial  solutions  is  made  for  E   (where 

x 

F  =  x  F  )  : 
x 

F«=  We;(wt^2V*         -ft.   ,,-W 

This  choice  is  similar  to  the  choice  made  for  the  homogeneous  case  [Eqs. 

B-(l)  through  B-(3)],  except  that  now  a  provision  is  made  in  the  form  of 

the  solution  for  F   to  allow  for  an  index  gradient,  i.e.,  instead  of 
xg 

+X2X  — 1R 

e     the  choice  is  e    ,  where 


c-(i) 


0     C-(2) 


3«3(x) 


C-(4) 


and  where  it  is  assumed  that  all  x  dependence  of  F   is  contained 

xg 

within  g(x).   The  function  [g(x)/k  ]  is  often  called  the  eikonal. 

From  Eqs.  A-(99)  through  A-(104),  the  choice  A  =  0  leads  to  the 
following  relations  for  the  E  and  H  fields  in  the  TE  case: 


C-(5) 
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C-(7) 
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-  (i/K) 


£5l 


C-(8) 


3  C-(9) 


«,-  -6V».)£&- 


<ted*  c-(io) 

The  boundary  conditions  on  E  and  H  (for  no  currents  in  the  case 
of  H)  require  that  the  tangential  components  of  E  and  If  always  be 
continuous  across  a  boundary.   This  statement  is  true  for  any  boundary; 
no  distinction  is  made  between  a  boundary  between  two  different  media 
or  a  boundary  between  two  different  portions  of  the  same  medium, 
includ'ng  an  inhomogeneous  medium.   The  following,  argument  is  useful. 
Eq.  A- (7)  and  A- (8)  from  Maxwell's  equation  are  integrated  over  the 
surface  t  =  Ax  Ay  shown  in  Figure  C-l: 


*  J   it 

From  Stokes'  theorem  for  any  vector  U, 


C-(ll) 


C-(12) 


[  Vx  0  dx  =    §   u  >dJl 


C-(13) 


where  I   is  the  contour  path  surrounding  t,  and  where  the  symbol  <£ 
indicates  that  the  integration  is  over  a  closed  path.   Hence  Eqs.  (11) 
and  (12)  become 
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n3U) 


® 


A 


v*,1*a*  <  j 


-^  (-  /?TM) 


Figure  C-l.   Rectangular  path  of  integration  within  an  inhomogeneous 
guide  to  show  behavior  of  tangential  TE  (and  also  TM)  component.   This 
view  of  the  guide  is  one  of  looking  "down  the  guide"  (in  the  direction 
of  propagation),  which  is  different  from  that  shown  in  Figure  II-l. 
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£-dJ:-±  \  ^.dr  c.(14) 


$tf.tf .  J-  Ji^l.^ 


C-(15) 

J?  T 

Referring  to  Figure  C-l ,  it  may  be  seen  that  the  dot  nroduct  p?  '  d£  will 
have  the  following  contributions  onlv : 


j  5-dt=  e^k)^  +  t^^o)^) 


C-(16) 


where  |ETE|  and  |ET]g|  are  the  amplitudes  of  E  at  the  levels  x  =  -  x 
and  x  =  -  x  .   A  similar  equation  results  for  the  TM  case: 


$■&!--      rT^-X^-f   ^6-Of-Atj)  c-(17) 

S. 

which  is  given  here  since  the  relationship  for  the  TM  case  will  be 
useful  in  the  next  section. 

Defining  the  direction  of  t"  as  t  =  "f/ |  t  |  ,  then  the  right-hand 
sides  of  Eqs.  (14)  and  (15)  can  be  written  as 

C  J  Jt         C      £t  J  C-(18) 

t 


£S^U-±^.S  6.-3 

X 

Using  Eqs.  (16)  -  (19),  Eos.  (14)  and  (15)  become 


C-(19) 


C--(20) 


161 


U;u)-  O-o]^,  x  A^ffh  Ak6 


C-(21) 


Since  Ay  is  arbitrary  and  is  common  to  both  sides  of  Eqs.  (20)  and  (21). 
then  in  the  limit  as  Ax  and  Ay  approach  zero,  Eqs.  (20)  and  (21)  become 

Assuming  3H/3t  and  3 (n  E/3t  as  well  as  H  and  n^E  remain  finite,  then 
the  right-hand  sides  of  Eqs.  (22)  and  (23)  go  to  zero  as  Ax  cr0es  to 
zero.   In  the  limit,  then, 

C-(24) 

ht;k)-  ^  W=o  c-(25) 

Since  -  x„  is  an  arbitrary  position  within  the  guide,  then  E   and  H 
(which  are  the  only  tangential  components  to  the  x-axis)  are  continuous 
within  an  inhomogeneous  guide. 

Because  of  this,  Eqs.  (2)  and  (24)  imply  that 


LU)t 

e 


At  the  level  z  =  0,  then  the  following  relationship  must  exist: 
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3kH^°) 


C-(27) 


i.e.,  g(x)  is  continuous  within  an  inhomogeneous  guide. 

It  is  now  useful  to  determine  the  relationships  among  8,3,  and 

a   g' 

0   using  the  boundary  conditions  on  E  and  H.   These  conditions  are 
given  by  Eqs.  B-(12)  through  B-(15),  and,  because  of  Fqs.  (6)  and  (10), 
by  Eqs.  B-(16)  and  B~(17)  too.   Substitution  of  Eqs.  (1)  -  (3)  into 
Eq .  (6),  together  with  the  boundary  conditions  given  by  Eqs.  B-(12) 
and  B-(14),  produces  the  following  equations : 


;/a»0|       =  (f%)\ 


na  n3  rv  ^         c_(2g) 


At  the  n  -  n  boundary,  Eqs.  (1),  (2),  and  (28)  produce 
i€  "e         L<3  3       J    C-(30) 


Since  F  ,  F?  e  .  and  F  e  are  constants,  then  for  arbitrary 
z  and  t,  the  following  relationship  must  exist  in  order  that  Ea .  (30) 
be  satisfied: 


(&.-/». 


C-(31) 
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Similarly,  at  the  n   --  n   boundary,  Eqs.  (2),  (3),  and  (29) 


produce 


<^-M  _iW    ;K-aFe;,H+F<-'/-«| 

f  <s-e(.a  3  J   c-(32) 


Eq .  (32)  can  be  satisfied  only  if 


?*~h 


C-(33) 


Since  3   is  the  wavevector  along  the  z-axis, 
g 


P*  =  ^°  ^GO^BC*)  C-(34) 


But  bv  assumption   all  the  x  dependence  of  F   is  contained  within 

xg 

g(x) ,  hence 


A  Cx)  =  C^W"  C-  ( 35 ) 

1  J 


It  is  convenient,  then,  from  Fqs.  (31)  and  (33)- (35)  to  let 


3  e  Jfe0  ru  6<)s^8&)  =fl3=  ftjA=  OuWr     c-o6) 


J     n  C~(37) 
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where 


m  fift) &*&(*)=- cjnsbuit 


C-(38) 


which  is  simnlv  another  way  of  expressing  Snell's  law. 

Returning  now  to  the  wave  equation,  Fqs.  (1)  -  (3)  are  substituted 
into  Fq .  A- (92),  producing 


C-(39) 


-(S'+s'+^^Jfi^O  C-«D 


where  g'  and  gr  are  the  first  and  second  derivatives  of  g  with  respect 
to  x.   Since  either  F„  or  F  can  be  set  equal  to  zero,  then  Eq .  (40) 
can  be  re-written  as  two  equations: 

F„e'3fV- ,    ,     . 

C-(42) 


C-(43) 


The  assumption  is  now  made  that,  for  a  slowly  varying  medium 
(i.e.,  for  a  refractive  index  that  varies  slowlv  with  x) ,  g';  is  small. 
(This  assumption  will  be  discussed  in  more  detail  below.)   Thus,  to  a 
first  approximation 
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Q     »  0  C-(44) 

For  this  situation,  Eqs.  (39)  -  (A3)  become  [with  the  aid  of  Eqs.  (36) 

2 

;h  the  trigonometric  identity  cos" 6 


2  2 

together  with  the  trigonometric  identity  cos" 6 =  1  -  sin  6] 


±  kfCHf^B-**) 


a 

U  C-(45) 

Mj  '       °    ^  C-(46) 


C-(47) 


Thus,  using  Eq.  (37)  in  Eqs.  (45)  -  (47)  results  in  the  following  relations; 

C-(48) 


'/a 


C-(49) 


C-(50) 


where 


<l603  njx)  *&>&(*)  c-(5D 


Integrating  Eq .  (49)  produces 


gM  =  A.  I  3d*  c-(52) 


=  A0  j  ^O^  <**W  ^  C-(53) 


'a 

6 
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It  may  be  noted  from  Fq .  (53)  that,  for  the  homogeneous  case  (n   cos P = 

constant),  g(x)  =  k  (n  cosG)x,  as  before;   further,  g"  =  0  is  then  an 
o   g 

exact  relation. 

Using  the  value  for  g'  given  by  Fq .  (49),  then,  to  a  second-order 
approximation t 

//_ 


3  =  M 

Substitution  of  Eq .  (54)  into  Eqs.  (42)  and  (43)  produces 


C-(54) 


3'«fcV  *iA.?0 


C-(55) 


^\%ti±^/K{) 


A* 


C-(56) 


where  the  +  and  -   refer  to  the  F2e  "  and  the  F_e    terms   respectively. 
Since  the  right-hand  term  in  the  bracket  of  Fqs.  (55)  and  (56)  is  small, 
then  g'  may  be  written,  to  a  good  approximation  as 


hence 


^\%0±  <t'/jVJ 


a'*  K%  tlv*i 


C-(57) 


C-(58) 


cix  C"(59) 
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Integrating   Eq.     (59)    produces 


qW;    6fcj    <*(y}  dj<     +     t^(<g 


C-(60) 


Using  Fa.  (60),  the  following  may  be  said 

X 


cr  J<K     r-   U    -'/air-   rA.Jtf**       :A  \  qf ct^  n 


C-(61) 


The  effect  of  the  second  approximation,  then,  is  to  cause  the  amplitude 

-1/2 
to  vary  by  the  factor  q     .   It  is  now  convenient  to  use  Eq .  (52)  for 

the  value  of  g(x),  writing  the  imaeinarv  term  in  Fq .  (60)  explicitly  as 

an  amplitude  factor.   Thus,  g'  =  k  q  and  g'  =  k  q'.   In  this  way,  the 

function  g(x)  will  be  real  and  thus  will  describe  the  phase  only  of  F 

From  Eq.  (?)  [and  Eq .  (36)],  then 

-'fa/    ,     v     Hut-Az) 


yrw-  ?  h^+^] 


C-(62) 


It  is  now  necessary  to  determine  under  what  conditions  Fq .  (62) 

is  a  good  anproximation  to  the  solution  of  the  wave  equation.   Since, 

ig      —is 
as  mentioned  previously,  e   and  e    are  orthonormal  functions,  then 

each  one  must  independently  satisfv  the  wave  equation  Fq .  ,A-(84).   Taking 

-ig 
the  first  derivative  with  respect  to  x  of  the  two  terms  e    in  F   leads 

to  the  following  two  relationships: 
ix    l  1  '"\   "3/LJ   '3,3   -    J*-  C-(63) 


C-(64) 
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±12; 

where,  for  the  sake  of  brevity,  the  term  F„  „e    represents  the  two 

is        ~ig 
orthonormal  terms  F9e   and  F  e    .   The  second  derivative  with  respect 

to  x  of  F   is  given  hv 
xg 


where  the  factor  in  braces  in  Fq .  (66)  is  valid  for  both  of  the  terms 


±1° 

F   e   ".   Further, 


C-(67) 


Substituting  Fqs.  (66)  and  (67)  into  Eq .  A-(92)  produces 

2    2   9  2   2    2 

or,  since  -(3+kn~   =  k.n   cos  6,  then 

o   g  o   g 


C-(68) 


C-(69) 


In  order  that  q  satisfy  the  definition  in  Fq .  (51),  then  in  Eq.  (69) 


fl.%al>^)  4sf/'f** - **/**)  -<*» 


1  -  3  H 
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«l 


C--(71) 


The  method  described  above  for  solving  the  wave  equation  for  an 

inhomogeneous  medium  is  called  the  WKBJ  method,     and  Eq .  (71)  is 

the  WKBJ  approximation  for  the  TE  case.     The  solution  given  by  Eq .  (62) 

[and  Eq.  (52)],  when  substituted  into  Eqs.  A-(IOO),  A-(102),  and  A-(104), 

is  the  TE  WKBJ-type  solution    for  an  inhomogeneous  medium.   It  may  be 

seen  that  for  large  wavenumbers  k   (i.e.,  small  X  ),  the  WKBJ  approximation 

o  o 

is  more  easily  satisfied.  However,  as  q  =  n  cosO  approaches  zero,  Eq. 
(71)  cannot  be  satisfied,  regardless  of  how  small  q'  and  q"  are.  This 
situation  will  be  discussed  in  Appendix  D. 

Proceeding  now  to  the  boundary  conditions,  substitution  of  Eqs.  (1), 
(3),  (36),  and  (62)  into  Eqs.  B-(12)  through  B-(17)  [or,  equivalently , 
Eqs.  B-(30)  through  B-(33)]  produces 


C-(73) 


C-(74) 


5) 
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From  Eq .  (52)  , 


C-(76) 


C-(77) 


Thus  Eqs.  (72),  (73),  and  (75)  become 


-<t.F.  = 


M*^fe-FJ 


-V, 


i\    ®%®[rv*\] 


C-(78) 


3- 


C-C79) 


')[F,e*>   H=« 


5 


3/ 


/Ac^aM. 


MWO))^"^-* 


The  assumption  is  now  made  that 


Vs 


4/%))»iirawiw 


C-(81) 


C-(82) 


which  is  compatible  with  the  WKB.T  approximation  in  Eq .  (71).   Using 
Eq  .  (82),  then  Eqs.  (79)  and  (80)  become 
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C-(83) 


C-(84) 


Eqs.  (78),  (83),  (74),  and  (84)  are  very  similar  to  Eqs.  B-(34) 
through  B-(37);  hence  these  equations  can  be  solved  simultaneously  in 
similar  fashion.   Multiplying  Fq.  (78)  by   +a  and  adding  to  Eq.  (83) 
results  in 


F. 


V: 


b%%)  +  «•*.  %W%b  i"^-^t4 


C-(85) 


-R  JV+c^o 


=  F3[a-  tA.^oj] 


C-(86) 


Multiplying  Eq .  (74)  by   -s  and  adding  to  Eq .  (84)  produces 

-F9[sf,/9(-W)-^4,,3M]e;9M 


C-(87) 


-^-^HK^^b^-*)] 


C-(88) 


Dividing  Fq .  (86)  by  (88)  results  in 


-a:3H       [o.-iKisMl     Ls-^.jj(-w)] 


C-(89) 
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Using   Fq.    B-(44),    Eq .     (89)    becomes 


Tf  ^r<r 


0       '  Ta3  > 


where 


and,    from  Eqs.     (48),    (50),    and    (51) 

2       ..  5  \'/a 


6L 


*r~  'C  > 


C-(90) 


C-(91) 


C-(92) 


C-(93) 


A^6)  *)a(0)cw^C6j  C-(94) 

A?H  Klj-WJQW^-Wj  C~(95) 


Using  Fq .  (53),  Fq .  (91)  becomes  the  waveguide  equation: 

o 


It  may  be  noted  that  when  n   (and  thus  cos  0  too)  is  constant  within 
the  guide,  then  Eq .  (96)  reduces  to  the  waveguide  equation  for  the 
homogeneous  case,  Fq .  B-(52).   Using  Eq .  (38)  the  integrand  in  Fq .  (96) 
can  be  written  as 


3 

Also, 
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C-(97) 


«.tw  "      b{^)-f] 


m 


C-(98) 


J-*.5 


M^        Lrfw-fl'' 


Yg~  C-(99) 


F.q .     (96)    becomes 

\    K  M  -)  J       -    3-  U*  f  /  ry 


Thus,  F.q .  (96)  becomes 

3k 


2 
'5 


C-(IOO) 


It  may  be  seen  that  Fa.  (100)  is  transcendental  in  both  rand  W. 
Writing 

d  C-CL01) 


then,  for  a  given  value  of  n  (0),  perhaps  the  easiest  way  to  solve 
Ea.  (100)  is  to  pick  a  value  for  0(0),  then  try  values  of  W  until  one 
is  found  which  satisfies  Eq .  (100).   By  using  different  values  of  6(0), 
a  set  of  curves  can  be  generated  for  the  different  modes  m. 

Fa.  (96)  [or  Fq .  (100)]  can  be  used  to  obtain  the  cutoff  thickness 

(VT   )    of  a  guide.   There  are  two  possibilities  for  waveguide  cutoff. 

min  TF 

The  critical  angle  at  the  n   -  n  boundary  is  defined  by 

g    s 

9.     s  4^~'  [**/"*] 

°c^s  V  C-CL02) 
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while  the  critical  angle  at  the  n  -  n  houndary  is  given  by 

a    g 


%=^"'[^K1 


C-(103) 


where 


Now ,  i  f 


«    =«a(-w 


Kf-^) 


M~v>)*  ©« 


0(104) 


0(105) 


C«S 


0(106) 


0(C)  4  %.?  0-(107) 

then  the  light  waves  will  escape  or  leak  out.   From  Fqs.  (102)  -  (105) 
the  cutoff  conditions  become 


s^  BM  £     ■■ 


Hcs  0(108) 


1         n. 


0(109) 


M}(-^W-*)ins  c_(110) 
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r\  (o)^ie(o)  ~  n^  c_(111) 

But,  by  Eq.  (38),  F.qs.  (110)  and  (111)  become 

V\   =  Y\^(K)  ^&X)    =   ns  C-(112) 

^  ^  K1  (x)  W\£&)  =  n^  C-(113) 

It  was  assumed  for  homogeneous  waveguides  that  n  >  n  .   Hence  as  9 
decreases  from  a  maximum  value  of  tt/2  (propagation  straight  down  the 

guide),  then,  for  a  given  functional  form  of  n  ,  the  cutoff  condition 

first  to  be  met  is  given  by  Eq.  (112).   That  is,  waves  will  start  to 

leak  out  at  the  n   -  n   boundarv  when  sin  G(x)  =  n  /n  (x) :  they  will 
g    s  s   g 

still  be  totally  internallv  reflected  at  the  n   -  n   boundary.   This 
result  is  identical  with  the  result  obtained  for  the  homogeneous  case. 
Substituting  Eq .  (112)  into  Eqs.  (92)  -  (95)  and  Eos.  (98)  -  (99) 


TE 
produces  6    =0  and 
gs 


"^'^-•m-oi  „m 


Hence  at  cutoff,  Eq.  (100)  becomes  the  following: 

V.  "»ru*i  /ye 

For  the  inhomogeneous  case,  the  situation  ng  <  ng  is  different  than 

n   >  n   (since  n   4   n   ).   Hence  if  n   <  n  ,  then  the  cutoff  condition 
s    a         ag    2S              s    a 

TF  TF 

is  eiven  by  En.  (113),  in  which  case  6        4   0,  d>   =  0,  and  Eq.  (100) 


becomes 


h^-^^l^-i^j 


-00™ 


C-(116) 
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In  order  to  find  the  final  solutions  for  E  and  H,  it  is  necessary 

to  calculate  the  coefficients  F. ,  F„ ,  F_,  and  F, ,  where,  as  before,  F„ 

1   z    3       4  z 

is  chosen  as  the  arbitrary  amplitude  coefficient.   Letting 


C-U17) 


;a^)" 

F  *  e 

h 

3 

"ji 

Then  Fqs.  (86)  and  (88)  become 

ere 


C-(118) 


-[s-cJ^jMje2^   ^c^  fs^cJp^-K/j]   c-(ii9) 


Dividing  the  complex  conjugate  of  Eq .  (119)  into  (118)  (and  re- 
arranging factors)  produces 

eiLf^\      a^WoJ        S'^.tf'*)     efltV'^ 


0.  -   -h0%(o)  5tdo  «(-*) 


Using  Fqs.  B-(44),  (92),  and  (93),  then 


C-(120) 


C-(121) 


f<p  rc"  =  -  s  ^  +  p  <p47  *  *j M+amir       c.(122) 


Using  Fq.  (91),  Eq .  (122)  implies  that 


<r-  - <p. 


C-(123) 
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Turning  now  to  F  and  F, ,  Eqs.  (78)  and  (74)  can  be  solved  in 
terms  of  F? .   Using  Fas.  (117)  and  (123),  Eq.  (78)  becomes 


yn=\ 


F,=  f>)F30.e"^) 


-'fc,    ^   -.     ^<C    .  ,  J- 


C-(124) 


=  g    ^)  FS  *""  ^  c3  <>*(&-  J  C-(125) 


Also,  from  Eq .  (74) , 


sW  -«/a/  \  c    r*  ^W     ^'c5 


\TZF 


^M"'WF»[tf  te'    ^e^J 


C-(126) 


=e    I 


7       /  rc  \  2  C-(127 


or,  from  Fa.  (91) 


*       Q   l   '  A        <-  ^  e  -i    C-(128) 

^SF-S5^^'1^    C*0  (fP«   *w)  C-(129) 


a 


Finally,  from  Fqs.  (62),  (117),  and  (123) 


__  C-(130) 
C-(131) 
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-1/2         -1/2 
It  may  be  noted  that,  except  for  the  factors  q     (0)  and  q     (-W) , 

Eqs.  (125)  and  (127)  are  identical  in  form  with  Eqs.  B-(68)  and  B-(72) 

TE      TE 
providing,  of  course,  that  the  appropriate  values  of  d>   and  <!>*  "  are  used. 

ag      gs 

Eq .  (131)  is  similar  to  Fq .  B-(74),  with  the  exceptions  that  the  term  gx 

is  replaced  by  g(x)  for  the  inhomogeneous  case:  also,  there  is  the 

-1/2 
additional  factor  q     (x) .   Further,  since  Eqs.  (5)  -  (10)  are  identical 

to  Eqs.  A- (59)  through  A- (64)  for  the  homogeneous  case,  it  is  possible 

-1       _* 

to  write  down  the  solutions  for  F   and  K   for  the  inhomogeneous  case 

I  L        J  t. 

by  using  Eos.  B-(75)  through  B-(86),  providing  of  course  that  the  term 

gx  is  replaced  by  g(x),  and  also  that  the  factor  (g/k  )  in  Eq.  B-(85) 

is  replaced  bv  [g'(x)/k  ]  =  q (x) .   [There  is  an  additional  term  for  H 

o  oz 

-1/2 
due  to  the  x  derivative  of  q     (x) ,  as  given  by  Fa.  (10).   However,  according 

to  Eq.  (81),  this  term  is  negligible  in  the  WKBJ  approximation.]   Thus, 

E   ,  H   ,  and  H   for  the  inhomogeneous  are 
ov   ox       oz 

y,  7j  L°J    ae    -w^r^yc    -frr  x?0        c-(i32) 


C-(133) 


C-(134) 


^)-  H'^  Fae'^  ^Cxl  +  <J^-W*^o 

UJx)  -    "   (f3/^)^Cx)             V  aSH  x  c_(135) 

ttol(x)*  -c(VA.Je03(x)        W     X^O  c-(i36) 

t\M-    -^WfcM,[gWf^]^(xJ  ^_W.X,0  c_(i37) 

^»6c)-  flfs/A^^W          W     X^-W  C-C138) 
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TM  Modes 
The  following  choice  of  trial  solutions  is  made  for  A^  (where 
A  =  x  A  ) 


A.   * 


^/i^e^'Ml 


-0L* 


-fev  X  ^o 


c^wt-^Oi 


V  ^W*       t^+A3e   J  V-^a^O 


K>~-(*Ji&^z) 


5X 


-W       x  =  -H/ 


C-(139) 


C-(140) 


C-(141) 


From  Eqs.  A-(99)  through  A- (104),  the  choice  F  =  0  leads  to  the  following 
relations  for  the  E  and  H  fields  in  the  TM  case: 


E    -[cA.n^)]^S 


az 


C-(142) 


V0 


JaA. 


er-bH.nMzZ 


C-(143) 


C-(144) 


Ux  =  c 


C-(145) 


H     a       d*. 


fl  *2- 


C-(146) 


tt2-  0 


C-(147) 


It  is  useful  to  first  determine  the  relationships  among  B  ,  B  , 

d   g 

and  6   using  the  boundary  conditions  on  F  and  H  given  by  Fqs.  B-(93) 
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through  B-(96).   Using  Eq .  (146),  then 

•f*  fl»)  Iryn.        *     (&"'») 'vn,  c"(149) 

Using  Eqs.  (139)  -  (140),  Eq.  (148)  and  (149)  become 


C-(150) 


-(151) 


Eqs.     (150)    and    (151)    can  be    satisfied    for   arbitrary   z   only  if 

(*■'(%  C-(152) 

f>s=A\  C-(153) 

Because  Hm,  (and  therefore  A  )  must  be  continuous  within  the  guide, 
TM  x 

even  if  inhomogeneous  [as  given  by  Eq .  (32)  for  two  arbitrary  adjacent 

points  x_  and  x..  within  the  guide],  then  6  must  be  a  constant  within 

the  guide.   It  is  convenient  once  again  then  to  let  @  be  defined  by 
Eq.  (36). 

Substituting  Eqs.  (139)  -  (141)  into  the  wave  equation  Eq.  A-(93) 
produces 


(-(»%«.*  fAJV)^-*? 


C-(154) 


is: 


-/S'+s'+A.VX^o 


C-(155) 


C-(156) 


Solving  Eqs.  (154)  and  (156)  leads  to  the  results  given  in  Fqs.  (48) 
and  (50).   Again,  since  e  8  and  e  lg  are  orthonormal  functions,  Eq. 
C-(155)  can  be  re-vnritten  as  two  equations: 

-f>  -Cg'J  +t3  "  — v*-^.  "4  =o 

1  5  C-(157) 


f»    I 


0 


C-(158) 


To  a  first  approximation 


r)r  v  C-(159) 


When  this  is  done,  both  Eqs.  (157)  and  (158)  have  the  result  given  by 
Eqs.  (49)  and  (52)  for  g'(x)  and  g(x).   Using  Eq.  (49),  then,  to  a 
second  approximation 

3"-^' -*.?'- ^4^ 

3  "H 


and  then  from  Eqs.  (157)  and  (158) 

C-(161) 


n      [Ax  dx 


C-(162) 


C-(163) 


C-(164) 


hence 


X 


C-C.165) 


Thus, 


Aoe  3  +  a  e"3=n 


-'/. 


^ 


■/v.  ^4^ „  -*4s* 


C-(166) 


which  is  a  result  similar  to  Eq .  (61)  for  the  E  vector  in  the  TE  case 

except  for  the  additional  factor  n  .   Again  it  is  convenient  to  use  Eq . 

g 

(52)  for  the  value  of  g(x),  letting  the  imaginary  part  in  Eq.  (165)  be 
explicitly  written  as  an  amplitude  factor.   Thus 


It  is  now  necessary  to  determine  under  what  conditions  Eq .  (167)  is 

a  good  approximation  to  the  solution  of  the  wave  equation.   Since  either 

ifi       —  ig 

A„  or  A_  can  be  set  equal  to  zero,  then  the  two  terms  A„e  3  and  A  e 

must  independently  satisfy  the  wave  equation.   Taking  the  first  derivative 
with  respect  to  x  leads  to  the  following  two  relationships 


^~-{^i^M^\^V^} 


C-(168) 
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HP 

where,  for  the  sake  of  brevity,  the  term  A   e  b   represents  A~ 


e+ig  and 


A  e   ,   The  second  derivative  is 


±c< 


C-(169) 


Also,  the  third  term  in  the  wave  equation,  Eq .  A- (93),  is 


n^        ^x 


C-(170) 


Adding  Eqs.  (169)  and  (170)  results  in 


where  the  factor  in  braces  in  Eq .  (171)  is  valid  for  both  of  the  terms 


A   e  lg.   Further 


3     .2 

x3         "0—     *  "«} 


C-(172) 


Substituting  Eqs.  (171)  and  (172)  into  Eq .  A-(93)  produces 


**iW-*r-'$'H#*s^, 


3) 
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In  order  that  q  satisfy  the  definition  in  Eq.  (51),  then  in  Eq.  (173) 


*.Yl » 


&-  £ '  3(? 


2  " 


C-(174) 


0 


3$f-3H 


v 


nc 


2  1*. 


l\* 


«\ 


C-C175) 


This  is  the  WK.BJ  approximation  for  an  inhomogeneous  medium  for  the  TM 

case.   It  has  the  two  additional  terms  with  n'  and  n"  that  TE  approxi- 

g      g 

mation,  Eq .  (71),  does  not  have.   The  solution  given  by  Eq.  (167) 
[and  Eq.  (52)]  is  the  TM  WKBJ-type  solution.     As  q  approaches  zero, 
Eq .  (175)  cannot  be  satisfied,  and  this  situation  will  be  discussed  in 
Appendix  D. 

Proceeding  now  to  the  boundary  conditions,  substitution  of  Eqs. 
(139),  (141),  (36),  and  (167)  into  Eqs.  B-(90)  through  B-(98)  [or, 
equivalently ,  Eqs.  B-(llO)  through  B-(113)]  produces 


K-  ^~ 


)^)[*/3%e-^)] 


C-(176) 


C-(177) 


-'/ 


,^^e-^l<KM{^k'(<^-%^ 


U(-ti)  ~c<*("V)l  185 


+- 


<^W^i^-Y^jHfr^]h«^+vj^ 


C-(178) 


V^=  VV^/^W^     c-( 


-(179) 


The  values  of  g(0),  g'(0),  g' (-W)  are  given  by  Eqs.  (75)  -  (77).   Also, 
the  assumption  is  made  that 


V 


ju'jj  i«;    ^W? 


■x)/      C-(180) 


Since  g'(x)  =  k  q(x),  then  Eq.  (180)  becomes 

Jl<2      I  ,,  n 


4- » 


n 


a 


/ 


n 


^-    -  -4— 


5^  °"h? 


C-(181) 


C-(182) 


which  is  the  same  order  of  magnitude  (or  better)  approximation  as  the 
WKBJ  approximation  given  by  Eq .  (175).  Using  Eq .  (181)  together  with 
Eqs.  (75)  -  (77),  Eqs.  (176)  and  (178)  can  be  greatly  simplified: 

C-C183) 

"^-"y^'v'Hc-o,, 


Using  Eq.  (75),  Eq.  (177)  becomes 


-'/. 


A,=  "/o)?-    (o)  [Ag+A2] 


C-(185) 


_2 
Multiplying  Eq.  (185)  by   +an   and  adding  it  to  Eq.  (183)  produces: 


C-(186) 


-&  [  a  r>>)/KUa  *-i  A.  J  ty]  -  A3  [a.  nfflnS\l(»)]  c-(187> 


Multiplying  Eq .  (179)  by   -sn   and  adding  it  to  Eq.  (184)  produces 


-A. 


C-(188) 


V 


:^^8-i*.|(-4^^l5^^^ 


C-(189) 


Dividing  Eq.  (187)  by  Eq.  (189)  results  in 


")%a  +fM^)  f-:3(-»)_  Mfc)h?-  ApM-») 


C-(190) 


187 


".'WV+Aj'0)     infH/nfrKtf-u)  c-(i9i) 


Using  Fq.  B-(44),  Eq .  (89)  becomes 


4-  a.h^'ls^l-ri/rt^K^-^n 


C-(192) 


-*$•») -W™-W™*amT 


C-(193) 


where 


C-(194) 


C-(195) 


Using  Eqs.  (48),  (50),  and  (51),  then  the  arguments  of  the  arctangents 
in  Eqs.  (194)  and  (195)  become 


ft  f\  *(o)        „       j^Yg)      (7   ~^ 


a\'/j 


o/^jfo)  V  *%(Q)tK)&(0) 


S^(-w) 


sn'/i 


¥-*)     ^-*/J 


's9\<5^^  ^s2        n,(-w)  <*>&£- w) 


C-(196) 


C-(197) 


TE 


TM      TM  TE 

It  may  be  noted  that  ;    and  i    are  identical  to  :    and  :    in  Eqs, 
ag      gs  ag      gs 

9       2       2       2 
(92)  -  (95)  except  for  the  factors  n~  (0)/n   and  n   (-W)/n   .   Using 

Eq.  (53),  Eq.  (193)  becomes 
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O 


'"""        C-(198) 
-  *l 


^j  *fi  *"&*)*«    "    Sf«q  ~  ^S  ^"^ 


which  is  the  waveguide  equation  for  the  TM  case.   Again  it  may  be  noted 
that  when  n  and  cos  0  are  constant  within  the  guide,  then  Eq.  (198) 
reduces  to  the  waveguide  equation  for  the  homogeneous  case,  Eq .  B-(129), 
Using  Eqs.  (97),  (194)- (197),  Eq .  (198)  becomes 


-  £mfr~  c-(199) 


It  may  be  seen  that  Eq.  (199)  is  transcendental  in  r\   and  W,  similar  to 
Eq.  (100). 

Following  the  arguments  preceding  Eq.  (114),  the  cutoff  condition 
is  given  by  Eq.  (112)  or  Eq.  (113).   Hence  at  cutoff  for  n  >  n  , 
Eq.  (199)  becomes 


For  n   <  n  ,  Eq .  (199)  becomes  at  cutoff 
s    a 

0 


Eqs.  (200)  and  (201)  are  similar  to  Eqs.  (115)  and  (116)  for  the  TE  case. 

In  order  to  find  the  final  solutions  for  E  and  H  for  the  TM  case, 

it  is  necessary  to  solve  for  A..  ,  A_,  and  A,  in  terms  of  the  arbitrary 
amplitude  A  .   Letting 
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A  l*f* 

h^    ~~&  A^  0(202) 


then   Eqs.     (187)    and    (189)   become 


C-(203) 

pi 

5  C-(204) 


Taking  the  complex  conjugate  of  Eq.  (204)  and  dividing  it  into 
Eq .  (203)  produces 


Using  Eqs.  B-(44),  (194),  and  (195),  then 

^      ai(.w)     -al<£      ♦.*# 


^--  a3(-w)-  o^  +-£cp™  +amir  c_(208) 


Using   Eq.    (193),    Eq.    (208)    implies   that 


?'-< 


C-(209) 
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Turning  now  to  A  and  A,,  Eqs.  (185)  and  (179)  can  be  solved  in 
terms  of  A^      Using  Eqs.  (202)  and  (209) 

A,-  r\      I     (0)  *a  (  I*  e        ^j 


C-(210) 


SS»V,<5'%)A,«"*»  oo(^)  c-(2ii) 


Also,    from  Eq.    (179) 


A,-  e     n 


-       -*WAie-Oe:WJ 


or,    from  Eq .    (193) 


-itgM^J] 


C-(212) 


\T7W 


At=esWrt  s^'^-w)/Vae"<'<fi3  <W<?J  ^)        c-<2i3> 


Finally,    from  Eqs.    (167),    (202),    and    (209) 


C-(214) 


It  is  now  possible  to  write  down  the  TM  solutions  for  E  and  H  by 

using  Eqs.  (142)-(147),  together  with  Eqs.  (139)-(141),  (211),  (212), 
and  (215). 

— k        A         K     _ 

CTM        *       &  C-(216) 

TM   "   t)   3  C-(217) 


Z     '  C-(219) 

.      l(uit-Az) 

VK'2')s  V  c-(22o) 
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C-(218) 


where 


*#  3      0  *  uc/  0J  c-i 


***>)*  ~5sii"  ^°)/Le"t  *a  <WC)  c"^  w  *±o 

0  oy  <?  c-(22i) 

j 
'3W  J"  '   0     '   ■       «  '  "-C  •'fl-'  C-(222) 

^-an^^J^e^cw^tmr)^   <*-n/  c.(223) 


From  Eq.  (142),  E  =  [i/k  n2 (x) ] 9H  /9z  for  all  x;  hence  from  Eqs.  (220)- 
(223) 

ox^  '    V^  /    3  C-(224) 


e,x  (x) ,  [jft/Ji.  n87x)]  W.a(x)        w   -w±x±o 


C-(225) 


C-(226) 


From  Fq.  (144),  F  =  [-  i/k  n2(x)]9H  /9x  for  all  x,  so  from  Eqs.  (220)- 
(223)  [recalling  that  g'(x)/k  =  q(x)], 

*  C-(227) 
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^  ;>)=  ~L  (*/****)    KM  kr    X±-W  C-(229) 


-1/2 
where,  in  Eq.  (228),  terms  of  order  d[n  (x)q     (x)]/dx  are  neglected 

according  to  Eq.  (180)  and  the  WKBJ  approximation. 


APPENDIX  D 
TOTAL  INTERNAL  REFRACTION  IN  INHOMOGENEOUS  WAVEGUIDES 

Introduction 

In  Appendix  C  it  was  shown  that  the  WKBJ  solutions  for  E  and  H 

are  valid  as  long  as  Eqs.  C-(71)  and  C-(175),  for  TE  and  TM  modes 

respectively,  are  satisfied.   When  either  Eq.  C-(71)  or  C-(175)  is 

not  satisfied,  it  may  be  seen  from  Eqs.  C-(132)  through  C-(135)  and 

Eqs.  C-(221)  through  C-(226)  that  at  least  two  field  components  (E 

and  H   for  the  TE  case,  H   and  E   for  the  TM  case)  become  infinite 
ox  oy      ox 

-1/2 
as  q  ->  0  due  to  the  factor  q    .   Since  this  is  physically  unacceptable, 

some  other  approach  must  be  used  to  solve  the  wave  equation  near  the 

level  q  =  0.   In  the  field  of  radio  ionospheric  physics,  Budden   and 

12 
Wait   have  developed  a  useful  approach  for  describing  the  propagation 

of  light  near  q  =  0.   This  method  will  be  adapted  to  optical  waveguides. 

First  of  all,  it  should  be  mentioned  that   -W  is  now  chosen  to 
be  the  value  of  x  for  which  0  =  tt/2  (and  q  =  0)  for  a  given  mode  m  and 
incident  angle  0(0),  as  shown  in  Figure  IV-1.   That  is,  at  the  level 
x  =  -W,  0(-W)  =  tt/2,  q(-W)  =  0,  and  total  internal  refraction  occurs. 

The  assumption  is  now  made  that,  for  a  slowly  varying  index  profile, 

q   is  approximately  linearly  proportional  to  x  near  the  level  q  =  0, 

i.e.,  near  x  =  -W.   This  is  nearly  the  same  thing  as  saying  that  n  (x) 

g 

is  linear  in  this  region,  as  may  be  seen  from  the  following  argument. 
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194 


It   is   convenient   to   let 


n5(x)=  n0+    bW"» 


D-(l) 


where 


rt3(°)-  no  +-^n 


at  the  top  of  the  guide,  and 


D-(2) 


n.C-t>)  =•  no  d-(3) 


where  -D  is  the  penetration  depth  of  the  index  profile  such  that  beyond 
x  =  -D  there  is  no  longer  an  index  gradient.   For  instance,  n  might 
be  the  index  of  the  film  before  ion  implantation,  and  An  the  maximum 
index  change  (where  An  <<  n  )  which  is  generally  a  measurable  quantity. 
The  function  b(x)  thus  gives  the  index  profile,  where 


b(0)=  I  D-(4) 

b(-D)*0  D~(5) 


Expanding  n  (x)  in  a  Taylor  series  expansion  around  x  =  -W,  produces: 

iy*)  =  n3(-w)  +•  (x4wj  r^(-w)    +    •••  D.(6) 

~    n04   An[M-w)  4-b((-Vt)[x+w)]  d-(7) 
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which  is  linear  in  x.   For  small  An 

ri^Me-  .vf-an.  An  [U-w)  +  b'(-»)(u-»)l       D.(8) 

At  the  level  q  =  0  (and  x  =  -W) ,  9  has  the  value  tt/2,  so 

=      n      (~Wj  D-(10) 

-   r)  A4.   £n04n    fcf-WJ  D_(11) 


hence    [from  Eq.    C-(97)] 

tfM^l^V-rj7  ]    ^    <3*0  A*   b'h)(K+w)       D-(i2) 


which  is  linear  in  x.   Eq.  (12)  is  equivalent  to  a  Taylor  expansion  of 

2 
q  (x)  in  the  region  near  x  =  -W. 


TE  Modes 

2  11 

Using  the  form  of  q   given  by  Eq.  (12)  and  the  substitution 

D  D-(13) 


r-f^J 


where 


(<?2)^  art6AwJ>Y^ 


D-(14) 


«*" 


)-o 
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D-(15) 


then 

|<5n0AnJ>f-W)J  7  '  D-(16) 


y/3 


D-(17) 


and 


^w^^^r6      D.(18) 


D-(19) 


The  wave  equation  A- (92)  can  be  written  as 


D-(20) 


*  D  0        0  0        0 

since  all  z  dependence  of  F   is  of  the  form  e     ,  and  (-8  +  k  n  )  =  q  . 
By  the  chain  rule 

-zf>      *$  d^  vw  D-(21) 


^  ;  D-(22) 


Hence,  Eq .  (20)  becomes 


r  ($'J       "f-  A3  ^  P,  -  O  D-(23) 
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which     by  Eqs.    (18)   and    (12)    is 


,Jcr,   Jfik.3  n.An  b'(-w)]   +  fa*.*''.**  b'hKx+4F*  =0  "-(24) 


dl5  U  *  D-(25) 


Using  Eq.  (17),  Eq.  (25)becomes  the  so-called  Stokes  differential 
11 


Its  solution   as  discussed  by  Budden   is  the  following: 

F^KSZ  K(S)  M27) 


where  K   is  a  constant  and  where  Ai(?)  is  the  Airy  integral  function 

defined  as  follows: 

*0 


*ifc)=  —  j  W^+  i^3;^ 


7T   ^  ^  D-(28) 

0 


Tables  of  Ai(Q  are  given  by  Miller.    A  plot  of  Ai(?)  versus  ?  is 
given  in  Figure  IV- 4.   It  may  be  seen  that  beyond  the  region  5  =  0 
(q  =  0,  x  =  -W)  the  solution  Ai(?)  is  an  exponentially  damped  wave — 
an  evanescent  wave. 
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There  is  another  solution  to  Eq.  (26)  since  it  is  a  second-order 
differential  equation.   It  is  another  Airy  integral  function  Bi(?). 
However,  beyond  the  level  x  =  -W  at  q  =  0,  the  function  Bi(?)  becomes 
indefinitely  large   as  £  increases,  i.e.,  as  x  ->  -»,   Since  it  is 
physically  unrealistic  for  F   (and  thus  E)  to  get  larger  and  larger  at 
distances  farther  and  farther  away  from  the  waveguide,  then  the  solution 
to  Eq.  (26)  cannot  contain  a  multiple  of  Bi(£). 

For  t,     sufficiently  large,  Ai(?)  can  be  given  by  its  asymptotic 

,   „,    11 

approximations 

and  , 


*  ^  ^r         ^  -f^^fT-  D- 


(30) 


where,  for  an  exponent  a, 

5    -     5  I      e  ^ 


From  Eq.  (17)  it  may  be  seen  that 

<5=  T      4-   (x-rW)  >0 


D-(31) 


<rq  $=  0       -V-    (x+wj^O 


D-(32) 


D-(33) 


where  Eqs.  (29)  and  (32)  correspond  to  propagating  within  the  guide 

2 
(where  q   >  0,  -W <  x  <  0) ,  and  where  Eqs.  (30)  and  (33)  correspond  to 

2 
the  evanescent  wave  region  (where  q  <  0,  x  <  -W) . 
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Referring  to  Figure  IV-5  and  following  the  arguments  of  Budden   and 

12  2 

Wait,   it  is  assumed  that  q   is  linear  over  the  entire  region  be. 

There  is  a  certain  range  cd_  where  |q|  is  so  small  that  it  violates  the  WKBJ 

approximation,  and  so  the  WKBJ  solutions  cannot  be  used.   Outside  this 

range  the  WKBJ  solutions  are  good  approximations.   In  the  regions  b£  and 

de  the  asymptotic  approximations  of  Ai(£)  can  be  used,  and  they  must 

be  matched  to  the  WKBJ  solutions. 

The  WKBJ  solutions  are  given  from  Eqs.  C-(52)  and  C-(62): 

It 

A 


v  f\^KJ-<"V,-!,'V-'^]  M», 


X 

r       i   f',a  r,hA\  '<■{&-&  r-     -K  J    1$'  **■ 


for  x  =  -W 


Eq .  (35)  follows  from  Fq .  (34)  using  the  following  argument.   From 
Eq .  C-(97)  for  q  (x) ,  it  may  be  seen  that,  because  n =  n  (x)sinG(x)  can 
be  set  equal  to  n  (-W)  [since  0(-W)  =  tt/2 ] ,  then  q  (x)  is  imaginary  for 
x  beyond  the  point  x  =  -W,  i.e.,  q  =  i|q|.   Since  the  E  field  cannot 
indefinitely  increase  as  x  ->  -  °° ,  there  can  be  only  one  term  in  Eq  .  (35) 
—  an  exponentially  decaying  one. 

In  order  to  show  the  relationship  between  Ai(0  and  the  WKBJ 
solutions,  Eq .  (17)  can  be  used.   From  Eqs.  (31)  and  (32) 


53/fl=  -c|9)3/a  V^S  =  t 
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Thus 


3/a 


=  -i  f  [a*>.a«  bY-wfJ[|  jtew)*«fc 

-W 

which  from  Eq .  (12)  becomes 

Eq .  (40)  can  be  expressed  as  follows: 

3  •>  I  -tf  o         J 

From  Eqs.    (17),    (31),    and    (32) 


%  =  ISlV*-^ 


,7y      i*7v 
5)      e 


=  )pJkD n0/UbMJ    (xfw)     e 


D-(37) 


D-(38) 


D-(39) 


D-(40) 


D-(41) 


D-(42) 


D-(43) 


D-(44) 


A. 


«3rt0AnbY-w)  a 


'A 


■y/fc 


'/r    tV/v 


^n0AnbT-i^w)J    e         d-(45) 
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which  from  Eq.    (12)   becomes 


*/—  n     /S    [        ko  / 

N5   "     ]         UpAnbW        C 


D-(46) 


Then  for  arg?  =  it,  substitution  of  Eqs.  (41)  and  (46)  into  Eq.  (29) 
produces 


^fe)=     ? 


-'& 


J--L  ,  _ 


3/\tA*b't~wj}  *-&/*  J      r..    ? 


•iUv** 


+ 


o  L       V 


M*i1 


^ 


"~         "" "VV  ri 


D-(47) 


The  integral  -ik  /  q  dx  is  a  constant  for  a  given  W,  and  from  Eq , 
C-(52),  can  be  written  as  ig(-W).   Letting 

*  =  J.  _L  [a*.  .A*  b'(-w)7  fc 


D-(48) 


then  Eq.  (47)  can  be  written  as 


-for-  d/i^£-7T 


For  arg?  =   0, 


'/•J        I       i  I"      ^ 


5^  )5l3/5 


D-(50) 


D-(51) 


and,    similar   to   Eq.    (40) 
> 
3/2 


3   ^  _yy 


-0 


D-C52) 
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Substitution  of  Eqs.  (50)  and  (52)  into  Eq.  (30)  produces 

2 
At  the  boundary  of  the  linear  q  region,  the  Airy  integral 

solutions  and  the  WKBJ  solutions  must  match.   Letting  x  =  -Xg  be  the 

location  of  the  boundary  between  the  regions  be  and  cd,  then 


[WKBsf]  ^      [d^] 


x=-x_  -         "-'x  =  -y. 


D-(54) 


From  Eqs.    (27),    (34),    and    (49) 


f 


From  Eq .  (55)  using  the  orthogonality  of  the  functions  e   ,  it  can  be 
concluded  that 


4.    LKe         e    D 


l#  '  D-(56) 

-rf-  ~  L  '  Te   I  D-(57) 


Tt"     ^  L  v  D-(58) 


D-(59) 


2  03 


It  should  be  stated  that  equating  the  derivative  of  the  WKBJ  and  the 
Airy  integral  solutions  at  the  level  x  =  -xQ  leads  to  the  same  results 
as  Eqs.  (58)  and  (59). 

At  the  boundary  between  the  regions  c_d  and  de,  symmetrically 
located  (with  respect  to  xQ)  on  the  other  side  of  x  =  -W,  the  WKBJ  and 
Airy  solutions  given  by  Eqs.  (35)  and  (53)  must  also  match.   Letting 
A  =   W  -  |Xq|,  and  letting  x-|_  =  -W  -  A  be  the  location  of  this  boundary,  then 


f  Ye  i 


(3  W     "  ~   'Ve,xi  n-(60) 


from  which  it  may  be  concluded  that 


fV-fS^K, 


or,  from  Eq .  (58) 


rH  <2  e  D-(62) 


At  the  n  -  n   (0)  boundary,  Eqs.  C-(78)  and  C-(83)  can  be  used  for 
a     g 

the  boundary  conditions.   They  yield  Eq .  C-(86)  which  can  be  written 
[with  the  aid  of  Eq.  B-(44)]  as 

F  -  F    e  3 

3    3  D-(63) 


TE 
where  4>   is  defined  in  Eq .  C-(92).   Dividing  Eq.  (59)  by  Eq.  (63) 

produces 


-  e  =    l  e    J      e        * 


D-(64) 
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33(-w)-hc3  <f) 


TH  IT 


D-(65) 


since  i  =  e     .   Using  Eq .  C-(53),  Eq .  (65)  becomes 

o 

24,  j  na  W  ^^)ix- <3  (^  -  x  ^^     D_(66) 

— w 

which  is  the  waveguide  equation  for  the  TE  case.   Comparison  of  Eq .  (66) 

TE 
with  Eq .  C-(96)  shows  that  the  only  change  is  to  replace  the  term  2d> 

gs 

with  tt/2  .   Thus  it  may  be  seen  that  there  is  a  phase  advance  of  tt/2 
associated  with  the  total  internal  refraction  phenomenon.   Further, 
from  Eq.  (66)  it  may  be  seen  that  the  WKBJ  solution  predicts  the  correct 
phase  change  (except  for  the  term  tt/2)  experienced  by  the  wave  within 
the  guide,  even  though  the  WKBJ  solution  itself  breaks  down  in  the 
region  q  =  0. 

Using  Eq.  C-(97)  and  C-(98) ,  Eq .  (66)  can  be  written  as 

1.    n4,  .  U  I  1   n«.  I     I     IT 

.a   /r  =  a/"T  D-(67) 


8*.  j  [rV*0-7*l  *-**-'  Mg 


Eq.  (67)  is  trancendental  in  n  ,  which  is  equivalent  to  saying  that  it 
is  transcendental  in  0(0),  according  to  Eq .  C-(101).   It  may  or  may  not 
be  transcendental  in  W  too,  depending  on  the  functional  form  of  n  (x) 
and  the  integral  in  Eq.  (67). 

Following  arguments  given  in  Appendix  C,  the  cutoff  condition  is 
given  when  either 


£■  n, 


D-(68) 
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*)   ~  ^ 


D-(69) 


For  n  >  n   (which  is  usually  the  case),  the  cutoff  condition  is  given 


by  Eq.  (68),  i.e.,  when  r=  n  .   When  n<n  ,  the  light  can  no  longer  be 

totally  internally  refracted, and  the  light  escapes  below  the  guide.   For 

this  cutoff  condition,  Eq.  (67)  becomes  an  integral  equation  for  (W  .  )   : 
/^  min  TE 

2 
where  the  (An)   term  in  the  denominator  of  the  argument  of  the  arctangent 

has  been  neglected  compared  to  2n  An. 

In  the  event  that  n  <  n   (a  much  less  common  situation)  where 
o    a 

(n  -  n  )  <  An,  then  the  cutoff  condition  is  given  when  n =  n  and  the 
a    o  a 

TF 
light  leaks  out  into  the  superstrate.   In  this  case  i>     '  =   0.   Thus 

ag 

Eq .  (67)  becomes 

o 


[  b{(*)-»Z]  Vx=  ad(w,)/s 


M      ^  D_(71) 

Since  sin  0(0)  cannot  exceed  unity,  then  from  Eq .  C-(38)  [or  Eq. 

C-(101)],  n  can  never  exceed  the  value  n  (0).   In  fact,  n  cannot  equal 

g 

n  (0)  either.   If  sin  0(0)  =  1,  then  0(0)  =  tt/2  ,  and  from  Snell's  Law 
[Eq.  C-(38)],  0(x)  =  tt/2  for  all  x.   Hence  the  integrand  in  Eq.  (66) 
would  be  zero  for  all  x,  and  in  order  to  propagate  a  light  wave,  only  an 
infinitely  thick  (W  ->  »)  guide  could  satisfy  Eq .  (66),  which  is  phys- 
ically unrealistic.   Thus  the  allowable  values  of  n  are 


n0-\-  k*  ">  i  >  n0 


D-(72) 
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in  order  that  a  wave  may  propagate  in  the  guide  (without  cutoff). 
Using  Eq.  C-(101),  this  says  that 


|  >  s^B{o)  > 


r^o-f-A"!  n_(73) 


which  gives  the  allowable  range  of  values  of  0(0)  for  solving  Eq.  (67) 

when  n  >  n  .   When  n  <  n  ,  then 
o    a        o    a 

f\n+br\  >  <n>  n* 

0  I  D-(74) 


no  +  &*  D-(75) 

The  effect  of  curvature  of  the  index  profile  must  be  investigated 
in  order  to  determine  the  validity  of  Eqs.  (12)  and  (17)  and  the  resultant 
conversion  of  the  wave  equation  to  the  Stokes  equation.   If  an  additional 
term  is  kept  in  Eq.  (6),  then 


VJ      3 


D-(76) 


which,  from  Eq.  (1)  becomes 


n  GO  =  n0t  An  [b6-w) f  b Y-w)$r+vy)  +  bT w)  (x+ W$i  J 


D-(77) 


•%  M  -  nf  f  a*.  Zln/b(-yvj+bMM+J)'^f^J    D_(78) 
3 
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and 


^  n*(-wh  i^  a  ne  A*  b(->v) 


D-(79) 


Hence 

3 


-(80) 


D-C81) 


where 


R  ~  r?DAn  b'Y-vv) 


D-(82) 


D-(83) 


Defining  £  according  to  Eq.  (17),  then 


5=-(^fB,)  J(*+>v) 


D-(84) 


and  the  wave  equation  given  by  Eq.  (20)  becomes  similar  to  Eq.  (24) 
except  for  the  additional  quadratic  term  from  Eq.  (81): 

dl-^  J  D-(85) 


9*[-5*sa(*.«,r,8,K.o  „.„ 


208 


If  Ai(?)  is  to  be  a  good  approximation  to  the  solution  of  Eq.  (86),  then 


>l»  I^CMiT*3^ 


D-(87) 


and  this  must  hold  for  values  of  £  large  enough  for  the  asymptotic 

approximation  of  Ai(f;)  (and  thus  the  WKBJ  solutions)  to  be  valid. 

According  to  Budden,    this  will  happen  when  |c|  ^  1.  So  Eq.  (87) 
requires  that,  for  |  X,  |  =  1,  the  linearity  criterion  is 


\%\«K*> 


\Ak 


D-(88) 


From  Eq.  (1) 


1^60=  An  b'U) 


D-(89) 


tf,%)  =  An  \o"(x) 


D-(90) 


hence  Eqs.  (82),  (83),  and  (88)  require  that 

,/3  ^ 


D-(91) 


Further,  from  Eqs.  (17)  and  (89),  the  condition  that  |c|  =  1  (at 

x  =  -x  and  x=  -W  -  A,  i.e.,  x=  -W±  A)  corresponds  to 
o 


D-(92) 
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It  is  now  possible  to  write  down  the  solutions  for  E   and  hL_ 

TE      TE 

using  Eqs.  B-(75)  through  B-(79)  together  with  the  results  of  this 
appendix  and  Appendix  C.   The  solutions  for  x  =  0  are  given  by  Eqs. 
C-Q32),  C-(135),  and  C-(136).   For  the  region  (-W  +  A)  =  x  =  0,  the 
solutions  are  the  WKBJ  solutions  given  by  Eqs.  C-(133),  C-(135),  and 
C-(137).   For  the  region  (-W  -  A)  =  x  =  (-W  +  A),  the  solutions  involve 
the  Airy  integral  function.   Using  Eq.  (27)  together  with  Eqs.  C-(5) 
through  C-(10),  then 


D-(93) 


H02(X)=  (i/K)^S.     ^2-  D-(95) 


D-(96) 


where  IC,  and  ?'  are  given  by  Eqs.  (58)  and  (18)  respectively,  and 

where 

00 

1$~~    ~        ^    0  D"(97) 

Values  of  d[Ai(^)]/d?  are  tabulated  by  Miller.  Finally,  for  the 
region  x  =  (-W  -  A) ,  the  solutions  are  the  WKBJ  solutions  [obtained 
from  Eq.  (35)]   and  Eqs.  C-(5)  through  C-(10): 

3  0      *  D-(98) 
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KM'    ~  (?*>)   ^g  W 


»„W  -  -  i 


Vx) 


D-(99) 


D-(100) 


where  F  is  given  by  Eq.  (62),  and  where,  in  Eq .  (100),  terms  order 
d{  Ul     J  /dx  are  neglected  according  to  Eqs.  C-(81)  and  C-(82)  and 
the  WKBJ  approximation. 

TM  Modes 


The  wave  equation  for  TM  modes  is  given  by  Eq.  A- (93)  and  can  be 


written   as 


^4         _l      X    ^[n/W 


a  j 


^kj    ~"    na6<)     ^x     "   ^. 


fif*^     D. 


(101) 


Making  the  substitution 


7,11 


U  =  Ax  In.  Oc) 


D-(102) 


then  Fq .  (101)  becomes 


<       ^x        ^x       '  l     3 


V 


D-(103) 


D-(105) 
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Thus,  re-arranging  terms,  Eq.  (105)  becomes 


f  wn,  i- 


f^X* 


u=o 


D-(106) 


Defining  an  "effective  refractive  index"   as 


NM~    3     MV  [TT    2fta«,a   <^*a 


D-(107) 


and  an  "effective  value  of  q"   as 


D-(108) 


D-(109) 


then  the  wave  equation  becomes 


D-(llO) 


From  Eqs.  (107)  -  (109) 


Q*(x)~  9aW  =  n360-  ^a(x)  D_(m) 


By  making  the  assumption  that  n  (x)  varies  linearly  near  the  level 
x  =  -W,  as  given  by  Eqs.  (6)  and  (7),  then 

D-(112) 


*X 


fW)_ 


&-0 


D-(113) 
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[<(*)]"'=  ^  (l-  ^  [k(-w)  +bW&rtj/>(U4) 


Hence,  from  Eq .  (109) 


4-   bY-w)  (>c-fW)]  f  D-ui5) 


2      2 
It  mav  be  seen  that  the  difference  between  Q  and  q   [and,  using  Eq . 

2       2  2 

(111),  between  N  and  n   ]  is  proportional  to  [An  b'(-W)/k  ]  ,  which 

is  a  very  small  quantity  for  optical  waveguides.   Hence,  near  the 

level  x  =  -W,  the  wave  equation  Eq .  (110)  becomes  approximately 

i5u    .    a,a   a. 


*  a:  <f(*)  u  -  o 


^xa     "  0  D-(ii6) 

and  all  the  techniques  described  for  the  TE  case  can  be  used  to  convert 
this  form  of  the  wave  equation  to  the  Stokes  equation,  and  to  then 
find  the  Airy  integral  solutions. 

Neglecting  the  right-hand  side  of  Eq .  (115)  is  equivialent  to  the 

2         2 
assumption  that  Q  (x)  and  q  (x)  have  a  zero  at  nearly  the  same  point. 

It  is  possible  to  calculate  this  slight  difference  between  the  TE  and 

TM  reflection  points  (i.e.,  the  points  of  total  internal  refraction). 

It  is  convenient  to  define  the  TM  reflection  point  as 


2 
where  V   is  the  value  of  x  at  which  q   is  zero,  and  where  6W  is  the 

difference.   Using  Eq .  (12),  Eq .  (115)  can  then  be  written  as 
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ft*M--a  +  U+Kfr^)oLi  D-(n8) 


rM     °    /  ^  | 


where  the  constants  0n  and  0,  are  defined  as  follows 


<*,*{>".">■<■•>> -ft  Mr*?] 


D-(119) 


D-(120) 


In  order  that  0(-W  )  =  0,  then  from  Eq.  (118)  it  is  necessary  that 


sw 


_  _a 


d,  D-(121) 


£  -  |  J^y).  „.(122) 

O     6 


It  may  be  seen  that  SW  is  indeed  small  for  optical  waveguides,  and 
Fq.  (116)  is  valid  for  TM  modes  to  a  very  good  approximation.   Further, 
since  6W  <  0,  then  VI       <   W,  i.e.,  a  TM  wave  reflects  at  a  value  of 
|x|  which  is  smaller  than  in  the  TE  case. 

In  Fq .  (114),  only  the  linear  term  was  retained.   If  instead  a 
quadratic  term  is  also  retained,  then 

+  (x+wf/o[6*  bY-w)/n0]    }  D_(123) 

where  it  is  assumed  that  terms  containing  the  factor  b"(-W)  are  small 
enough  according  to  the  linearity  criterion  given  by  Eq.  (88).   Eq.  (115) 
[or,  equivalently,  Fq .  (ID0*)]  then  becomes 
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+  (*4w)(ane  am  b'l-vd  +  ia  [An  b' (-*)]' 

1  w M 


D-(124) 


qY*)=-&0  +  L*+%c^&i-k+\cMfQLA      D-d25) 


D-(126) 


where  the  constant  Q   is  defined  as 

a^  17  t— «: — J 


Collecting  terms  in  Eq .  (125)  leads  to 

Using  Eq.  (121)  together  vzith  the  fact  that  6W  is  small  and  Q   <<  Q 
then  Eq .  (127)  becomes  to  a  good  approximation 


ast*) *  (w-O a,  -  (*+ wr«)=  aa  d_(i28) 


2 

which  is  of  the  same  form  as  Eq .  (81)  for  q  (x) .   Now,  it  is  required 

2 
that  0  (x)  be  linear  near  the  region  x  =  -W  [or,  nearly  the  same, at  the 

2 
level  x  =  -W   ] ,  just  as  it  was  required  that  q   be  linear  for  the  TE 

case.   The  quadratic  term  in  Eq .  (128)  is  subject  to  the  linearity 
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criterion  given  by  Eq .  (88),  i.e., 
Thus,  from  Eqs.  (120)  and  (126) 


D-(129) 


D-(130) 


where  the  second  term  of  Q,  i-n  £<!•  (120)  has  been  neglected  in  Eq. 
(130).   Re-writing  Eq .  (130)  leads  to 


n0  ^>lpgl  a„  b'(~*)/h  D-(i3D 

Since  the  minimum  value  of  n  is  n  before  cutoff,  then  Eq .  (131) 

o 

requires  that 

2 


,f  V>  ^f   An   h' (-*)/*< 


D-(132) 


Budden   defines  the  quantity  B  as  follows: 


[where  optical  waveguide  notation  has  been  used  for  the  right-hand  side 
of  Eq .  (133)  instead  of  Budden's  radio  ionosphere  notation].   Using 
Eq.  (8),  Eq.  (133)  becomes 


S  -  >f  L  <V<' -'  '-  °"  °[""^  D-C134) 


=  [JkJs».i»  bU)f 
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Again  setting  n  =  n,  then  Eq.  (134)  becomes 


=r  B       \Z   bn  b't-ri/Aj 


D-(135) 


Using  Fq.  (135),  Eq.  (132)  becomes 

4\ 


D-(136) 


B>*>   I.  '7 


D-(137) 


For  optical  waveguides,  Eq.  (134)  almost  always  satisfies  Fq .  (137). 
That  is,  if  the  linearity  criterion  given  by  Eq.  (88)  is  satisfied  for 
TF  modes,  then  it  is  satisfied  for  TM  modes  also,  to  an  extremely  good 
approximation.   For  instance,  if  the  index  profile  is  exponential  such 
that 

then  Eq.  (138)  satisfies  Eqs.  (1)  -  (3).   If  the  following  values  are 
chosen:  An  =  0.1,  D  =  lym  (which  is  a  fairly  steep  profile),  Aq  = 
0.63ym,  d  =  2/D  =  2ym    (so  that  the  e~  point  of  the  index  profile 
occurs  when  x  =  -D/2  =  -0.5um);  also, letting  n  =  uq   =  1.5  (i.e. ,  the 
wave  is  propagating  nearly  at  cutoff),  then  b'(-W)  =  b '  (-D)  = 
-d/(l  -  edD)  =  0.313;  for  these  values,  then, from  Eq.  (134),  B  =  51, 
which  certainly  satisfies  Eq.  (137).   For  waveguide  index  profiles  with 
An  =  0.01  over  a  distance  D  =  lOum,  B  has  the  value  of  1089. 
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Returning  now  to  Eq.  (125),  the  solution  U  has  the  form  given 
by  Eq.  (27).   Hence,  from  Eq.  (102) 


This  has  to  be  matched  to  the  WKBJ  solutions  at  the  levels  x  =  -W  ±  A, 
as  in  the  TE  case.   [Eq.  (92)  defines  the  value  of  A.]   The  WKBJ 
solutions  are  given  from  Eqs.  C-(52)  and  C-(167):        ^ 


for  -W  =  x  =  0 


i  \-'h/  u.^^fi    -•*•  i ,  ft'1** 


D-(140) 


D-(141) 


for  x  =  -W 

Eqs.  (139)  -  (141)  differ  from  Eqs.  (27),  (34),  and  (35)  only  by  the 
term  n  (x) .   Hence,  upon  matching  the  boundary  conditions  at  x  =  -W  ±  A 
[and  using  the  relations  given  by  Eqs.  (58),  (59),  and  (62)], 

KrM-   (Aa/Fa)  Krc- 

A3=  (^/Fj  F3  D_(143) 


D-(144) 


«r(rtA)F, 


At  the  n   -  n  (0)  boundary 
a    g 

"3  ~  As  e  D-(145) 
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TM 
similar  to  Eq.  (63).   [Fq.  C-(194)  defines  the  quantity  <j>   J.   Thus  the 

waveguide  equation  for  the  TM  case  differs  only  slightly  from  the  wave- 
guide equation  for  the  TE  case,  i.e., 

0 


D-(146) 


similar  to  Eq.  (66).   Using  Eqs.  C-(97)  and  C~(194),  Eq.  (146)  can  be 
written  as 


similar  to  Eq.  (67).   The  cutoff  condition  for  n   >  n   is  given  by 
Eq .  (68),  in  which  case  an  integral  equation  for  (W  .  )™w  results: 


— /w  ■  ^  ..  air  t 

4-  W  / 


a\'/. 


^      I  5nflAn 


similar  to  Eq .  (70).   In  the  event  that  n   <  n  ,  then 
^  o    a 


D-(148) 


D-(149) 


which  is  identical  to  Eq.  (71).   Again,  the  allowable  values  of  n  and 
0(0)  are  given  by  Eqs.  (72)  -  (75). 

It  is  now  possible  to  write  down  the  TM  solutions  using  Eqs.  B-(149) 
through  B-(153),  Appendix  C,  and  this  appendix.   For  the  region  x  =  0,  the 
solutions  are  given  by  Eqs.  C-(221),  C-(224),  and  C-(227).   For  the  region 
(~W  +  A)  =  x  =  0,  the  solutions  are  the  WKBJ  solutions  given  by  Eqs.  C-(222), 
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C-(225),  and  C-(228).   For  the  region  (-W  -A)  =  x  =  (-W  +A) ,  the  solutions 
involve  the  Airy  integral  function  [using  Eqs.  C-(142)  through  C-(147)]: 


e« 


3   T   "=      '  D-(150) 

(x)=  j^/A.n^x)]    rt    frj  N1S1) 


^M-I'/^f^^^i,]^. 


For  the  region  x  =  (-W  -  A) ,  the  solutions  are  the  WKBJ  solutions 
obtained  from  Eq.  (141):  X 


(152) 


D-(153) 


D-Q54) 


where  in  Eq.  (155),  terms  of  order  d  |  n  |q|     r  /d  are  neglected 
according  to  Eqs.  C-(181)  and  C-(182)  and  the  WKBJ  approximation. 


APPENDIX  E 
THE  PRINCIPAL  AXES  AND  INDEXES  OF  CUBIC  ZINC  SULFIDE 

The  problem  is  to  find  the  angle  fi  through  which  the  y  and  z  axes 

a       A 
must  be  rotated  to  the  new  axes  y'  and  z1,  since,  from  Eqs.  VI-(5)  and 

A       /V  A        * 

VI-(6),  x  =  x  and  n  =  n  .   Referring  to  Figure  E-l,  y  and  z  can  be 


written  as 


3  E-(l) 


£  -  £  '  d^o-TL  -  y  '  &«*»  -A- 


E-(2) 


Substituting  Fqs.  (1)  and  (2)  into  Eq .  VI- (3)  results  in 

'-2    ~         '  /   -  a  *       i 

E-(3) 


2 
Subtracting  Eq.  VI- (4)  from  Eq.  (3),  and  using  the  fact  that  cos  fi  + 

2 
sin  £2  =  1, 
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X,  X 


Figure  E-l .   Angular  relationships  between  the  (x,  y,  z)  set  of  axes 

A     A     A 

and  the  (x  ,  y  ,  z  )  set. 
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3 

*  2 

-fy'2'2  r^  e,  (  c^s,-0.-s^2I2_)=  O  E_(4) 


The  angle  Q   is  that  angle  which  causes  the  y'  z'  terms  to  vanish,  i.e., 


uV  5  rfj  ex  (ce^^-s^ta.)  ~° 


2       2 
For  arbitrary  e  ,  this  requires  that  cos  Q  =  sin  ft;  thus 


E-(5) 


-TL  =rr/i  +  j  V5  ,     J 


=  0  /  2  ••• 

J    J      '  E-(6) 


It  is  convenient  to  use  j  =  0  in  Eq.(6)  with  no  loss  of  generality.   The 

new  axes  x' ,  y',  and  z'  —  the  principal  axes  —  are  shown  in  Figure 

VI-1. 

A       A  ,2.2 

Because  the  y'  and  z'  axes  are  orthogonal,  the  y   and  z   terms 

in  Eq.  (4)  must  each  be  zero.   Using  Q  =  it/4,  this  requires  that 
/a/  _i_  _ 


j      I  "*         n5a    ^'* 


-^  **  -0 


E-(7) 


~r:  eK)=0 


Z     I    n;-'    "    n2,a  '   r4ic*y~"u  E-(g) 


Solving  Eqs.  (7)  and  (8)  then 
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r,  e. 


n(3   ~        fl   3       +  rv,  CX  E-(9) 

_JL_     _    — L     __  r    p 

n  <  =>    ~       yn   3  4<       *  E-(10) 


■2 


O 


Letting 


n^  s  n0  +A*' 


n2  s  n0  +  Aoz 

then  Eqs.  (9)  and  (10)  can  be  written  as 
I  \ 


E-(ll) 


E-(12) 


'  -^  -  r;  e,  E_(14) 


n. 


•A^)A 


a 


In  general,  the  index  changes  An  and  An  are  very  much  less  than 

y       z 

n   such  that  An' /n   <<  1  and  An' /n   <<  1.   So  Eqs.  (13)  and  (14)  can  be 
o  y  o  z   o 


E-(15) 


written  as 

£•(' 

- 

M 

l+    "*rHteKj 

(        /l 

- 

I 

(l-    ^\^ 

fl*  (    ' 

^ 

E-(16) 


Thus 

3 


3  E-(17) 


'       3      / 

An2-   -f  h0  >;,  ex/3  e-(is) 


Substituting  Eqs .  (17)  and  (18)  into  Eq.  (11)  and  (12)  [and  repeating 
Eq.  VI-(6)],  the  principal  indexes  are 


V\J  =  KQ 


nM'  ^  n0-  n0    rH[  €x/b 


E-(19) 


E-(20) 


■Z  °  o         <f/  r  E-(21) 
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